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Crystal Pattern Synthesis by an Approximate Summation of Fourier Series 


JoserH S. LUKESH 
Crystallographic Laboratory, Department of Geology, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received February 10, 1947) 


A faithful representation of the locations of the peaks of electron density maps can be 
obtained rapidly by determining for each point in the projected unit cell the algebraic sum of 
the amplitudes of those reflections which are contributing either a maximum or a minimum 
at each point. Since maxima and minima do not, in general, fall at the specific points of the 
grid, an approximation must be introduced. Using the semi-polar form of Fourier series, the 
approximate summation can be accomplished with the use of one hundred and fifty strips. 
The method was tested by determining the peak locations in the structure of diopside projected 
on 001. Potential applications of the approximate summation include rapid checking of the 


validity of assumed crystal structures. 


INTRODUCTION 


RYSTAL patterns can be synthesized con- 

veniently by determining the projected 
electron density through the summation of a 
Fourier series of the following general type: 


p(x, y) = (1 ‘A)d en Fu 0) cos(27r N)(hx+ky), 
(1) 


which gives the projection normal to the ¢ axis. 
Other axes are treated by appropriate choice of 
the indices, h, k, or 1. Equation (1) is valid only 
when the projection is centrosymmetrical. The 
non-centrosymmetrical case will not be con- 
sidered here; however, it can be treated in a 
similar manner by applying the same reasoning 
to the sine terms. The form of the contribution 
of a given reflection, hk0, to the total electron 
density is invariant, regardless of the actual 
crystal structure. It consists of a 
alternation of positive and negative density. 
Actually it is a cosine surface whose periodicity 
in-each axial direction is given by the appropriate 
index. The different distributions of electron 
density in different crystals arise from varying 


sinusoidal 


amplitudes of the cosine surfaces and from 
different arrays of contributing reflections. Ob- 
viously, in order that an accumulation of electron 
density occur at any point in the projected unit 
cell, the sum of the positive solutions of the terms 
of Eq. (1) must be greater than the sum of the 
negative solutions. In order that a high peak 
occur, not only must a large number of reflections 
give positive solutions, but also they must be 
near their maxima. Similarly, in order that 
there be little or no accumulation of electron 
density, either most reflections are contributing 
very little, or those near their maxima are 
compensated for by others near their minima. 
This suggests that if one were to list for each 
point in the unit cell the amplitudes of all 
reflections that are contributing a maximum and 
all that are near a minimum, one might obtain a 
faithful representation of the electron density 
distribution merely by subtracting the negative 
terms from the positive. If this hypothesis is 
valid, the labor of computation is reduced to the 
listing and totalling of the appropriate ampli- 
tudes. Of course, detail will be lost in the regions 
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of low electron density, but these regions are 
unimportant in at least the early stages of 
crystal structure analysis. In many cases, loca- 
tion of the strongest peaks will be of material 
aid in the solution of a crystal structure. 


THEORY 


In an earlier paper the writer' developed a 
semi-polar form of the series of Eq. (1). This 
concept, while not essential to the approximate 
method to be discussed, offers a number of 
advantages in computational simplicity and will 
be followed in this report. If the point, x, y, is 
specified by the coordinate, x, and the slope, c, 
of the line from the origin to the point, Eq. (1) 
may be written: 


p(x, c)=(1/A) don Doe Fonco) cos(24/N)x(h+ck). 
(2) 


The half-unit cell for which the summation is to 
be made is divided into four segments, and the 
form of the cosine term varies for each. For a 
given plane and a given radial line, Eq. (2) 
reduces to 


p(x, c)=(1/A) don Doe Fone) cos(24/N)Cx, (3) 


where C is a constant and x is an integer. In the 
article cited' it was pointed out that the con- 
ventional rectilinear grid is inconvenient, since 
the slopes of the lines from the origin to the 
various points are not suited for easy calculation. 
Furthermore, only a few points occur along any 
one line. For these reasons, a semi-polar grid 
was described in which the points at which Eq. 
(3) is to be evaluated are defined by the inter- 
sections of radial lines and either the x or y 
coordinates. If the axes are divided into thirty 
points (a sufficient number since the projection 
is only an approximation), the number, JN, in 
Eq. (3) is thirty. For this spacing of points, it is 
convenient to choose radii whose slopes are .1, .2, 
.3, etc., and the constant, C, can be determined 
from the indices of a plane merely by inspection. 
Thus, for the reflection, 320, the constant along 
the radius whose slope is .4 is equal to 3 plus 
4(2), or 3.8. 

If Eq. (3) is plotted as a function of x for a 
given amplitude, F, and a given constant, C, a 


' Joseph S. Lukesh, J. App. Phys. 18, 321 (1947). 
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cosine curve results. This is the form of the 
density contribution of a reflection, hkO, along a 
radial line of slope, c. The amplitude is deter- 
mined by F and the periodicity by C. It is 
obvious that maxima and minima occur at 
various points along the radius. Since the posi- 
tions of the maxima and minima depend only 
on the indices of the reflection and on the slope 
of the radius, they can be determined in advance 
and tabulated for use in all structure determina- 
tions. Because, as was pointed out in the article 
previously cited,' there are only one hundred and 
fifty independent values of C for the grid chosen, 
the table will not be extensive. 

In general, the maxima and minima will not 
occur at integral values of the coordinate, x. 
This necessitates the introduction of an approxi- 
mation, since it is at the integral values of x that 
the summation is being made. For the purposes 
of this report, it has been assumed, quite arbi- 
trarily, that a plane contributes to the electron 
density only at those points where its amplitude 
times the appropriate value of the cosine is fifty 
or more percent of the maximum. Similarly, a 
plane subtracts from the density only when the 
product is fifty or more percent of the minimum. 
In other words, a reflection is assumed to be 
ineffectual when the product of its amplitude 
times the cosine is between zero and either plus 
or minus one-half of F. Furthermore, it is as- 
sumed that a plane contributes or subtracts the 
entire value of F when it is effective. This latter 
assumption eliminates all calculation, and the 
summation consists of tabulation of the ampli- 
tudes of the effective planes at each point and 
their summation. 


EXAMPLE OF THE APPROXIMATION 


In order to test the hypothesis developed 
above, computations have been made of the 
electron density of diopside projected on (001). 
Data for the amplitudes of the various reflections 
and a map of the conventional summation have 
been reported by Bragg.? Because of the high 
symmetry of diopside in this projection, only 
one-eighth of the cell need be considered. How- 
ever, in order to test the approximation, compu- 
tations were made for one-quarter of the cell. 


2 W. L. Bragg, Zeits. f. Krist. 70, 475 (1929). 
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TABLE I. Coordinates, x, for which the function: cos(#/15)Cx has maxima and minima. 
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Fic. 1. Peak locations in the 001 projection of diopside 
determined by the approximate summation, Case I 
(equal amplitudes). 


Since the quarter-cell is in itself centrosym- 
metrical, the amount the major peaks depart 
from centrosymmetry will be a test of the value 
of the approximation. Three electron density 
maps have been computed. In the first, only the 
signs of the amplitudes have been used, all 
magnitudes being considered equal. In the 
second, estimated magnitudes have been intro- 
duced; and in the third, the actual measured 
values were used. 

In Table I are listed the one hundred and 
fifty values of the constant, C, of Eq. (3) and the 
corresponding values of x for which maxima and 
minima, as defined previously, occur. To com- 
pute the summation it is only necessary to list 
for each radial line the amplitudes of the reflec- 
tions which are maximum and minimum at each 
point. As a typical example, the maxima and 
minima of the plane 570 along the radial line 
whose slope, c, is .4, will be determined. The 
constant, C, is equal to 5 plus .4(7), or 7.8. 
From Table I, opposite 7.8 it is found that 
maxima and minima occur at the following 
points: 


Maxima—0, 4, 8, 11, 12, 15 
Minima— 2, 6, 9, 10, 13, 14 


If the amplitude of the reflection, 570, is plus 
forty, the number forty is recorded in the col- 
umns of the worksheet corresponding to the 
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Fic. 2. Peak locations in the 001 projection of diopside 


determined 


the approximate summation, Case II 


(estimated amplitudes). 


coordinates, x, for which the contribution is a 
maximum. Similarly, minus forty is recorded in 
the columns corresponding to the minima. When 
the amplitude is negative, one must, of course, 
make appropriate sign changes. 


Case I—Equal Amplitudes 


This example, which represents an exploration 
of the possibility of determining atomic locations 
from knowledge of the array of reciprocal lattice 
points, but not their relative weights, is of little 


TABLE II. F values of the reflections of the type hk0.* 


hkO F hkO F 

200 0 040 0 
400 15 240 16 
600 100 440 — 50 
800 76 640 29 
110 0 150 — 84 
310 76 350 —55 
510 — 63 550 0 
710 41 750 —96 
020 19 060 —94 
220 —5i1 260 31 
420 —23 460 0 
620 —19 170 0 
820 —17 370 0 
130 13 570 40 
330 50 080 0 
530 37 280 — 32 


* After W. L. Bragg, Zeits. f. Krist. 70,475 (1929). 
Note: Because of symmetry, reflections of the types Ak0, hkO, and 
* hkO have the same F values as those of the type Ak0. 
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practical importance, since one can always esti- 
mate the relative weights from a photographic 
film. It is included here because the results are 
rather striking. 

The data used in this and succeeding cases 
were taken from Bragg? and are listed in Table I]. 
The peaks of the “electron density’’ map drawn 
from the computations are represented in Fig. 1 
by circles. (The actual map is, of course, dis- 
torted. In the interests of clarity, only the peak 
positions are shown.) Figure 1 should be com- 
pared with Fig. 4, which shows the peak locations 
taken from the conventional summation of 
Bragg.” It is important to note that the major 
peaks, representing heavy atoms or two or more 
atoms superimposed, are correctly positioned ; 
small peaks are slightly displaced. 


Case II—Estimated Amplitudes 


In this example, the amplitudes listed in 
Table II have been assigned weights of from 
one to five. Those from one to twenty-five, for 
instance, have been given the weight of one; 
those from twenty-six to fifty, two, etc. This 
was done in an attempt to simulate amplitudes 
that might be calculated from intensities esti- 
mated from an x-ray film. Figure 2 shows the 
locations of the peaks in the summation. Im- 
provement is noted in the positions of the 
smaller peaks. 























l@ 
WZ 
a 
Fic. 3. Peak locations in the 001 projection of diopside 


determined by the approximate summation, Case III 
(measured amplitudes). 
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Case III—Measured Amplitudes 


In Fig. 3 are shown the peaks, with their 
relative heights indicated by contours, resulting 
from the approximate summation in which 
measured amplitudes were used. There is no 
appreciable change in the positions found in the 
previous case. 


DISCUSSION 


The parameters of the peaks shown in Figs. 1 
through 4 are listed in Table III. In all cases, 
the values obtained from the approximate sum- 
mations are within .02 of those from the con- 
ventional summation. (The values for the con- 
ventional summation may or may not be in 
agreement with the actual parameters of diop- 
side. Those listed were measured from the elec- 
tron density map in order to present a fair 
comparison with those derived in the same 
manner from the approximate summations.) On 
the basis of this one test, it would appear that 
the electron density peaks can, indeed, be located 
with reasonable accuracy by considering only 
the maximum and minimum contributions at 
each point in the cell. One would, however, 
hesitate to recommend the approximate method 
until independent checks have been made. The 
choice of diopside may have been a fortuitous 
one. 




















Fic. 4. Peak locations in the 001 projection of diopside 
determined by the conventional summation (after Bragg). 
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TaB_e III. Parameters measured from Figs. 1 through 4. 
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10 
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ae 
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40 


One of the potential uses of the approximate 
method is in the checking of an assumed struc- 
ture. If a logical structure can be deduced from 
crystal chemical knowledge, the phases of the 
strong reflections can be computed. An approxi- 
mate summation, using only those reflections, 
might suffice to demonstrate whether or not the 
assumed structure is consistent. 


COMPUTATION BY MEANS OF STRIPS 


Although the labor of computation is reduced 
in the approximate method to the addition of 
columns of numbers, it is still necessary to look 
up data in a table and to record them. This 
effort can be eliminated by preparing a series of 
strips, one for each value of the constant, C. 
Each strip has a space for each value of the 
coordinate, x, and those spaces where the func- 
tion is a maximum are colored blue. Minima are 
indicated by red. The labor of summation in- 
volves only the selection of the strips, their 
arrangement on a suitable board, and the addi- 
tion of the coefficients. 

In checking an assumed structure, it is not 
necessary to perform the entire summation. Only 
those regions which should be high in density 
need be treated. Using the approximate method, 
the summation can be done in a very short time. 
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Single Crystal Electron Diffraction by Micro-Crystalline Materials* 


NorRMAN DavipsoNn** AND JAMES HILLIER 
RCA Laboratories, Princeton, New Jersey 


(Received March 11, 1947) 


It is shown that the double lens probe-type electron diffraction camera described earlier is 
capable of producing useful electron diffraction patterns of single microscopic crystals in the 
size range 200-2000A; the formation of an opaque carbonaceous deposit, which occurs when 
an electron probe is focused on a specimen, is prevented in the diffraction camera by previously 
exposing the specimen chamber to an electron spray. A number of examples of two-dimensional 
grating patterns and Kikuchi-line patterns obtained from such crystals are shown. From 
either type of pattern it is possible to obtain accurate values for some of the lattice spacings, 
while from the Kikuchi-line patterns it is possible to deduce accurately the orientation of the 
crystals relative to the electron beam. The results are shown to agree reasonably well with the 
predictions of the dynamical theory as presented by MacGillavry. It is concluded that a 
two-dimensional grating pattern is to be expected from a crystal oriented with an important 
lattice vector approximately parallel to the electron beam independent of slight disorders or 
warping of the crystal. Kikuchi-line patterns are to be expected for most other crystals which 
are not excessively thick or disordered. The theory also shows that both types of patterns 





may be obtained simultaneously. 





I. INTRODUCTION 


HE study of the diffraction of fast electrons 

by single crystals has been pursued in 
several different ways up to date. Because of 
the limited penetration of the electron beam, 
work with macroscopic crystals which are readily 
manipulated and oriented is restricted either to 
reflection diffraction by a surface or, for trans- 
mission electron diffraction, to the relatively 
rare case of materials, such as mica, which can 
be obtained as thin sheets. 

There would be many useful applications of 
an effective technique for obtaining single crystal 
electron diffraction patterns and thus deter- 
mining the lattice constants of crystals of di- 
mensions less than a micron—that is, crystals of 
dimensions appropriate to electron diffraction 
and electron microscopy. Some of these applica- 
tions are: 

(a) Identification of a material which has a complex 
powder pattern or which is present as a minor constituent 
of a mixture of crystals. 

(b) The determination of crystal structure when single 
crystals of sufficient size for x-ray diffraction are not 
available. 

* Presented before the American Society for X-Ray and 
Electron Diffraction, June 10, 1946, at Silver Bay, New 
York. 

** Present Address: Gates and Crellin Laboratories of 


Chemistry, California Institute of Technology, Pasadena, 
California. 
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(c) Further study of the fundamental phenomena of 
electron diffraction by crystals. 

Various methods of observing the electron 
diffraction patterns from single micro-crystals 
have been reported in the literature. The ob- 
servations of crystalline reflections in electron 
micrographs? have been used to confirm the 
dynamical theory of electron diffraction, to study 
the warping of crystals, and in general to provide 
information about strong reflections by properly 
oriented crystals. The method is limited to 
strong reflections. 

Boersch has obtained electron diffraction pat- 
terns in which the collimated electron beam was 
restricted by an aperture to have a cross section 
of 0.005 mm at the specimen.’ The spot patterns 
obtained for the various films examined were 
such that only one or a few crystals were re- 
sponsible for the diffraction, but no correlation 
was possible between the morphology of the 
irradiated portion of the specimen and _ its 
diffraction pattern. 

The same author has suggested a scheme in 
which a fine,aperture moves in the image plane 
of an electron microscope.‘ When the aperture 





1R. D. Heidenreich and L. Sturkey, J. App. Phys. 16, 
97 (1945). 

2 J. Hillier and R. F. Baker, Phys. Rev. 61, 722 (1942). 
References 1 and 2 refer to other papers on this topic. 
3H. Boersch, Zeits. f. Physik 116, 469 (1940). 

*H. Boersch, Ann. d. Physik 27, 78 (1936). 
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is placed on the image of a particular crystal, 
the angular distribution of the transmitted rays 
provides information about the diffraction pat- 
tern of the particular crystal. Boersch used this 
technique to obtain diffraction patterns of 0.08- 
mm portions of gold foil. This method would be 
difficult to develop into a practical method for 
true micro-crystals because of the effects of 
spherical aberration on the diffracted rays, and 
of the decrease, caused by magnification, in 
angular aperture of an imaging pencil. 

Ardenne has suggested several schemes for 
modifying the electron optical system of an 
electron microscope so as to obtain diffraction 
patterns from fields 1-10 microns in diameter.® ® 

The purpose of the present article is to report 
the phenomena observed using a mode of opera- 
tion of the electron diffraction camera described 
by Hillier and Baker’ wherein a beam of electrons 
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Fic. 1. Schematic diagram of lens system and ray paths 
in the electron diffraction camera as used for shadow 
microscopy and single crystal diffraction. 


5M. von Ardenne, Zeits. f. Physik 117, 515 (1941). 

® M. von Ardenne, E. Schiebold and F. Gunther, Zeits. 
f. Physik 119, 352 (1942). 

7 J. Hillier and R. F. Baker, J. App. Phys. 17, 12 (1946). 
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reduced to a probe of ca. 200A cross section, and 
of maximum angular aperture ca. 7.51073 
radian is diffracted by single micro-crystals. It 
is believed that this method offers a practical 
approach to the problem of obtaining electron 
diffraction patterns from single micro-crystals, 
and of correlating the morphology and orienta- 
tion of a particular crystal with its diffraction 
pattern. 


Il. EXPERIMENTAL TECHNIQUES 


The double lens system of the electron diffrac- 
tion camera can be used in the two following 
ways (Fig. 1): 

(a) Using the first lens as a strong lens to 
form a reduced image, V, of the electron source, 
and the second lens to re-focus the image on the 
photographic plate or fluorescent screen, a 
focused electron diffraction pattern is obtained 
of an area of specimen determined by the size of 
the limiting aperture, A, of the second lens. 

(b) When the focal length of the second lens 
is reduced, the image V of the source is re-focused 
to an electron probe P, somewhere between A 
and the image screen. A bright field shadow 
image of the specimen is formed on the image 
screen within the central cone of illumination, C. 
The rays diffracted out of C by crystalline 
material in the irradiated field are the subject 
matter of this article. 

From the point of view of obtaining a con- 
ventional powder diffraction pattern, one is 
producing, in this mode of operation, a ‘‘de- 
focused”’ electron diffraction pattern. The degree 
of defocusing may be expressed by the numerical 
value of the angular aperture of the convergent 
beam of electrons passing through crystalline 
material in the field of irradiation. It will be 
seen, however, that under some circumstances 
it is possible to obtain precise values of inter- 
planar spacings and lattice constants from de- 
focused patterns. 

For the sake of clarity, some of the well-known 
facts about the operation of a shadow micro- 
scope** should be mentioned before describing 


8 These factors are discussed in some detail and other 
references given in reference 7. 

®V. K. Zworykin, G. A. Morton, E. G. Ramberg, J. 
Hillier, and A. W. Vance, Electron Optics and the Electron 
Microscope (John Wiley and Sons, New York, 1945), 
Chapters III and XIX. 
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Fic. 2. Two-dimensional grating pattern from graphite, 
kho, reproduced from Hillier and Baker, reference 7, Fig. 
15 (c). The bright field shadow image of the crystal has 
been printed on from a negative given a shorter exposure. 
Note the variation of reflecting power of different areas 
of the crystal for different reflections. 


the diffraction results. The magnification and 
the resolving power of a shadow microscope are 
dependent on the probe position with respect to 
specimen and electron-recording screen and the 
probe size. To obtain a small probe, the focal 
length of the lens should be as short as possible. 
To obtain high magnification the probe should 
be reasonably close to the specimen. In practice 
a probe of minimum cross section of ca. 100 

200A is obtainable with lenses of 2-3 mm focal 
length. While no precise tests of the resolving 
power of the present instrument have been 
made, resolutions less than 1000A and probably 
somewhat more than 200A have been obtained 
with magnifications of from 4X 10* to 10‘. With 
the pole piece configuration employed the focal 
point of the second lens for a probe of ca. 200A 
cross section was beyond the lens opening; it 
was, therefore, easy to mount the specimen close 
to the probe position. Under these conditions 
and with a 25-» aperture for A, the maximum 
angular aperture of irradiation of the specimen 
was ca. 7.5X10-* radian. The angular aperture 
decreases rapidly as the separation of the probe 
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Fic. 3. Two-dimensional grating pattern, BeO (hexa- 
gonal), hol, dioo=2.33A, door =4.38+0.02A. (ASTM, 
2.34A, 4.38A, respectively.) The bright field shadow image 
of the crystal has been printed on from a negative given 
shorter exposure. The 104 reflection is_relatively much 
brighter than for the x-ray case; the 301 and 303 are 
evident, although the calculated structure factor for these 
reflections is zero. 





Fie. 4. BeO, hhl, dii0= 1.33;+0.01A, doo2 = 2.20;+0.01A. 
(ASTM, 1.350A, 2.19A.) Note the intensity of the hh4 
and hh4 reflections. The central image of the crystal has 
been printed on from a shorter exposure. 


and specimen is increased, and was 10-* for 
most of the work described below. 

When the probe is focused directly on the 
specimen, the magnification is infinite and no 
specimen image is observed, but rather a pro- 
jection of the limiting aperture, A, provided the 
specimen is partially transparent. Actually, be- 
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cause of asymmetries in the lens,!° the image will 
first go to infinity in one direction, and then in 
another perpendicular to the first as the focal 
power is varied. Images observed at near infinite 
magnification were somewhat distorted because 
of this asymmetry. 


The correlation of diffraction pattern with 
‘morphology and orientation of the crystal can 
be made on the basis of the electron shadow 
micrograph. Alternately, the shadow micrograph 
may be used for control purposes. With a little 
practice, the same crystal can be located in the 
field of view of a conventional magnetic electron 
microscope, which provides a higher resolution 
image of the crystal. One difficulty with the 
electron shadow microscope as used here is that 


J. Hillier, J. App. Phys. 17, 307 (1946). 
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Fic. 5. Tip of aAl,O;-H.O crystal (béhmite) (ortho- 
rhombic) 2h02/, dioo=1.86+0.01A, doo2=1.42+0.01A. 
(ASTM, 1.86A, 1.43A.) (a) Short exposure showing image 
of crystal and one strong diffracted image, 202. (b) Longer 
exposure showing full two-dimensional pattern and some 
2h+1, ?, 21+-1 reflections. (c) Probe focused on specimen ; 
each diffraction spot is an image of the defining aperture, 
some containing the striations discussed in the text. The 
Kikuchi-line pattern shows that the strongly diffracting 
202 plane is oriented at its Bragg angle with respect to 
the incident beam. 


the magnification, which varies with the exciting 
current of the magnetic electron lens, the position 
of the pole-piece material on its hysteresis loop, 
and the specimen-lens distance, was unknown 
for any particular photograph. It is estimated 
that most of the crystals shown here for which 
conventional micrographs are not shown were 
ca. 2000A—5000A in dimensions. The aAl,O;H,O 
crystals are usually 5—10 microns long. 


III. EXPERIMENTAL RESULTS 


The diffraction results are very similar to 
those observed by Kossel and Méllenstedt for 
the diffraction of a convergent beam of electrons 
by thin sheets of mica.!! The present work is an 
extension of their results to single crystal 
analyses of arbitrary micro-crystals. 

Two principal types of diffraction patterns 
were observed under the experimental conditions 
described above: (a) two-dimensional grating 
patterns and (b) Kikuchi-line patterns. Figures 
2-7 are illustrations of two-dimensional grating 


11 W. Kossel and G. Mdllenstedt, Ann. d. Physik 36, 113 
(1939). 
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(a) 


Fic. 6. Diffraction pattern from a pair of AlxO;-H2O 
crystals. In the short exposure (a), there is one strong 
diffracted image (002) from one of the crystals. One 
portion of this crystal diffracts much more strongly than 
the other. It is evident from the position of the Kikuchi 
line in the longer exposure (b) that the strongly diffracting 
region of the crystal makes the Bragg angle with the 
incident beam. A two-dimensional grating pattern is 
developing in (b) and became much more evident in 
longer exposures not shown here. Because the two crystals 
are not exactly parallel, the diffracted images are separated. 
Note the variation in reflecting power of neighboring parts 
of the crystal for a particular reflection, especially in the 
002 and 202 reflections. This was noted in the other images 
on longer exposure. An electron micrograph of the same 
pair of crystals is shown, for comparison, in (c). Because 
of the large field of view required to include both crystals, 
the shadow micrograph in (a) was not obtained with the 
highest usable magnification and the smallest probe size, 
and is not intended to illustrate the ultimate in performance 
of the shadow microscope. 
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patterns. In these cases, the probe was focused 
between the specimen and the electron-recording 
screen, so that a bright field shadow image of the 
crystal is formed in the central cone of illumina- 
tion. 
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Fic. 7. ZnO (hexagonal), hh2, dio =1.62+0.01A, 


The field of view was adjusted so that only 
one crystal was present. Through any one point 
of the crystal passes a pencil of rays of very 
small angular aperture (ca. 10~* radian, depend- 
ing on the probe size and the specimen to probe 
separation). If only the rays passing through one 
portion of the crystal were diffracted, there would 
be an uncertainty in the diffraction angle pro- 
portional to the dimensions of the shadow image. 
Actually it is frequently observed in the two- 
dimensional grating patterns that each diffracted 
spot is a dark field shadow image of the diffract- 
ing crystal. Furthermore, the distances between 
corresponding points of a diffracted image and 
the central image correspond precisely to the 
angles of scattering deduced from the inter- 
planar spacings of the crystal; the azimuthal 
directions of scattering are those of the normals 
to the corresponding planes of the crystal. The 
two-dimensional grating patterns may be con- 
sidered as being formed by linear combinations 
(with integer coefficients) of two vectors, Bax: 
and 6by-%-1, of the reciprocal lattice of the crystal. 
The total deviation of any ray that is diffracted 
must be twice the Bragg angle for a particular 
reflection. 

When such two-dimensional grating (cross- 
grating) patterns are obtained, there is an 
apparent breakdown of the Bragg law for the 
diffraction of electrons by the crystal. Although 
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the total deviation of a ray is invariably twice 
the Bragg angle for a particular reflection, the 
beam need not be incident on a crystal at pre- 
cisely the Bragg angle with respect to the re- 
flecting plane. It is sufficient that the ray be 
approximately parallel to the reflecting plane. 

It is evident that in order for a cross-grating 
pattern to be possible, two of the shorter vectors 
of the reciprocal lattice must be approximately 
perpendicular to the electron beam; an equiva- 
lent statement is that an important lattice vector 
must be approximately parallel to the electron 
beam. 

The method of specimen dispersion generally 
used'"—spreading a drop of an amyl acetate 
solution of collodion containing a fraction of a 
milligram of the powdered specimen between 
two microscope slides by sliding them against 
each other—favors orientation of the crystals 
with well developed faces parallel to the collodion 
film. For the crystals used, this condition is 
equivalent to having an important lattice vector 
approximately parallel to the electron beam. 

It should be emphasized that the intensity 
diffracted onto most of the spots observed in the 
cross-grating patterns is quite low—from 0.1 
percent to 3 percent of the incident intensity. 
Sometimes one or two of the reflections are 
very strong containing practically all of the 
scattered intensity. The crystal is presumably 
oriented so that the beam makes an angle close 
to the Bragg angle with the plane responsible 
for these reflections. Figures 5 and 6 provide 
examples of this phenomenon. 

As will be discussed in more detail later on, 
calculations based on the dynamical theory of 
electron diffraction indicate that diffracted in- 
tensities of the order of magnitude of from 0.1 
percent to 3 percent are to be expected for each 
of a set of spots of a two-dimensional grating 
pattern when an electron beam is incident 
approximately parallel to an important lattice 
vector of a perfect crystal. In addition, any 
disorder or warping of the crystal will contribute 
to the development of a two-dimensional grating 


2 This technique was used in all the examples cited in 
this article except for the case of aAlxO3-H2O. An aqueous 
suspension of this material dries down on a collodion film 
with most of the crystals lying flat on the film. 
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pattern, a factor which is probably important 
for most of the patterns exhibited here. 

Qualitatively, it is observed that the relative 
intensities of the various diffraction spots of a 
two-dimensional grating pattern are not the same 
as the x-ray diffraction intensities for the given 
material. In particular, reflections which are very 
weak for the x-ray case may have a much greater 
relative intensity for electron diffraction. Of 
course, these anomalous intensities have been 
observed frequently before in electron diffraction 
patterns. A reflection which would be weak or 
vanishing by the structure-factor criterion, but 
which can be formed by the vector sum of two 
strong reflections will have appreciable intensity 
because it will be formed by the process of 
multiple scattering. This is possible for electron 
diffraction because a beam need not be incident 
at the Bragg angle with a plane to be strongly 
diffracted by that plane. 

The most striking examples of this phenom- 
enon of anomalous intensities among the samples 
studied occurred for BeO. The 404 reflections of 
Figs. 3 and 4 are as strong as the other reflec- 
tions, although the hk4 reflections of BeO are 
extremely weak in the case of x-ray diffraction." 
The 301 and 303 reflections are weak but evident 
(Fig. 3) although the structure factor for these 
reflections is zero. The 301 reflection, for example, 
may be considered as resulting from the combi- 
nation of 200 and 101, or 300 and 001, etc., all 
of which are allowed and strong. However, in 
Fig. 4 there are only hh 2] reflections, and no 
hh 21+-1. In this case there is no way of com- 
bining two allowed reflections to obtain a for- 
bidden reflection. 

It should be recalled at this point that Bethe 
in his fundamental paper on the dynamical 
theory of electron diffraction’ suggested the 
possibility of these ‘“‘anomalous’’ reflections, re- 
sulting from the combination of two strong 
reflections." 

In the diffraction patterns from 
aAl,03:H.O, the symmetry of the orthorhombic 


x-ray 


3 e.9., L. H. Germer, Phys. Rev. 61, 309 (1942). 

4 W. H. Zachariasen, Zeits. f. Physik 40, 637 (1927), 
lists as the relative intensities of the 002 and 004 reflec- 
tions for BeO, 2.5 and 0-0.5. 

‘6H. Bethe, Ann. d. Physik 87, 55 (1928). 

16 Thid., pp. 97-99. 
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crystal prohibits all but the 2h 0 2/ reflections of 
the h0O/ series.’ However, in longer exposures 
images appear in the middle of the rectangles 
that form the orthorhombic pattern. The lattice 
constants of aAl,O3-H.O area =3.72A, b=12.5A, 
c=2.83A. Because of the large value of b, the 
hil reflections occur through nearly the same 
angle, and at approximately the same angle of 
incidence as the hol reflections. Therefore, the 
indices of the ‘‘anomalous’”’ spots are probably 
2h+1, 1, 2/+1. 

The second general class of phenomena ob- 
served are the production of Kikuchi-line pat- 
terns. In general, three requirements must be 
fulfilled in order for the Kikuchi-line phenomena 
to occur: 


(a) The crystal must be sufficiently thick so that the 
primary process of incoherent scattering—to produce a 
cone of rays of much larger angular aperture than the 
incident cone of illumination—take place with sufficient 
probability to produce a general blackening of the entire 
photographic plate. In order that all the background 
scattering should contribute to the patterns, the electron 
probe is focused on the crystal in these experiments. 

(b) The Bragg law must be satisfied for the diffraction 
of electrons by the crystal; that is, electrons incident upon 
a reflecting plane at the Bragg angle must be diffracted 
with much greater intensity than electrons incident at 
another angle. 

(c) Some planes must be oriented approximately parallel 
to the axis of the electron beam. This requirement will 
usually be fulfilled because of the great number of reflecting 
planes in a crystal. 


The distance from a ‘“‘depleted’’ line to the 
corresponding ‘‘enhanced”’ line of a Kikuchi 
pattern corresponds to an angle of scattering of 
twice the Bragg angle. The positions of the 
various lines of a Kikuchi-diffraction pattern 
allow one to infer not only some of the inter- 
planar spacings of the crystal, but also the 
precise orientation of the crystal with respect to 
the electron beam. Where a ‘‘depleted”’ line runs 
through the central area of irradiation, some of 
the rays of the incident radiation must make an 
angle equal to the Bragg angle with this partic- 
ular plane, and the entire cone of illumination 
makes an angle quite close to the Bragg angle. 
There will then be a weak diffracted image of the 
defining aperture containing the strongly dif- 





17S. Goldstaub, Bull. Soc. Franc. Mineral. 59, 348 (1936). 
See also C. Gottfried, Strukturbericht, 4, 129 (1936). 
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Fic. 8. (a) Kikuchi-line pattern from a BeO crystal, 
with diffracted images of aperture of illumination. Note 
that there is a 110-diffracted image and a broad 110- 


fracted beam caused by the rays that were 
incident at the Bragg angle. In the circular 
images, there may be a set of striations perpen- 
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(d) 


Kikuchi line that does not go through this image. (b) 
Deduced orientation of the reciprocal vectors of the BeO 
lattice with respect to axes in space. The electron beam 
is moving in the direction of the +2-axis. (c) Calculated 
Kikuchi-line pattern from the transformation equation. 
All the calculated lines _are_present in the photograph 
except in cases (such as 3 1 1 or 0 4 1) where the + and 
— lines are almost symmetrically disposed with respect to 
the central unscattered beam. (d) An electron micrograph 
ofjthe crystal responsible for this pattern. It seems likely 
that the crystal is oriented with the hexagonal axis perpen- 
dicular to the collodion film; because of proximity to the 
supporting mesh, the film is not perpendicular to the elec- 
tron beam direction. 


dicular to the direction from the central image 
to the diffracted image. (These parallel striations 
will be bent and otherwise perturbed when 
another Kikuchi line runs through them.) This 
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(b) 


Fic. 9. (a) and (b) show the Kikuchi-line pattern from 
an MgO micro-crystal. The deduced orientation is ex- 
hibited in (c) and (d). In this case, the unit vectors 
perpendicular to the faces, rather than the reciprocal 
vectors, are used in the transformation equations. The 
specimen was prepared by dispersing MgO smoke in 


phenomenon was observed for mica by Kossel 
and Modllenstedt, and interpreted by MacGil- 
lavry'® in terms of the dynamical theory of 
electron diffraction. The spacings of these stria- 





18C, H. MacGillavry, Physica 7, 329 (1940). 
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collodion; this procedure promoted the isolation of single 
MgO crystals. 

Note the difference in the (420) striation patterns of 
(a) and (b), two different exposures, indicating a very 
slight change in orientation between exposures. The 
calculated and observed Kikuchi patterns are in agreement 
except that the reflections with odd indices are absent 
except for the case of (5 3 1), where the plane is at its 
Bragg angle with respect to the incident beam. 


tions permit a determination of the thickness of 
the crystal (assuming that it is of uniform thick- 
ness in the direction of the electron beam), and 
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in favorable cases, of the structure factor for the 
particular reflection. 

Figures 8, 9 and 5c are examples of the Kikuchi 
lines produced when the electron probe is concen- 
trated on a micro-crystal. Calculated patterns 
and the inferred orientation of the crystal are 
also illustrated for Figs. 8 and 9. The orientations 
were deduced by requiring that two of the planes 
whose negative Kikuchi lines run through the 
central cone of illumination make an angle equal 
to their respective Bragg angles with the z-axis 
(the direction of the electron beam in Figs. 8b 
and 9c). The predicted patterns are calculated 
by finding by inspection all the linear combina- 
tions (with integer coefficients) of the unit 
reciprocal vectors that have small negative com- 
ponents along the z-axis. In all cases, the meas- 
ured spacings and angles between Kikuchi lines 
were correct within experimental errors (+0.2—5 
percent for the spacings, depending on the sharp- 
ness of line ; +ca. 1°-2° for the angles, depending 
on the length and sharpness of the lines). 

In the early stages of the present work, when 
the electron probe was focused on a specimen or 
a membrane, a contamination “‘spot’’ formed 
rapidly on the irradiated area. Figure 6 (c) is 
a conventional electron micrograph of a pair of 
Al,O;-H.0 crystals that had been exposed to an 
electron beam of rather large area for several 
minutes in the diffraction unit; the contamina- 
tion spot, defining the area of irradiation, ap- 
proximately 2.5 microns in diameter, is clearly 
visible. In practice, the area irradiated with a 
concentrated probe would become completely 
opaque to electrons after 30-60 seconds of bom- 
bardment. It was not possible, therefore, to 
study systematically the production of Kikuchi 
lines by micro-crystals until a method of pre- 
venting contamination was developed. 

_ This difficulty in “probe” analysis was previ- 
ously encountered in the electron microanalyzer.'® 
The carbonaceous nature of the deposit is 
demonstrated in this instrument; the inelastic 
scattering spectra of thin films of silver or silica?® 
showed no peak caused by, the K-level of carbon 
when first exposed to a probe of ca. 2000A 
diameter. In later exposures the carbon K-line 


‘8 J. Hillier and R. F. Baker, J. App. Phys. 15, 671 (1944). 
20 Silica experiment by R. F. Baker. 
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developed, and by the time the specimen was 
opaque only the carbon line was observed. 

It has been discovered that, in the diffraction 
unit or in the electron microanalyzer, contami- 
nation can be avoided if the specimen chamber 
of the unit is exposed to an electron spray 
(60ma at 800v) for a few minutes before concen- 
trating the probe on the specimen. Contamina- 
tion will be entirely negligible for 10 minutes or 
so, and not serious for several days after this 
treatment. It should be mentioned that this 
technique is not effective in preventing contami- 
nation of electron microscope specimens. Here 
contamination appears to be a somewhat more 
complicated phenomenon. 

In early stages of the work, before the method 
of avoiding contamination was discovered, pat- 
terns of a few spots were frequently obtained 
when the probe was focused behind or in front 
of the specimen. If the diffracted spot does not 
form a complete image of the crystal, and if a 
complete two-dimensional array is not obtained, 
it is difficult to interpret the pattern. One reason 
for this is that the particular portion of the 
crystal that is diffracting the beam is unknown, 
so that a precise measurement of the angle of 
scattering is not feasible. It is believed, however, 
that most of such crystals would give useful 
Kikuchi patterns which provide precise informa- 
tion about the inter-planar spacings and angles. 

In concluding the description of the experi- 
mental work, it should be pointed out that the 
work to date has been confined to the study of 
the patterns produced by crystals of known 
structure. It is evident in the Kikuchi patterns 
of Figs. 8 and 9 that even a highly symmetrical 
and simple crystal can give rise to a quite 
complicated: electron diffraction pattern. The 
application of the method to the identification 
of unknown crystals or the determination of 
their lattice constants, would be greatly facili- 
tated by the development of a stage which 
allows precise adjustment of the orientation of 
the crystal with respect to the electron beam 
without moving it from the field of view. 


IV. THEORETICAL INTERPRETATION 
The spacings of the striations in the 420 
reflection of MgO in Fig. 9(a) were compared 
with the predictions based on the dynamical 
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TaBLE I. Observed and calculated displacements between 
maxima of striations of Fig. 9a. 





Maxima involved 1-0 0-1 1-2 2-3 3-4 4-5 
Observed displacement ( 10+) 13.0 13.45 7.3 7.8 69 i 
Calculated displacement 13.1 13.1 7.9 74 7.0 7.0 


theory of electron diffraction. The formula given 
by MacGillavry'® for the angular variation of 
diffracted intensity is: 


1 Dope 4r°A’K? \3 
sin : (1+ ) 
I(A) 2K neiVnne 
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where A=deviation of angle of incidence from 
Bragg angle, in radians 


K=22/\=wave vector of electron beam, 
D= thickness of crystal, 
d)x:=inter-planar spacing, 
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In this last relation, R is the volume of the unit 
cell, and >>; the electronic structure factor for 
the hkl reflection. The total deviation of the 
rays in this treatment is always 263. 

The formula (1) is applicable to a crystal slab 
of uniform thickness in the direction of the 
electron beam. It is assumed that only one 
strong reflection occurs. There are clearly two 
strong reflections in Fig. 9. In Fig. 9(a), the 
striations are sensibly parallel and perpendicular 
to the direction of scattering so that this pattern 
was chosen for analysis. It is interesting that in 
the next exposure (a minute or so later) (Fig. 
9(b)) the 420 pattern was highly perturbed by 
interaction with the 262 reflection; evidently a 
very slight change in the orientation of the crystal 
took place. The next exposure showed parallel 
striations again. If in formula (1) we set A,=d/D, 


; | A’ ume) 
ain*{ e§ —+——_—_- | 
T(A) A? 47n°A?K? 


i 422K? 
It - 


Dnyxr*Vrxi? 











VOLUME 18, JUNE, 1947 


The quantity A; is evidently the limiting value 
of the angular displacement in radians between 
the intensity maxima of the striations as A 
increases. Table I lists the measured values of 
the displacements (in radians) between maxima. 
Taking A,;=6.7+0.5X10-', and Ajp—A;=13.2 
+0.2x10-, one calculates that D=1400+100A 
and v429= 7.742.010" cm~. The value calcu- 
lated from the definition of v;,. is 7.85 X 10% cm. 
Using this value and the value 1400A for the 
thickness of the crystal, the calculated displace- 
ments between maxima of Table | are obtained. 
Evidently the greatest distance is between the 
central maximum and the next one; all the other 
distances are rather close to the limiting value. 
For a strong reflection the method provides a 
reasonably accurate determination of the struc- 
ture factor for a reflection ; for a weak reflection 
the percentage accuracy would be very low. 

In the particular example cited here, no elec- 
tron micrograph was obtained to substantiate 
the assumption of uniform thickness of the 
crystal. The agreement of the results with theory, 
however, supports this assumption. 

We may now consider why in some cases 
cross-grating patterns, in which the Bragg law 
of incidence is apparently not obeyed, were 
observed and in other cases, Kikuchi-line pat- 
terns, which depend on those electrons incident 
at or very close to the Bragg angle being much 
more strongly diffracted than those incident at 
other angles. We first note that there is no 
logical contradiction between the occurrence of a 
strong reflection of rays incident at the Bragg 
angle, and weak reflections of rays incident at 
other angles for a particular crystal. Figure 5 is 
an example of a two-dimensional grating pattern 
and a Kikuchi-line pattern obtained from the 
same béhmite crystal; from the pattern of Fig. 8 
for BeO we see that there is a 110 reflection of 
the entire aperture of illumination and displaced 
from that a broad Kikuchi line because of 
electrons inelastically scattered onto a direction 
making the Bragg angle with the 110 plane. 

It is beyond the scope of the present paper to 
attempt to assess exactly the effect of disorder 
in a crystal on its diffraction pattern. In general 
the effect of such disorder will be to increase the 
intensity of diffraction of rays not incident at 
the Bragg angle. 
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TABLE II. ‘“‘Rocking crystal” curve for MgO, 200 reflection. 


5 1.285 5 
4 0 10) = X10 10°? KO K 10? 
I/Io 1/2 0.485 0.29 0.127 0.085 0.068 


Do (period of 
sin? term, A) 360 354 274 254 181 149 


The dynamic theory of electron diffraction, 
however, does indicate that even for perfect 
crystals the kinds of diffraction pattern reported 
in this paper are to be expected. MacGillavry’s 
formula (1) for the intensity of a reflection gives 
for the half-width of incidence of a reflection 
(neglecting the sin? term) 


Ay ~~ Varidnkt, 
Ay~ (Le acidan’, 


so that wide angle reflections give the best 
Kikuchi lines, whereas the inner reflections give 
broad Kikuchi lines and readily contribute to 
cross-grating patterns. This relation is evident 
in the 110 and 220 Kikuchi lines of BeO in Fig. 8. 

As an example of the intensity of diffraction 
for angles of incidence that are considerably 
different from the Bragg angle, it is of interest to 
apply MacGillavry’s formula to calculate a 
“rocking curve”’ for the 200 reflection of MgO, 
using A=5.4X10-'© cm. In this calculation 
(Table II), the average value of the sin? term is 
taken as }. 

Even at normal incidence (A=6,) the dif- 
fracted intensity is great enough to give notice- 
able reflection. One similarly computes for the in- 
tensities of the 220 and 400 reflections at A = 6229 
and A= 64009 the values 0.011 and 0.0013, respec- 
tively. 

Miss MacGillavry’s formula (1) considers only 
the interaction between one diffracted ray and 
the incident ray, the treatment of the general 
case of the interaction of all the diffracted rays 
being very complicated. It is of interest, however, 
to present calculations for the symmetrical case 
of normal incidence on a cube face of MgO 
considering the interaction of all four of the 200 
reflections and all four of the 220 reflections. 
This is a simple computation and serves to 
indicate if formula (1) is seriously wrong or if it 
gives results of the same order of magnitude as 
those obtained by more elaborate calculations. 
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Only an outline of the methods of calculation 
and the results need be presented.*! 
Bethe’s dynamical equations for this case are: 


(K? — Ro007) Wooo t4v200¥200 + 40220220 =0, 


4r° 
troobout ( K*— hoes! ey +t 20a am 


doo" 
+ (2v200+ 20420) P220=9, 
Vo20Wo00t (2v200+ 27420) W200 


; 4x 
+ («: — Roo0? rad 20+ 000 Wan =U. 
220° 


In these relations, Rooo is the magnitude of the 
wave vector of the undiffracted wave in the 
crystal; its three possible values are determined 
by equating the determinant of the linear equa- 
tions to zero. The y’s are the amplitudes of the 
various waves in the crystal. With numerical 
coefficients inserted, the cubic characteristic 
equation is solved, and the relative amplitudes 
of the internal waves determined for each solu- 
tion. The customary boundary conditions for 
the Laue case lead to the intensity formulae: 
(D in Angstrom units) 


Too = 1 — 0.904 sin*?(r#D/ 146) —0.032 sin*?(#D/125) 
— 0.060 sin?(7D/67.4), 

Too9 = 0.192 sin?(x#D / 146) +0.0088 sin?(#D/125) 
— 0.0083 sin*(#D/67.4), 

T2299 = 0.037 sin?(#D/146) —0.0077 sin?(#D/125) 
+0.0098 sin?(#D/67.4). 


The average value for the intensity of the 
diffracted waves is Io992¥0.10, I2292+0.02. 

The precise numerical values obtained are of 
no great significance because of the very special 
case treated; what is important is that these 
values are of the same order of magnitude as 
those obtained from formula (1) so that this 
formula may be used to indicate reliably the 
extent of a cross-grating pattern for a perfect 
crystal. 

An outstanding deficiency of the dynamic 
theory as presented is that it neglects the effect 
of incoherent scattering on the propagation of a 


21 See references 15 and 18 for a complete discussion. 
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wave system in a crystal. The energy of the 
electron beam in the crystal is pictured as oscil- 
lating back and forth between the incident and 
the diffracted beams without damping. The 
effect of incoherent scattering will be to increase 
the general background, to tend to destroy the 
diffracted rays—particularly for rays not incident 
at the Bragg angle,—and to cause the develop- 
ment of Kikuchi lines. 


V. CONCLUSIONS 


One can conclude then that two-dimensional 
grating patterns are to be expected for crystals 
oriented with an important lattice vector approx- 
imately parallel to the electron beam, provided 


the crystals are not too thick and independent of 
whether or not the crystals are slightly disordered 
or warped. Kikuchi-line patterns are to be ex- 
pected for most other crystals that are not 
excessively thick or disordered; the wide angle 
reflections will give the sharpest and best Kikuchi 
lines. Because of the shape of the “rocking 
crystal’ reflection curves, crystals which give 
rise to two-dimensional grating patterns may 
also give rise to Kikuchi-line patterns. 
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Microsecond Measurement of the Phosphorescence of X-Ray Fluorescent Screens* 


Fitz-HuGH MARSHALL 
Westinghouse Research Laboratory, East Pittsburgh, Pennsylvania 


(Received February 4, 1947) 


Microsecond square-wave x-ray pulses, repeated a thousand times a second, were obtained 
by operating an x-ray tube directly from a radar pulser. The phosphorescent build-up and 
decay of light emitted by fluorescent materials when irradiated by the x-ray pulses was observed 
with a multiplier phototube connected to an oscilloscope. Oscillograms of CaWO, (radio- 
graphic intensifying screen) and BaSO, indicate simple exponential response, with decay to 
1/eth intensity in 6 and 0.8 microseconds, respectively. MgW0O, and Patterson B (fluoroscopic) 
and D (photofluorographic) screens were found to complete their main response to microsecond 
excitation in 10 to 100 microseconds, but their decay is not exponential. 


SUSTAINED interest in the subject of 

phosphorescence during the last decade 
has led to extensive experimental work, though 
rarely with x-ray excitation, and to various 
theoretical proposals to account for the many 
cases in which the observed decay does not 
follow a simple exponential law as would be 
expected for a simple monomolecular process.'~4 
Attempts have been made to interpret these 
more complicated cases as composites of super- 
imposed exponential curves, as hyperbolic decay 
associated with a bimolecular mechanism, or as 
due to emission delays by electron traps. The 
present experiment is not sufficiently compre- 
hensive to test the validity of these theories, for 
it had as its immediate purpose only the deter- 
mination of the speed of response of standard 
types of commercial fluorescent x-ray screens. 
But the experimental method employed has 
proved so simple and effective that others may 
wish to develop it for a more detailed treatment, 
either for the study of x-ray screens as such or 
for a more fundamental investigation of phos- 
phorescence. The microsecond timing technique 
used, instead of the more common millisecond 
timing, is clearly necessary to resolve rapid 
phosphorescence effects, which should not be 
disregarded in a general study. The use of 





* Reported at the Chicago meeting of the Am. Phys. 
“Soc., June 22, 1946 (Phys. Rev. 70, 114 (1946)). 

1 J. T. Randall, M. H. F. Wilkens, and G. F. J. Garlick, 
Proc. Roy. Soc. A184, 347-433 (1945). 

2 J. W. Strange and S. T. Henderson, Proc. Phys. Soc. 
London 58, 369-401 (1946). 

3R. T. Ellickson and W. L. Parker, Phys. Rev. 70, 
290-299 (1946). 

*W. B. Nottingham, R. B. Nelson, and R. P. Johnson, 
J. App. Phys. 10, 335-342 (1939). 
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x-rays, instead of the more usual ultraviolet or 
cathode-ray excitation, not only provides new 
test conditions but also simplifies absorption cor- 
rections when they are necessary, as in hyper- 
bolic decay: in general the x-ray correction is 
less, can be determined accurately and easily, 
and does not involve taking into account a large 
change of intensity within individual crystals. 
The equipment used in this study of the short 
persistence phosphorescence of x-ray fluorescent 
screens is illustrated in the block diagram of 
Fig. 1. An x-ray tube was operated directly from 
a radar pulser. The pulser provided microsecond 
square-wave potential pulses at 10 to 30 kv, 
repeated a thousand times a second. The x-ray 
tube current during the pulse was as high as 
500 ma, but the corresponding average current 
was well within the rating of the air-cooled 
fluoroscopic tube used. With a water- or oil- 
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Fic. 1. Equipment for microsecond measurement of the 
phosphorescence of x-ray fluorescent screen materials. 
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cooled tube, much larger currents could have 
been obtained if needed, though probably with 
some reduction in the filament life because of 
operation at abnormally high temperatures. The 
x-rays irradiated: the fluorescent sample, which 
for simplicity was taped directly against the 
R.C.A. multiplier phototube used to measure 
the light emitted, the phototube then being 
wrapped in black paper to exclude room light. 
The distance from the x-ray target to the fluo- 
rescent sample was about 15 cm. Stray high- 
frequency pickup from the pulser was practically 
eliminated by (1) covering the phototube with a 
grounded hood of aluminum foil to complete the 
system of standard grounded circuit shields, (2) 
introducing a cathode-follower stage mounted in 
pre-amplifier position directly under the photo- 
tube, and (3) grounding various extra points of 
the equipment for optimum results as determined 
by trial and error. The phototube and cathode- 
follower stage were mounted on top of the pulser 
along with the x-ray tube and its filament isolation 
transformer. The multiplier phototube was oper- 





Fic. 2. 
phosphorescent response of barium sulphate. The total 


time of the trace is 5 microseconds. 


Oscillogram of the growth and decay of the 


ated at about 90 v per stage from a d.c. supply, 
which was well filtered and well regulated. 

The cathode-follower stage was connected 
directly to the amplifier input of a DuMont 248 
oscilloscope, no additional amplification being 
required. This oscilloscope, which has a number 
of special auxiliary circuit features, was used as 
a synchroscope to trigger the pulser a micro- 
second after the beginning of the rapid sweep of 
the cathode-ray beam. The sweep time could be 
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varied and accurate time-marker intervals super- 
imposed on the trace. From the phototube 
through the oscilloscope the response of the 
circuit was checked to be nearly flat to about 
5 megacycles per second. The response of the 
phototube itself is good to about 100 megacycles 
per second. 

The phosphorescent response of the fluorescent 
sample during and after the microsecond of 
irradiation by x-rays thus registered on the 





Fic. 3. Response to direct x-rays. 


oscilloscope screen and was photographed. Addi- 
tional exposures on the same film were used to 
introduce reference traces. Figure 2 is an oscillo- 
gram obtained in this way for barium sulphate, 
the fastest fluorescent material studied. The 
time represented by the complete curve is 5 
microseconds. Trace @ is the main response 
curve, trace 6 is a base line recorded with x-rays 
off and displaced downward to prevent confusion 
from overlapping of the small initial irregularities 
caused by stray pick-up, and traces ¢ are similar 
lines above and below upon which microsecond 
time markers have been introduced for accurate 
time measurement. 

As will be shown, the response for barium 
sulphate is exponential. Following through the 
succession of events recorded in Fig. 2, we see 
that during the first microsecond there are no 
x-rays. Then the x-ray pulse starts, and the 
phosphorescent response begins to rise as the 
material is irradiated at constant intensity for a 
microsecond. During excitation active phospho- 
rescent centers are created at a constant rate. 
But as the number of active centers increases, 
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YOTAL TiME TOTAL TIME TOTAL TIME 
5 MICROSECONDS 3) MIGROSECONDS 100 MICROSECONDS 






BARIUM SU PATE 


. CALCIUM TUNGSTATE 


PATTERSON 6 


PATTERSON D 


MAGNESIUM SUL PHATE 


Fic. 44 Comparative oscillograms of the phosphorescent response of various fluorescent materials to a microsecond 
x-ray pulse. 


so does the rate at which they disappear with equilibrium level d at which the active centers 
the emission of light, this rate being always disappear as fast as they are created. But in 
proportional to the number remaining. If the this case the x-rays are turned off after a micro- 
excitation were to continue indefinitely, the second of irradiation, and the light intensity 
response would exponentially approach that decays exponentially toward zero as the active 
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centers disappear. For barium sulphate the 
phosphorescent period (time for decay to 1/eth 
intensity) is only 0.8 microsecond. 

A direct test of the resolution and general 
effectiveness of the equipment is provided by 
Fig. 3. It so happens that the phototube responds 
weakly to direct x-rays with no intervening 
fluorescent screen. The oscillogram of Fig. 3 is 
for direct x-rays of extremely high intensity. 
The curve differs from a square wave because of 
the presence within the phototube of a weakly 
fluorescent material which phosphoresces expo- 
nentially with a period of 1.2 microseconds. 
When this superimposed response, similar to the 
one in Fig. 2, is subtracted, a reasonably square 
wave remains, resolved to about 0.2 microsecond. 

The broadening of the main trace in all oscillo- 
grams recorded has been studied carefully and 
found to be due to noise of the statistical 
fluctuations in the number of x-ray quanta that 
are effective in producing phototube response. 
The noise is exceptionally high in Fig. 3 for the 
following reason: although very few x-ray quanta 
are absorbed at the photosurface, each quantum 
which is absorbed produces a shower of photo- 
electrons. 

In Fig. 4 is presented a complete set of oscillo- 
grams for the various fluorescent materials 
studied. The first, at the upper left, is a repetition 
of the 5-microsecond curve for barium sulphate 
shown in Fig. 2. On the 31-microsecond oscillo- 
grams, 10-microsecond time markers have been 
introduced in addition to the 1-microsecond 
markers. On the 100-microsecond oscillograms 
only 10-microsecond markers are included. This 
set of oscillograms covers the fastest and slowest 
materials studied and the three most commonly 
used types of x-ray screens. The x-ray intensity 
has been adjusted to bring the maximum re- 
sponse to approximately the same value for all 
curves. 

The barium sulphate (which included 5 percent 
of lead sulphate) and the magnesium tungstate 
were selected from an assortment of available 
experimental powders which had been prepared 
for an unrelated purpose in another laboratory; 
the fluorescent intensity of each of these powders 


6 Fitz-Hugh Marshall, J. W. Coltman, and L. P. Hunter, 
Rev. Sci. Inst. tentatively scheduled for publication in 
July, 1947 issue. 
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was very low in comparison with that of com- 
mercial x-ray screens. The calcium tungstate, 
which fluoresces blue, was a standard commercial 
radiographic intensifying screen; the shape of 
the response curve appeared to be the same for 
the several brands, grades, and thicknesses of 
screens tested, though this has not been checked 
by careful measurement. The Patterson B screen 
fluoresces with the familiar yellow-green of 
fluoroscopic screens, of which it is a standard 
example; it may be classed as a zinc-cadmium 
sulphide screen. The Patterson D screen, which 
fluoresces blue, is one which has been developed 
recently for photofluorography; it may be 
broadly classed as a zinc sulphide screen. 

The simplest kind of phosphorescent response 
is an exponential one, the equations for the 
intensity, J, during”excitation or growth and Ja 






INTENSITY ON LOGARITHMIC SCALE 


0 \ 2 3 4 
TIME IN MICROSECONDS 


Fic. 5. Response curves plotted from measurement of 
5-microsecond oscillograms. 
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during decay, being of the form 
[,=A(1-—e-"T) and J,=Be-"’, 


where ¢ is time, T is the phosphorescent period, 
and A and B are constants. While there is no 
particular reason for expecting exponential re- 
sponse, it is a simple matter to check whether a 
response is exponential from either the growth 
or the decay portions of the curve. In _ this 
experiment the decay portion yielded much 
greater accuracy than the steep, poorly resolved 
growth portion. 

In Figs. 5 and 6 the decay curves for the 
various materials studied are plotted on semi- 
logarithmic paper, on which a simple exponential 
decay would appear as a straight line. The data 
are taken from careful measurement of the 
oscillograms. Intensities have been multiplied by 
arbitrary factors to separate the curves conveni- 
ently. An attempt has been made to adjust the 
time scale so that the decay is shown starting at 
zero time. In Fig. 5 portions of the earlier growth 
curves are included as dotted lines for reference; 
the initial drop-off of the curve for direct x-rays 
is also dotted, for it represents the poorly resolved 
response of the circuit to the cut-off of the square 
wave rather than a phosphorescent decay. The 
spread shown for lower points on some of the 
curves indicates error due to uncertainty in 
measurement of small ordinates on the oscillo- 
grams; the deviation shown is about a thirtieth 
of the width of the recorded trace. Somewhat 
larger errors may very well be present from 
non-linearity in registration of the cathode ray 
tube. 

Within the limits of experimental error, the 
response appears to be exponential for calcium 
tungstate, barium sulphate, and the fluorescence 
of the phototube under direct x-rays, the phos- 
-phorescent periods being 6, 0.8, and 1.2 micro- 
seconds, respectively. The several small but 
clearly discernible deviations from linearity in 
the curves for those materials can be explained 
qualitatively in the following way. A slight bend 
for barium sulphate during the first microsecond 
of decay is attributable to a small component of 
direct x-ray action on the phototube. A bend in 
the calcium tungstate curve at 5 microseconds 
must be ascribed to an overswing in the ampli- 
fying circuit, for at the corresponding time a dip 
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through a minimum is clearly observed on the 
31- and 100-microsecond oscillograms for barium 
sulphate and also on the oscillograms in Fig. 7. 
As for the deviation of the points from a straight 
line at the lower end of the calcium tungstate 
curve, other experiments‘ indicate that calcium 
tungstate is so rapid in its decay that it would 
be uncalled for to suspect a progressive upswing 
in the curve on the basis of these few points of 
large experimental error; however, Strange has 
reported two approximately exponential cathodo- 
luminescent periods for calcium tungstate.’ 

A detailed analysis has not been attempted for 
the magnesium tungstate and sulphide screens. 
Their decays are obviously not exponential. The 
comparatively sharp bends near the start of 
these curves are evidently real, and they suggest 
the presence of superimposed decay curves. The 
later portions of these curves fit simple hyperbolic 
decay a little better than exponential decay, but 
the difference is not considered significant in the 
absence of additional supporting evidence and a 
more complete knowledge of the composition of 
the commercial screens involved. No change in 
the shape of the curves was noticed on the 
oscilloscope with variation of x-ray potential or 
intensity, as might be expected for hyperbolic 
decay. If a more detailed experiment is under- 
taken by others, this point should be checked 
more carefully, along with temperature depend- 
ence. Other studies which might be of interest 
include the response as a function of wave-length 
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Fic. 6. Response curves extended through 100-micro- 
seconds of decay. 
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of light emitted by the screen and the effect of 
infra-red stimulation of the screen.’ 

In interpreting the oscillograms and curves in 
this paper, it is important to note that the 
unusual technique of microsecond excitation has 
altered the shape of all curves which are not 
simple exponential responses in such a way as to 
accentuate rapid phosphorescence effects and to 
suppress slow persistence effects. Thus the results 
may not immediately appear to be consistent 
with observations made under more conventional 
conditions. This was forcibly brought out during 
the experiment by an attempt to observe the 
phosphorescence of willemite, a screen with a 
persistence comparable to the flicker time of the 
eye. Although the willemite could be seen to 
glow brightly when irradiated with x-rays, an 
oscillogram obtained according to the normal 
procedure of this experiment (i.e., one made with 
X-ray intensity adjusted to bring the peak re- 
sponse to the usual amplitude) revealed only the 
response of the phototube to direct x-rays and 
was indistinguishable from an oscillogram for 
direct x-rays without the intervening willemite. 
A difference was finally obtained, as illustrated 
in Fig. 7, by greatly increasing the x-ray intensity 
(and the amplification), past the point of over- 
loading the amplifier at the peak. In the second 
part of Fig. 7, the main evidence of the lumi- 
nescence of the willemite is the broadening of the 
curve by the high noise level of the light, which 
is intense but which does not decrease appreci- 
ably during the short interval between pulses. 

Let us assume for simplicity that the reci- 
procity law holds perfectly for each of the 
components involved in the oscillogram at the 
right in Fig. 7, that is, that the total area under 
the curve (extended to infinite time) associated 
with each component is independent of the time 
of x-ray excitation so long as the total quantity 
of x-rays remains unchanged. It is evident that 
if the same quantity of x-rays had been intro- 
duced during a pulse time of 10 microseconds, 
instead of one microsecond, the willemite portion 
of the response would not have been appreciably 
altered, but the direct x-ray portion of the 
response would have been spread over approxi- 
mately ten times as much time and would 
therefore have risen only a tenth as high above 
the willemite response. With a pulse time of a 
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Fic. 7. Comparative 100-microsecond oscillograms for 
direct x-rays and willemite. The amplifier gain and the 
x-ray intensity have been so greatly increased as to cause 
overloading at the peaks. 


thousandth of a second, such as might ordinarily 
be used for studying willemite, it is doubtful that 
any direct x-ray component could be detected, 
even with the x-ray intensity increased further 
to bring the peak response to the usual level. 

This dependence of the shape of the response 
curve upon the length of time of excitation 
might provide the basis of a method for isolating 
component decays in composite curves. A pulser 
with variable pulse time would permit observa- 
tion of the change of shape of the response curve 
with pulse time. A further useful innovation 
might be that of Strange of balancing out the 
main part of the response with counter expo- 
nential condenser discharge curves and meas- 
uring only the difference.? It should be noted 
that hyperbolic decay curves would behave 
differently from composite exponential curves, 
there being an actual energy shift from long- 
persistence effects to short-persistence effects as 
the pulse length is shortened. 

It should be mentioned that in the present 
experiment all persistence after a thousandth of 
a second is rendered completely indistinguishable 
by the beginning of the next sweep cycle and the 
advent of a new x-ray pulse; the zero ordinate 
point seen at the beginning of each oscillogram 
trace is merely the response level at the end of 
the previous millisecond cycle. 

Let us finally consider the practical importance 
of the results of this experiment with respect to 
the speed of résponse of commercial x-ray screens. 
It is necessary to introduce supplementary infor- 
mation from other sources, inasmuch as long- 
persistence effects of relatively high energy might 
be present without detection by the present 
experiment. 
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Under ordinary conditions calcium tungstate 
does not show any detectible persistence. When 
a calcium tungstate screen is whirled as a disk 
before a narrow beam of x-rays, there is no 
noticeable blurring of the edge of the primary 
spot on the screen covered by the x-ray beam. 
We may therefore conclude that the curves for 
calcium tungstate recorded in this experiment 
include all of the response which is of practical 
importance and that the response is almost 
complete in 10~° second. 

In contrast, the sulphide screens, especially 
Patterson D, are commonly known to have long- 
persistence effects which’ are visible for some 
seconds or minutes after exposure to intense 
x-rays. Since these are not detected in the 
present experiment, it is necessary to consider 
whether a large fraction of the total light emitted 
by the screen is associated with the long- 
persistence effects. The case for Patterson D is 
of greatest practical importance, for it is usually 
found to be the brightest screen as measured 
photographically or photoelectrically. A rotating- 
disk test shows the primary spot on the screen 
to be bright and sharply delineated from the 
phosphorescent trail. On the basis of this and 
other tests, it appears safe to say that the main 
useful part of the response (perhaps three-fourths 
for Patterson D, more for Patterson B) is 
included in the oscillograms of this experiment 
and is complete in 20 microseconds or so. Thus 
the amount of long-persistence light is too small 
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. obtained by 


to be especially important from the standpoint 
of efficiency; on the other hand, it is large 
enough to be a disturbing factor in fast response 
applications for which it is essential that there 
be no detectible residual from earlier x-ray events. 

As a practical case in point, we may consider 
X-ray motion pictures. These are now generally 
photographing the fluorescent 
screen. When Patterson D is used with conven- 
tional x-ray and camera equipment, noticeable 
blurring or multiple-image effects occur, especi- 
ally at the leading edge of a dark area moving 
into a bright area. This blurring does not occur 
when calcium tungstate screen is used, but the 
effective intensity is less. 

These conclusions also apply to the photo- 
multiplier x-ray detector. This detector consists 
of a fluorescent screen wrapped around a multi- 
plier phototube and is therefore equivalent to 
the x-ray detector used in this experiment. Speed 
of response is important in the analysis of fast 
x-ray phenomena. Here again the rapid response 
and high intensity of Patterson D make it very 
suitable for many applications. - However, the 
long-persistence effects of this screen cannot 
always be disregarded. For example, when the 
detector is used as an integrating exposure meter 
to terminate a radiographic exposure after the 
incidence of a given quantity of x-rays, the 
long-persistence effects result in an objectionable 
dependency upon exposure time and upon recent 
exposure history. 
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X-Ray Diffraction Studies of Chrome-Steel Slags' 


G. P. CHATTERJEE AND S. S. SipHu 
University of Pittsburgh, Pittsburgh, Pennsylvania 
(Received March 21, 1947) 


X-ray diffraction studies of the crystalline constituents in the acid open hearth chrome- 
steel samples showed that chromium exists in these samples as a chrome-iron spinel of the 
form, FeO-Cr2O;. Its density was 5.109 g/cc and its melting 2160°C. It crystallized in the 
face-centered cubic structure with ap>=8-348A and Z=8. 


I. INTRODUCTION 


HE essential major elements in most of the 

common alloy steels are chromium (Cr), 
nickel (Ni), and manganese (Mn). Nickel is not 
oxidized during steel making. Of the other two 
elements, Cr is more costly, and its loss in slag 
through oxidation or any other combination 
increases the cost of steel to a greater extent than 
an equivalent loss of Mn. The nature of control 
methods, which may lead to decreased losses of 
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Fic. 1. The schemetic diagram of the small indirect arc 
furnace for preparing synthetic spinel. T=top electrode; 
C=rotating crucible; F =fused charge; B = unfused charge; 
R=variable resistance; and A =ammeter. 


' The paper was presented, in part, at the joint meeting 
of the American Society for X-Ray and Electron Diffraction 
and the Electron Microscope Society of America, Pitts- 
burgh, Pennsylvania (Dec. 1946). 
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Cr, depends on the knowledge of the different 
phases in which Cr is present in the slag. The 
object of this investigation was to study by 
x-ray diffraction methods if Cr exists in the acid 
open hearth chrome-steel slag as (a) free Cr.O; 
or CrO, (b) solute in FeO solvent, (c) silicates 
of chromium, and/or (d) spinels of the type 
FeO -Cr,O; or CrO-FesQ3. 


II. EXPERIMENTAL PROCEDURE 


Acid open hearth chrome-steel slag samples 
were obtained from regular heats of steel and 
poured into Fluidimeter.? The slag solidified and 
cooled rapidly inside the mould to preserve 
nearly the original conditions in which the 
different constituents existed in the molten slag. 
From the chemical analyses of these samples, 
samples of synthetic slag were prepared. A 
mixture of pure FeO and Cr.QO; (mole proportion 
by weight) was put in graphite crucible C of a 
small 250-volt d.c. arc furnace as shown in Fig. 1. 
A small amount of borax was added to facilitate 
rapid melting. The top electrode T was placed 
slightly eccentric with respect to the bore of the 
crucible. An indirect arc was struck on the side 
of the crucible while the crucible was rotated. 
In one complete rotation of the crucible, the 
entire mixture was melted beneath the arc. The 
time of rotation was 20 seconds and the current 
through the electrodes was 10 amperes. During 
the arcing the electrode T was kept negative and 
the crucible positive. It was found by optical 
pyrometer that the positive electrode was always 
at higher temperature than the negative elec- 
trode. 


* Fluidimeter designed by the Acid Open Hearth Re- 
search Association under the direction of Dr. G. R. 
Fitterer, University of Pittsburgh, Pittsburgh, Pennsyl- 
vania. 
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To remove borates and other contaminants 


formed during melting, the synthetic samples 
were powdered and boiled in 1:2 HCI solution 
for a few minutes in a test tube. The solution 
with the residue was filtered and washed several 
times with hot water until the filtrate was free 
from iron salts and chlorides. The residue was 
dried in a muffle furnace at 120°C for one hour. 

X-ray diffraction patterns of both the natural 
and synthetic slag samples were made by the 


TABLE I. Powder diffraction data. 


d =interplanar spacing in A units, //Jo =estimated relative intensity. 


= Very; S =strong; M = medium; W =weak. 
Chrome steel slag Synthetic slag Faylite 
0 ‘ I/Io d 1/Io d 
M 4.82 M 4.83 
M 3.70 - - S 3.71 
. M 3.53 
S(broad) 2.93 M 2.95 S 2.85 
V.S. 2.54 V.S. 2.51 M 2.55 
M M 2.29 
M (broad) 2.08 M 2.08 W 2.06 
S 1.91 S 1.91 V.W. 1.91 
M 1.72 W 1.71 S 1.76 
S 1.61 S 1.61 M 1.61 
W 1.54 M 1.51 
S 1.492 S 1.488 
W 1.442 M 1.42 
W 1.352 V.W. 1.328 W 1.33 
M 1.281 M 1.279 V.W. 1.277 
W 1.216 W 1.210 
V.W. 1.197 
W 1.176 V.W. 1.170 V.W. 1.176 
M 1.111 V.W. 1.115 
W 1.079 M 1.083 
~ W 1.045 
M 0.968 M 0.965 
— W 0.933 
W 0.878 W 0.878 
M 0.858 M 0.854 


usual powder method with filtered MoKa radia- 
tion and a cylindrical camera of 171.9 mm 
effective diameter. 


III. DISCUSSION OF RESULTS 


The x-ray diffraction data obtained from the 
powder patterns are given in Table I. It is 
apparent from the data that the chrome-steel 
slag contains the same compound or compounds 
that are present in the synthetic slag and also 


_ §.109 g/cc X (8.348)* x 10~*4 cc X 6.023 X 10%/mole 


small amounts of faylite, (FeO) SiO». The analy- 
sis of the diffraction data of the synthetic slag 
sample shows that it contains no free FeO, 
CroO3, nor a solid solution of FeO and CrsQs, 
but a compound formed from the two oxides. 
That this compound is a spinel similar to 
chromite (FeOQ-Cr.Q3), one of the natural spinels 
reported by Dana* may be shown from its 
crystal structure as determined below: 
For cubic crystal system 


dpi =do .(h?+k?4+/7*)}, 
or 
1 
(h?+k?+/*?) =ao?X ue (1) 


Ann” 


where d,,:= interplanar spacing corresponding to 
plane hkl and ay=side of the unit cubic cell. 
For face-centered cubic crystal, the Miller indices 
of the planes that given reflections and the sum 
of the squares of the Miller indices are given in 
Table Il. By plotting (4?+?+/) as function of 
1 /d*, one gets a straight line as shown in Fig. 2. 
The slope of this line as given in Eq. (1) is ay’. 
An accurate value of da») thus determined is 
8.348A. 

The density of the synthetic slag sample 
determined by the usual picnometer method was 
5.109 g/cc. The molecular weight of the com- 


TABLE II. Interplanar spacings and the corresponding 
Miller indices of the first eight lines of the synthetic 
spinel. 


Miller indices 
hkl 


I/Io dinA 1/d? in A~ (h? +k? +12) 
M 4.83 .043 111 3 
M 2.95 530 110(2) 8 
V.S. 2.51 .158 113 11 
M 2.08 231 100(4) 16 
S * 1.91 .274 133 19 
W 1.71 .342 112(2) 24 
S 1.61 .386 111(3), 115 27 
S 1.48 457 110(4) 32 


pound FeOQ-Cr.O 3 is 223.87 g/mole. Therefore 
the number of molecules per unit cell of the 
compound is 





o) 


= 8.00. 


223.87 g/mole 


3E. S. Dana, A Text-Book of Mineralogy (John Wiley and Sons, New York), p. 418. 
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If, on the other hand, it is assumed that the 


Ye in A? 


weight is 227.71 g/mole the number of mole- 


compound formed is CrO-Fe2O; whose molecular cules per unit cell is 


4 


_ 5.109 g/cc X (8.348)* K 10-*4 cc X 6.023 X 108 /mole 





— = 7.86. 


227.71 g/mole 


Since for a compound the number of molecules 
per unit cell must necessarily be whole number, 
it is evident that the chemical formula of the 
compound formed is FeO-Cr.O; with Fet+ as 
the divalent element and Cr*** as the trivalent 
element. 


IV. CONCLUSIONS 


1. The crystalline constituents in the acid open 
hearth chrome-steel slag samples are (a) chrome- 
iron spinel of the type FeO-Cr.O;' and (b) small 
amounts of (FeQ)sSiQs». 


* More exactly (FeO, MnO)-Cr2QO;, since the acid open 
hearth slag sample contains some amount of MnO which 
is usually present as (MnQO).SiOz dissolved in (FeO) 2SiOxe. 
It is not unlikely, however, that some MnO may replace 
FeO in the spinel FeO-Cr.0;. The atomic weight of Mn 
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2. The density of the spinel as determined by 
the picnometer method is 5.109 g/cc and its 
melting point about 2160°C. 

3. The spinel crystallizes in the face-centered 
cubic system with lattice constant ad» as 8.348A 
and the number of molecules per unit cell is 8. 
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differs only by 0.92 from that of Fe and the atomic scat- 
tering factor of Mn does not differ much from that of Fe, 
the effects of replacing FeO to a small extent by MnO is 
therefore hardly detectable. 
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Spherical Aberration of Compound Magnetic Lenses 


L. Marton* anp K. Bor 
Stanford University, Stanford University, California 
(Received January 15, 1947) 


A reduction of the spherical aberration of strong electron lenses can be achieved by a strong 


lens as a virtual image former and by a transformation of the image in a real one by means 


of one or more weak lenses. Calculations are carried out for bell-shaped magnetic fields of the 


axial field distribution 


H(s)= 


Ho 


1+(s/a)’ 


and numerical values of the achieved reduction of the aberration are given for different lens 


strength, magnification, and image distance of the compound system. 


ARLY in the development of electron 
optics, it was recognized that spherical 
aberration is one of the most serious limiting 
factors in the performance of electron lenses. 
From the beginning, efforts have been directed 
toward a reduction of this defect. One of the 
earliest attempts was that of Scherzer' who 
demonstrated that for a thin and weak electro- 
static lens an axial field distribution of the form 
(z) =goe8*” gives the minimum spherical aber- 
ration. The first similar calculation for a mag- 
netic lens was carried out by Glaser.? Later 
Recknagel* and Plass‘* contributed to the same 
subject. All these calculations deal with very 
weak lenses; the only practical lenses with 
somewhat reduced aberration were designed by 
Gray® and Law.® 
Scherzer demonstrated in another paper’ that 
the spherical aberration of electron lenses can 
never be reduced to zero. Rebsch and Schneider® 
showed furthermore that no combination of 
electric and/or magnetic lenses can have less 
aberration than any of the components alone 
would have. Already at that time, however, 
_Rebsch® demonstrated that, theoretically, the 





- aw; at the National Bureau of Standards, Washington 
| Ee ce 

1Q. Scherzer, Zeits. f. Physik 101, 23 (1936). 

2W. Glaser, Zeits. f. Physik 109, 700 (1938); 116, 19 
(1940). 

3A. Recknazel, Zeits. f. Physik 117, 67 (1940). 

*G.N. Plass, J. App. Phys. 13, 49 (1942); 13, 542 (1942). 

5 F. Gray, Bell Sys. Tech. J. 18, 1 (1939). 

®R. R. Law, Proc. I.R.E. 25, 954 (1937). 

70. Scherzer, Zeits. f. Physik 101, 593 (1936). 

8 R. Rebsch, and W. Schneider, Zeits. f. Physik 107, 
138 (1937). 

® R. Rebsch, Ann. d. Physik 31, 551 (1938). 


522 


aperture defect may be reduced below any pre- 
scribed limit if in a compound lens the physical 
size of the component lenses is reduced propor- 
tionately. For obvious reasons this method is not 
very practical. 

A somewhat different approach to the problem 
of reduction of the spherical aberration of short 
focus, thick electron lenses was briefly reported 
eight years ago by one of us.'® It consists of com- 
bining a strong lens, forming a virtual image, 
with one or more weak lenses which convert the 
virtual image into a real image. This procedure 
is universally adopted for the construction of 
light microscope objectives and other light 
optical systems, and a more thorough inves- 
tigation of its application to electron optical 
systems seems desirable. In what follows, a 
quantitative evaluation of the reduction of the 
spherical aberration of strong magnetic lenses is 
given. For this purpose the theory of strong 
magnetic lenses, developed by Glaser,'' must be 
somewhat extended. Since this discussion is 
limited to the reduction of spherical aberration, 
the reduction of other aberrations by the use of 
compound lenses, which was recently reported,” 
will be completely omitted. 

Glaser has shown" that for practical fields both 
the spherical and chromatic aberrations have 
minima when the object is near the center of the 
lens, for which position the image is virtual for 
practically all lenses. 


TL. Marton, Phys. Rev. 55, 672 (1939). (Abstract: 
more detail is given in U. S. Patent No. 2,233,264.) 

''W. Glaser, Zeits. f. Physik 117, 285 (1941). 

2 J. Hillier, J. App. Phys. 17, 411 (1946). 
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Fic. 1. Representation of 
symbols and coordinates used 
in Case I, a single objective 
Z4 lens. 











Our problem, then, is to determine how small 
is the additional aberration introduced by the 
weak lenses. If it is sufficiently less than the 
decrease achieved by the formation of a virtual 
image by the first lens, our purpose is accom- 
plished. 


BASIC THEORY FOR FIELDS OF FORM 
Ho/(1+ (z/a)?] 
It was first shown by Glaser' that the spherical 
aberration produced by a magnetic lens with 
field distribution H along the z-axis is given by: 


Ar e 7i 7 2e 
—=a' f ( H'+SH!— HH" )yds. 
M 96m V 42, \mV 

(1) 


Here y; represents the path equation of an elec- 
tron satisfying the initial condition: 








yil(zo) =0, yi’ (Zo) = 1. (2) 


As for the other symbols: Ar is the divergence 
in the image plane, z;, of a beam of aperture a 
at zo, M is the magnification, V the electron 
energy in terms of the accelerating voltage, and 
e/m the ratio of charge to mass of the electron. 
If H is given by 


H = Hy/(1+(z/a)*], (3) 
then the path equation is 


y=a(1+(z/a)?)'{D, sinf(1+?)! arccots/a | 
+ Dz cos[(1+?)' arecotz/a]}, (4) 


D, and Dz being arbitrary, while k? = eH;?a?/8m V. 
a is the so-called half-width of the magnetic field. 
k® is a measure of the refracting power. The 
equation is put in a more palatable form if one 
sets z/a=cot¢ (see Fig. 1) and (1+)!=w. One 
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then obtains 


sinw(y~— ¢o) 
y =aC.————_ 


: (5) 


sing 
one constant having been chosen so _ that 
y(zo) =0. For convenience we quote a few more 
formulas from Glaser: 
Zr,=acot(mm/w), 2r;= —acot(nr/w), (6a) 


a ==@ 
f=———, f= (6b) 


sin(n7 /w) sin(n2/w) 





N<w<N-+1, (6c) 


N being the number of focal positions. Here f is 
the focal length and zy is the location of the 
focal point. 

The aberration for a single lens of this type 
has been thoroughly studied, whence we give 
just a resumé of some of Glaser’s work. 

Our problem, to be sure, is to study the aber- 
ration produced by a ¢ompound objective. This, 
however, we can immediately do by extending 
Eq. (1) over the whole system, as long as we 
obey the restriction that the radial distance of a 
paraxial ray be zero at the endpoints of the path 
of integration. 

Case I. Single lens. The equation for the 
spherical aberration is obtained by integrating 
Eq. (1). From it follows that: 





—=a' —-K,, (7) 


M 6w1* sintyo 
where K;, is given by 
¢i 
Ki={ [2(4k,?—3) sin?g+7 | sintwi(y— go)de. 
7 (8) 
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- Spherical Aberration for High Magnification 
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Fic. 2. Spherical aberration vs. lens strength for high magnification and optimum object position. 
(From Glaser.) 


Since we need other upper limits later on, we leave it variable and obtain for the complete 


expression : 





3(4k,?—3) 





Ki(¢o, ¢) = $w17(¢— go) —w1 sin2wi(¢y— go) +41 sin4dwi(e— ¢o) = sin2d¢ 
3 w,! 4k,?*-—3 F 4k,*-3 4k,*-—3 
-— — sin2 ¢9>+————# sin2w;(y¢— ¢o) cos2¢— — cos2w;(y— ¢o) sin2d¢ 
4 k,? 4k,?+3 4k,’ 1? 
14k’?-3 1 4k,?-3 
—— ——-w sin2wi(¢— go) cos2¢-+-— — cos4wi(yg—¢o) sin2dy. (9) 
8 4k,?+3 16 4k,?+3 


For high magnification, however, go=¢r while 
¢i~0, and from the lens equation (6) wer =z, 
so that the expression is much simplified. The 
aberration for go= gr as a function of k? is given 
by Fig. 2. 

The aberration is, to be sure, not a minimum 
for go=¢gr. The presence of the factor sin~*go 
suggests immediately that the object should be 
close to the center, or go>~ 2/2, as was mentioned 
in the introduction. Indeed, Glaser also plots the 
aberration for the optimum object position, 
which we give in Fig. 2. 

For high k*® the aberration is not appreciably 
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altered by picking the optimum position; for 
the lower k® that are actually used in practice 
(1<k?<2), however, improvement is possible. 

There are two objections to having the object 
right at the optimum position : the magnification 
is almost unity for practical values of k, and the 
image is virtual. 

In order to accomplish anything we must 
therefore sacrifice some improvement in_ the 
aberration to gain magnification, and we must 
put in at least a second lens to attain a real 
image. We must, therefore, invesigate next the 
aberration introduced by the second lens. 


JOURNAL OF APPLIED PHYSICS 











St rong Lens 
Hy 











Assumed Field-Free 
Plane in which Paths Mato hed 


Electron Path 





Weak Lens ee 
H 

—— 

he Hoo 




















Fic 3. Representation of symbols and coordinates used in Case II, for a doubly compound objective lens. 


Case II. Two lenses. We first define our axial 
field H as follows: 


(i= Ho [1+(z a)? | for — x <2<383, 


H= (10 
) 


H.=Hy2/\1+[(s—D) /a2 ]?} for 23<2< @. 


D is the distance between lens centers. 
Quite arbitrarily we choose 23 so that 


H,(23)=H2(z3). 


We also assume that the two fields do not interact 
to an appreciable extent. This is done so that 
we can use the path equations of the previous 
section for each field, the constants of the second 
being adjusted te match the first at 23. Only in 
the absence of such interaction are the calcula- 
tions strictly valid." ” 
Equation (1) then becomes 


Ar e 23 24 

—=a' \f f(A) yitde+ f(tts)y.ds 

M 96m VL¥,, 23 J 
(11) 





The integrals can be transformed into that of 
Eq. (8) by adopting a new coordinate system for 
the second lens. We take 2’=z—D and 2'/a2 
=coty’ (Fig. 3). Applying furthermpre the 
initial condition (2) to Eq. (5), we obtain for y; 
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and ye 


ad;  sinw:(¢go— ¢) ) 


’ 





_ 
W1 SiN go sing 

and -(12) 
sinwe(y’ — ¢;’) 
Y2=2C2 ’ 


sing’ 








where ¢;’ is chosen in accordance with Newton’s 
image equations, which hold for this type lens, 
and where C2 is to satisfy the matching con- 
dition. 

Slightly modifying Glaser’s results we then 
obtain 











Ar ayk;? 1 doko” 
=a ne" cK, | 
M 6w14 sintgo 
=a*C,p. (13) 
For K, and Kz we write (see Eq. (8)) 
$3 
Ki=[ f(ki, 9, gode, 
0 (14) 


¢3’ 
K2= -{ f (Re, 9’, oi dg’. 


As we have already seen, (Eq. (9)), the exact 
expressions for K, and Ke are exceedingly cum- 
bersome ; but fortunately they may be approxi- 
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mated—for the values we use—to within a few percent by the following equations: 


3w'4k?—-—3 








K, = $w1?(¢3 — go) + - sin 2 ¢0; 
4k, 4k?+3 
(15) 
4k.?—3 4k.?—3 
Ke= —| 3w2?(¢3’ — ¢:’)+ e we SiN2we( 3’ — ¢;’) cos2¢3' — oe COs2we(¢3’ — ¢;’) sin2¢;’ 


Now we but need an evaluation of C2 to have an analytic expression for the aberration constant, Cyp. 
To do this it is convenient to write y2 as 


- , 
SsINnwe(¢g — ¢o ) 


¥2™a302—————_, (16) 


go being the position of the virtual image formed by the first lens as seen by the second, and again 
being determined by the Newtonian image relations. 

Although we must match both y; and ye and their slopes, it is obvious, go’ having been determined, 
that if C2 satisfies either condition it will automatically satisfy the other. We thus pick as most 
convenient the condition that the path distances (see Eq. (12)) be equal. We have then: 





‘ , ; e 
SINWe(¢3 — ¢o ) a, Sinw3(¢o— ¢3) 
aols Se a (17) 
sin ¢3' W1 SING) SIN g3 
Fortunately, the arguments are such that we can approximate as follows: 
wo(ys3’ — go’) ay T—w3i(go— ¢3) 
a2.0. ——___—__ = —___- -___—_—__- (18) 
; 
T— £3 @) SIND 3 
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— ad WwW 
a= — — — Virtual Image 2 
Pa ‘ 
— — - — Infinite Inmace YO at 5 
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-_ 
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OBJECT POSITION in Terms of §-= Easzi.0,,=fr-\fo 


Fic. 4. Spherical aberration vs. objective position in first lens, for a second lens of varying strength, 
but fixed half-width. An infinite image is formed when the first lens forms a virtual image on the 


focal position of the second lens. Above this boundary, the system forms a real final image, below, 
a virtual one. 
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Fic. 5. Spherical aberration vs. the final image distance for a typical set of conditions. 


Now setting go= ¢r —6, and applying the fact that wer =m (Eq. (6)), gr being the focal position, 


we obtain: 


ay, 1 1+5/¢3 T— 03 


C,.=— 











. (19) 
a2singog we ¢3'—¢0 
Thus we have obtained for the aberration constant of a two-component compound lens: 
Cw 1 ki? ai\*_ ke? ™—~3 \4 
NG XE kd SO bo EN al ek (20) 
a, 6sin‘*¢golw;' a2 we! go — $3 


DISCUSSION 


We have already seen that gp should be close 
to 2. The question remains how the second 
term varies with go. As the beam leaving the 
first lens is converged less, or go—7/2, one can 
expect it to increase rapidly, since the aberration 
increases with the cube of the aperture. The 
behavior with go for different ke is shown in 
Fig. 4. 

For the same reason the two lenses should be 
close together. It is unfortunate that we have not 
succeeded so far in treating the case of interacting 
fields, so that we could find the exact behavior 
when D becomes small. It is certain, however, 
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that practical compound lenses will need to be 
close enough together to make field interaction 
inevitable. 

One reason for the need is the final image 
distance involved. Figure 5 gives an idea of the 
increase in image distance for only a 20 percent 
reduction in the aberration. Decreasing D is the 
only way of decreasing the aberration and 
keeping a reasonable image distance. 

Another way of decreasing C,, is by increasing 
the ratio of the half-width of the second lens to 
that of the first. (See Fig. 6.) But here again we 
increase the final image distance, as is also 
indicated on the graph, in a manner approxi- 
mately proportional to a2/a,. 
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Fic. 6. Spherical aberration vs. object position, for a second lens of varying half-width. The 
total magnification along the curves is infinite, but to gain an idea as to how the image distance 
q’ varied, we took the lens strength indicated by the aberration minimum and calculated q 


for a magnification of 150. 


Thirdly, the contribution of the second lens is 
almost proportional to k2?(k.<1), shown in 
Fig. 4. There again, however, another considera- 
tion enters. The focal length of the second lens 
is inversely proportional to k,”. The greater this 
focal length, the higher must be the magnification 
by the first lens in order to produce a virtual 
image outside the second focal point and thus 
secure a real final image. 

But for high magnification by the first lens one 
achieves almost no reduction in the aberration. 
The result is that there is an optimum value for 
ke, as is shown in Fig. 7 (which is merely Fig. 4 
plotted against total magnification). 

This optimum value, we note, is not neces- 
sarily the best from the point of view of image 
distance. Referring once more to Fig. 5, it is seen 
that especially for high magnification the aber- 
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ration minimum is very flat. One can thus gain a 
considerably shortened image distance, without 
appreciably increasing the aberration, by taking 
a value of k2? somewhat less than the optimum 
value. 

We saw from Fig. 4 that some decrease in the 
aberration is possible using only two lenses, but 
that greater with the 
formation of a virtual image by a very weak 
second lens. The question then arises whether 
one cannot use three lenses to good advantage. 


even decreases come 


Our next problem is to decide therefore 
whether the contribution of the third lens can be 
made sufficiently small to make its addition 
worth while. 

Case III. Three lenses. Again proceeding as in 
the previous section we easily obtain for the 
aberration: 


ky? a, . 6 ‘ k.” © — £3 4 
“nal (CY (2) 8) 
w 4 a2 3 we'X\ ¢; F 


d2\* Ks ks? sintwo(¢s’ — go’) sintgs” 
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Fic. 7. Aberration vs. total magnification. This is a plot of graph II vs. magnification instead of 
object position. 


The symbolism is exactly analogous to that for 
Case II except that the subscript 4 refers to the 
field-free plane between the second and _ third 
lenses. All double primes refer to the coordinate 
system of the third lens. 

Equation (21) leads us to the following quali- 
tative conclusions: 

(1) By making a2z/a3 small enough (say }) we 


can obtain Csp, lenses real image - Csp, lenses virtual image* 


In other words, the last term in Eq. (21) must 
be <1. 

(2) The image distance, g, depends almost 
entirely on the spacing between the last two 
lenses. If we wish to satisfy the condition of no 
interaction between the lenses, g will be greater 
than for the two lens combination. 

Both aberration and image distance, for the 
combination of two or three lenses, show a con- 
siderable improvement only for conditions ex- 
ceeding the limits imposed by our method of 
derivation. It is hoped that a more rigorous 
derivation and/or experiments will enable us in 
the near future to substantiate and to extend 
the results of this preliminary calculation. 
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CONCLUSION 


By combining two lenses one can approach 
the lower limit to the aberration set by the first 
lens, used as a virtual image former. For such 
lens separation for which the interaction of ‘the 
lens fields is negligible, the image distance is con- 
siderably greater than that of the first lens used 
as a real image former. 

How closely one can approach this lower limit 
cannot be answered without a study of the aber- 
ration for very close, interacting fields. 

Certainly it seems desirable to have the lenses 
close together. 

In addition it has been found that there is an 
optimum value for the strength of the second 
lens in a two-lens system, and that the half- 
width of the second field should be several times 
greater than that of the first. 

Three lenses may be used to improve the 
aberration, but only at a cost of a great image 


distance. 
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Physical Properties of Calcium Tungstate X-Ray Screens 


J. W. Cottman, E. G. EBBIGHAUSEN,* AND W. ALTAR 
Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 


(Received January 6, 1947) 


The absolute energy of the light from a fluorescent calcium tungstate x-ray screen, its 
spectral distribution, and the efficiency of conversion from x-rays to light have been determined. 
The absolute energy and its spectral distribution were obtained spectrographically by com- 
parison with a standard tungsten filament lamp. Equations are given for determining the 
distribution of light in the diffusing screen material, and calculation of the light losses are 
made for varying screen thickness and x-ray absorption. The data for energy distribution of 
the ingoing x-rays were taken from measurements obtained in this laboratory by a new method 
to be published elsewhere. Allowance for scattering and secondary x-radiation was made and 
a conversion efficiency of 5.0 percent was obtained for this calcium tungstate phesphor. 


INTRODUCTION 


HE primary purpose of this research was to 


obtain experimentally the absolute effi- 
ciency of the transformation of x-ray energy into 
light energy by calcium tungstate. Of the various 
fluorescent materials which might have been 
examined, calcium tungstate appeared to be the 
one whose investigation would be most fruitful, 
first, because of its wide utilization in radiog- 
raphy, and second, because the fluorescence of 
this phosphor is not due to the presence of a 
foreign material serving as an activator. Thus 
calcium tungstate phosphors produced by dif- 
ferent manufacturers should the same 
spectral distribution, and probably will not vary 
greatly in their efficiency. As a basis for theo- 
retical studies of fluorescence, knowledge of the 
behavior of calcium tungstate is probably of 
more value than that for the sulphide phosphors 
which are very sensitive to minute quantities of 
impurities. 


have 


The problem of determining the conversion 
efficiency divides itself into three main parts. 
First, the absolute energy and spectral distribu- 
tion of the light emitted from the screen must 
be determined. Second, the light losses in the 
screen must be determined in order that the 
total light generated may be calculated. Third, 
-the energy of x-rays absorbed in the screen must 
be measured. The resultant efficiency will be a 
characteristic of the calcium tungstate itself and 
of the wave-length of the exciting x-rays. With 
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this information and a knowledge of the func- 
tional dependence of the light losses on the 
structure of the screen as presented herein, the 
behavior of any calcium tungstate screen may 
be predicted. The method of analysis developed 
to obtain the light losses are applicable to the 
transmission of light in diffuse media in general 
and, in particular, to phosphors other than cal- 
cium tungstate. 


Ia. LIGHT ENERGY AND ITS SPECTRAL 
DISTRIBUTION 


The spectrographic part of the problem con- 
sisted of the determination of the spectral energy 
distribution of the CaWO, screen when excited 
by 80 kv x-rays, filtered by one half millimeter 
of Cu and one millimeter of Al, and with a tube 
current of 20 milliamperes. The procedure in- 
volved a photographic comparison between the 
known spectral energy distribution of a tem- 
perature-calibrated tungsten ribbon filament 
lamp and that of the CaWO, screen. 


The Experimental Apparatus 


The standard source was a lamp with a tung- 
sten ribbon filament and quartz window loaned 
by the National Bureau of Standards. The ribbon 
filament is about 1} inches long and } inch wide 
and is supported in the bulb by heavy nickel 
leads. 

The temperature calibration of the lamp was 
made by Mr. M.S. Van Dusen, Acting Chief of 
the Pyrometry Section, who gave the value of 
the brightness temperature outside of the quartz 
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window as 2614°K at a current of 39.20 amperes 
and for a point about midway along the filament. 
The position is marked by a small nick on one 
edge of the filament itself. Values of the spectral 
transmission of the quartz window were also 
supplied by Mr. Van Dusen. 

In order to compute the spectral distribution 
of the radiation from the filament, it was neces- 
sary to obtain the true temperature of the fila- 
ment inside the bulb and to have a knowledge of 
the spectral emissivity of tungsten over the 
intended working range of wave-length (3000A— 
7000A). A compilation of spectral emissivities 
was furnished by Dr. W. W. Coblentz, retired 
head of the Radiation Laboratory of the National 
Bureau of Standards. The conversion from the 
brightness temperature to the true temperature 
is accomplished by the following formula, 
derived from Planck’s radiation formula: 


1 1 r 


T, Teo 








4343C2 





logio(ty: e), (1.1) 


where 


7,=true temperature of the filament 
Tsgo=brightness temperature of the filament outside the 

window as determined by the N.B.S. optical 
pyrometer 

\=mean effective wave-length in cm of the filter used 
in the N.B.S. optical pyrometer 

C,=the second constant in Planck’s radiation formula 

‘,=transmission of the quartz window at the above- 
mentioned wave-length 

€, =the spectral emissivity of tungsten at X. 


With Tg0=2614°K, \}=665010% cm, Co2= 
1.435, 4: =0.927, e =0.426, the true temperature 
of the filament is 2947°K. 

The spectral energy distribution outside the 
lamp was obtained by substitution of the above 
value of the true temperature in Planck’s radia- 
tion formula, and correction of the result for the 
spectral emissivities of tungsten and for the 
spectral transmission of the quartz window. 

When the image of the central portion of the 
filament was focussed on the slit of the spec- 
trograph, the resulting spectrum was a con- 
tinuous band with a narrower and darker band 
running down its length. This latter band was 
caused by the region of the filament at the nick 
where the filament was narrower, and its re- 
sistance and hence its temperature higher, than 
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regions of the filament above and below the nick. 
This narrow band was not wide enough to 
measure in the microphotometer, so it was 
necessary to use some other portion of the fila- 
ment where the temperature variation with 
length was small. Measurements of temperature 
with an optical pyrometer were made covering 
the 3” range from a point { inch above the nick 
to within § inch of the top of the filament. The 
true temperature fell only about seven degrees 
over this range, and at its center was 41° lower 
than that at the nick. The portion of this region 
used for projection onto the length of the slit 
(2 mm) of the spectrograph was about # inch. 
Hence the true temperature variation over this 
range was of the order of about 4°. From 
Planck’s radiation law one derives the percentage 
change in energy at a given wave-length for a 
change in true temperature to be 


(dE) / Ey) percent = 100(C2/AT?)dT. (1.2) 


For a true temperature of 2906°K the percentage 
change of Ey per degree is about 0.57 percent at 
3000A and 0.28 percent at 6000A. One may take 
0.4 percent as a mean. Hence, a change of 4° 
(true) corresponds to less than a 2 percent change 
in Ey. Furthermore, since the quartz lens used to 
focus the filament on the slit of the spectrograph 
does not give a perfectly sharp image, the effect 
is to smooth out variations in E, and produce a 
spectrum of uniform density across its width. 
Instead of recomputing the energy distribution 
for the adopted temperature of 2906°K, the one 
for 2947°K was used with the correction factor 


C=exp(C2/d)(71— T2/T1T2) 
=exp(C2/A)(AT/T?), (1.3) 


where AT=41° and T is the mean of 2947°K 
and 2906°K. It appears that the combined error 
from all sources in the computed energy varies 
from about 6 percent at 3000A to about 3.0 
percent at 6500A with an average for the middle 
of the spectrum of about 4.5 percent. 

One of the most difficult problems was the 
reduction of the filament brightness to a value 
closely approximating that of the fluorescent 
screen. This was necessary because of reciprocity 
failure of the photographic emulsion which 
demands that the exposure times on the standard 
lamp and the fluorescent screen be of quite the 
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Fic. 1. Schematic arrangement of reflection-type 
attenuator. 


same order of magnitude. A rough computation 
on the basis of available data indicated that the 
lamp might be brighter by a factor of the order 
of one million. Only by the combined use of three 
methods could the filament brightness be reduced 
to the point where the exposure times became 
nearly equal. Reduction by transmission filters 
had to be rejected as impractical, because trans- 
mission coefficients were not known with suf- 
ficient accuracy at the various wave-lengths. Of 
the following three methods used, one will be 
discussed in detail now and the second and third 
in later sections: 

1. The method of successive reflections. 

2. Rotating sectors. 

3. Stops in the collimator of the spectrograph. 


In the first method the brightness was reduced 
by successive reflections from four air-glass 
surfaces. The apparatus designed and 
structed for this purpose was called the “at- 
tenuator.”” The of four 
borosilicate crown glass prisms arranged as in 
Fig. 1. The arrows indicate the path of the light 
through the system. The system was enclosed in 
a brass box, and the prisms were carefully 
aligned so that parallelism and perpendicularity 
of the faces was assured. The angle of incidence 
of the central ray of the beam closely approxi- 
mated 45° at all of the four surfaces of reflection. 
At each surface a small fraction of the light was 
reflected, and the remainder was transmitted 
into the prism. In order that as little as possible 
of this light should be returned to the beam by 
internal reflections, all prism surfaces except 


con- 


attenuator consisted 
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the ones used for reflection were coated with a 
solution of propyl methacrylate mixed with 
lampblack, and the prisms were baked in an 
oven at 120°F for two hours. The result was a 
non-reflecting and absorbing coating which 
served to remove to a very high degree all light 
transmitted into the prism. 

According to Fresnel’s formulas the amplitude 
ratios of the reflected light are: 


sin(J — J’) 


sin( J+’) 


tan(J—1’) 
tan(J+T’) 





and 





for light polarized in, and, respectively, at right 


angles to, the plane of incidence. The intensity 
of the parallel component after four successive 
reflections becomes: 


sin(I—I’)\ ‘7? 
Rw) =| .707(———-) | 
sin(J+1’) 


1 /sin(J—1’)\§ 
= (= ) , (1.4) 
2\sin(J+ 1’) 


and similarly for the perpendicular component: 


1 /tan(J—J’)\$§ 
Ri) =-(——— ) (1.5) 
2\tan(I+/’) 





Computation shows that for the range of spec- 
trum to be studied, J; is about 9000 times fainter 
than J,, and hence the light may be said to be 
completely polarized. 

With J fixed at 45°, J’ depends upon n (Snell’s 
Law) and since n varies with the wave-length, R 
itself is a function of the wave-length. 

In order to determine the variation of m with 
wave-length, one of the prisms was mounted on 
a spectrometer and n was determined for each of 
the three mercury lines of wave-length 5460.7A, 
4358.3A, and 4046.6A. The resulting indices 
were, respectively, 1.5260, 1.5349, and 1.5389. 
Each determination had a standard deviation of 
unity in the fourth place of n—1. With these 
three values of the wave-length and their re- 
spective indices, the constants of the Hartmann 
formula were obtained: 


83.416 
2 = 1.5041 -+-——_____. 
A— 1648.0 


(1.6) 
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For several wave-lengths between 3000A and 
6500A, n was determined from (1.6), then J’ 
from Snell’s Law and finally the reduction 
factors from (1.4). The curve of reduction factor 
vs. wave-length is given in Fig. 2. These factors 
were applied to the spectral energy distribution 
curve of the filament outside the lamp and 
result in the distribution curve of the light as it 
leaves the attenuator. 

Attention must be given to the fact that the 
light from the attenuator is practically com- 
pletely polarized, while the light from the 
fluorescent screen is unpolarized. Thus the loss 
in the spectrograph may not be the same for both 
sources. In order to test this, two spectrograms 
of equal exposure times were taken of the lamp 
through the attenuator. For the first exposure, 
the plane of the normals to the reflecting surfaces 
was held perpendicular to the length of the slit 
and parallel for the second exposure. From the 
resulting densities of the spectra it appeared that 
the effective exposure was only about half as 
great in the first case as in the second. In order to 
avoid the difficulty the attenuator was mounted 
so that the plane of the normals made an angle 
of 45° with the slit and, therefore, as far as the 
spectrograph was concerned, it was receiving 
isotropic light. 

The errors in the computed reduction factors 
are evaluated as follows: 


1. As already stated nm has been determined with a 
standard deviation of one part in the fourth place of »—1. 
For the middle of the spectral range, this value results in 
an error of about 0.1 percent in the value of R— a negligible 
quantity. 

2. Since m was determined for only three wave-lengths, 
interpolation had to be used for shorter and longer wave- 
lengths using formula (1.6). The error in R to be expected, 
particularly in the short wave-length region where n is 
changing rapidly, is estimated to be of the order of 0.3 
percent in the neighborhood of 3200A and of 0.1 percent 
for wave-lengths longer than 4000A. 

3. The beam through the attenuator, as defined by the 
lens which focusses the filament on the slit and by the 
angular size of the collimator, has'an angular diameter of 
about two degrees, of which only the central ray will 
strike the reflecting surfaces at an angle of incidence of 45°. 
The limiting values at the first surface will be 44° and 46°. 
However, it will be seen that (by design) if J is 46° at the 
first surface, it will have the same value at the second, 
44° at the third, and 44° at the last surface. Calculation 
shows that the error made in R by assuming J=45° for 
the whole beam is about 0.2 percent and thus is negligible. 
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Fic. 2. Reduction factor of attenuator as a function 
of wave-length. 


4. Particular care was used in the alignment of the 
attenuator so that the resulting error in R is probably less 
than 0.4 percent. 


The total error is probably not in excess of 
1 percent. 

The spectrograph was an all quartz optics, low 
dispersion instrument. The arrangement was 
such that the plate holder could be moved at 
right angles to the spectrum, and thus several 
spectra could be obtained on one plate. The 
length of spectrum from 3000A to 6500A was 
about } inch. The spectrograph was mounted on 
a turntable which permitted the slit to be turned 
to either the standard lamp, the mercury lamp 
(for a comparison spectrum), or the fluorescent 
screen. By means of lead sheets, the whole 
instrument was very well shielded from the 
direct or scattered radiation of the x-ray tube. 
On the turntable was mounted a motor to turn 
cardboard sectors, used to reduce the lamp 
brightness by a given amount, and also a mount 
for the quartz lens to focus the lamp filament on 
the slit. Inside the prism case was an arrange- 
ment to hold stops for the further control of the 
lamp brightness. The length of the slit was 2 mm 
and the width 0.2 mm. 

The fluorescent material, a sample of Eastman 
Ultra Speed back screen, consisted of a layer of 
CaWO, covered by a film of lacquer and mounted 
on a sheet of white cardboard four inches square. 
It was centered over a two inch square hole in a 
sheet of } inch lead and the whole unit was 
mounted over the open end of the protective cone 
of the x-ray tube. 
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The Experimental Procedure 


Generally speaking, our procedure was a com- 
parison, at a number of wave-lengths, between 
the brightness of the screen and the standard 
lamp filament as seen through the attenuator. 
Given the magnitude of the latter and its spectral 
distribution, it was possible to compute the light 
output of the screen on an absolute scale. 

One of the most difficult problems, the reduc- 
tion of the brightness of the standard lamp to a 
level comparable to that of the fluorescent 
screen, was accomplished by the use of the 
attenuator (already described), rotating sectors, 
and collimator stops. This permitted the use of 
equal exposure times for the lamp and the 
fluorescent and thereby avoided the 
danger of reciprocity failure. In order to provide 
a calibration on each plate, five exposures were 
made using five cardboard sector disks, eight 
inches in diameter; the first, A, providing a 
reduction factor 12.5 and the others, B, C, D, 
and E approximately 25, 50, 100, and 200. When 
using these disks, one must avoid the inter- 
mittency effect because of which the photo- 
graphic density is a measure of the total light 
impinging only under certain conditions. A con- 
tinuous exposure and an intermittent one of the 
same total energy will produce the same density 


screen 


on the photographic plate, provided the total 
time of exposure is the same in both cases, and 
the frequency of interruption is above a certain 
critical rate. The Eastman Kodak Company’s 
publication, ‘‘Photographic Plates for use in 
Spectroscopy and Astronomy,” states that this 
critical rate varies from about 100 flashes per 
second for a one-second exposure to a few. flashes 
per second for a one hour exposure. The motor 
used to drive the sectors operated at 2850 r.p.m. 
and since no exposure was less than twenty 
minutes, no error caused by the intermittency 
effect was incurred. 

The final adjustment of lamp brightness was 
obtained by placing one of several stops between 
‘the prism and the collimator lens at a position 
close to the latter. The first stop defined the 
beam leaving the collimator lens to the extent 
that all of the beam entered the prism and passed 
through the camera lens to the photographic 
plate. This stop was used when photographing 
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the spectrum of the fluorescent screen, and its 
reduction factor was designated as unity. Six 
more stops were made with reduction factors 
varying from 2 to 37. With this last stop, the 
total reduction of the lamp brightness, including 
the average of that of the attenuator and that of 
sector A, reached a factor of nine million. This 
much was needed when working in the region 
around 6000A where the screen fluorescence is 
weak, and the lamp output high. 

All photographic plates used were coated with 
Eastman 103-F emulsion. On each plate eight 
exposures were made; one from the fluorescent 
screen, five from the standard lamp using sector 
disks A to E, and two exposures from a mercury 
lamp to obtain wave-length scales at the top and 
bottom of a plate. 

The densities of the spectra were compared in 
an Arl Dietert microphotometer. At each position 
the density was determined for the fluorescent 
spectrum and for each of the five comparison 
spectra. 

Since different portions of the spectrum re- 
quired widely varying exposure times to give 
densities on the linear portion of the emulsions, 
it was necessary to use several plates to cover the 
complete spectrum. For each wave-length the 
transparencies of the five comparison spectra 
were plotted on semi-logarithmic paper against 
the reduction factors of the sectors, taking that 
for sector A as unity. From this curve was deter- 
mined the factor by which the fluorescent spec- 
trum was fainter than the spectrum belonging to 
sector A. 

This factor was combined with reduction 
factors for the collimator stop, the factor 12.5 
for the sector disk A and the factor for the 
attenuator, to give a total reduction ratio 
between intensities at corresponding wave- 
lengths from the standard lamp and from the 
CaWO, screen. A total of seven plates were 
measured. 


Experimental Results 


Near the end of the investigation it was found 
that the portion of the CaWO, screen which had 
been exposed to x-rays was somewhat fainter, 
showing that it had suffered a certain degree of 
“fatigue.”” Measurements indicated that the 
exposed portion of the screen was about 8 percent 
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fainter than the outer portion that had been 
protected from x-radiation by the lead. 

In order to deduce from the measurements on 
a fatigued screen the spectrum of a fresh screen, 
measurements on each plate were corrected by 
an amount proportional to the total x-ray ex- 
posure time of the screen up to the time when 
the plate in question was taken. 

The total time of exposure of the screen to 
X-rays was about six and one-half hours. This 
includes exposure times for test plates and 
several others that for one reason or another 
were not used. The results from the individual 
plates are given collectively in Fig. 3. The 
ordinates are 10-° ergs/cm® sec. A unit. The 
assumption is made that both the lamp filament 
and fluorescent screen radiate 
Lambert’s cosine law, and a 
departures will be added later. 

The curve given in Fig. 3 is not final, but must 
be corrected for slit broadening. The effect is to 
pull in both ends of the curve to a certain extent 
without changing the total area. Integration of 
this curve yields the result that the total energy 
output of the screen was 2.54 ergs/cm?/sec. 

A correction factor of 0.92 must be applied to 
this value to correct for the departure from 
Lambert’s law in the angular distribution of the 
light from the screen. The experimental deter- 
mination of this factor is described in Section Ib. 
The total energy thus reduces to 2.34 ergs/square 
centimeter per second, for a measured x-ray 
intensity of 0.149 r-units per second incident on 
the screen. The final corrected spectral curve, 
Fig. 4, is given in terms of light energy generated 
in the screen per r-unit absorbed. The method for 
measuring the light losses is given in Section 2. 
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Fic. 3. Spectrum of CaWO, as obtained from seven plates 
(uncorrected). 
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Fic. 4. Fluorescent spectrum of CaWO, excited by 


80 Kv constant potential x-rays. Filter used was } mm 


copper+1 mm aluminum. Screen contained .109 gm/cm? 
of CaWQ,. 


An estimate of the accuracy of this deter- 
mination involves a consideration of such sources 
of error as the temperature of the lamp, the 
spectral emissivities of tungsten, the brightness 
reduction factors, and errors of the photographic 
emulsion. It appears that the final total energy 
has been determined to within an accuracy of 
about 12 percent. 


Ib. ANGULAR DISTRIBUTION OF THE LIGHT 


In order to calculate the total light emitted 
from the screen from the measured brightness 
in a specific direction, it is necessary to know 
the angular distribution of the emission from the 
screen. In the case of the filament, this factor was 
taken care of by using the normal emissivity of 
tungsten for calculation. For the screen, measure- 
ment of the angular distribution was made by 
means of a photo-multiplier cell illuminated 
through a collimator consisting of a simple lens 
and a small diaphragm stop at its focal point. 
The whole assembly was capable of rotation 
about an axis lying in the plane of the fluorescent 
screen. 

If the emission of light followed a Lambertian 
(cosine) distribution, the light received by the 
photocell would be independent of angle. Ac- 
tually a falling-off at large angles was observed 
as shown in Fig. 5. Integration of this curve 
combined with suitable angular variables yields 
a constant which relates the total light emitted 
from the screen to that which would have been 
emitted from a Lambertian source with the same 
normal intensity. This factor varies slightly with 
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Fic. 5. Angular distribution of light from CaWQO, screen 
relative to a Lambertian surface. 


the wave-length of the exciting x-rays, since the 
latter affects the distribution of light-generating 
centers within the screen. Under the conditions 
of the spectroscopic experiment the CaWQ, 
screen had a factor of 0.920. 


II. TRANSMISSION OF LIGHT IN 
DIFFUSE SCREENS 


The light produced by the x-rays in the 
interior of the fluorescent layer is eventually 
transmitted out of the screen by means of a 
series of reflections and refractions from indi- 
vidual crystals. Part of the light is lost by ab- 
sorption in the crystals and in the binder, and 
part finds its way out of the back surface of the 
layer. Some of the latter may be recovered by 
means of a reflecting backing, a white mounting 
card being commonly used. 

Since the particles involved are small com- 
pared to the area of interest, and the lateral 
dimensions of the screen are large compared to 
the thickness, the problem can be treated as if 
the material were homogeneous and as if the 
light were traveling in one dimension only, i.e., 
perpendicular to the screen. 

In the following discussion the words “light 
forward” and “light backward” refer to direc- 
tions in the same sense as that of the x-rays. 
Thus with the screen located between the x-ray 
_ tube and the observer, the light forward is that 
seen by the observer and comes from what will 
be termed the front of the screen; the light 
backward from the back of the screen is directed 
toward the x-ray tube. 

The first step in the analysis is to derive the 
reflection and transmission coefficients of a 
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diffuse layer as a function of its thickness. A 
layer of thickness x (Fig. 6) has by definition a 
reflection coefficient R(x). If a layer of thickness 
dx is added to this, a new reflection coefficient 
will result. The reflection coefficient of the 
elementary layer itself will to the first order be 
proportional to its thickness. Let this reflection 
coefficient be kdx. Similarly, the amount not 
transmitted will be proportional to the thickness 
to the first order; let this coefficient be adx. It 
will be evident that a>k, being equal to k if 
there is no true absorption,of light in the layer. 
(In Fig. 6 the layer dx is shown separated from 
the layer x for convenience in tracing the course 
of the light. It is to be understood that in this 
and later derivations the layers are actually con- 
tinuous.) Light is incident on the combination 
from the left. A fraction kdx will be reflected 
from the first layer, a portion 1—adx will be 
transmitted to the second layer. The amount 
reflected from the second layer is (1—adx)R, of 
which a quantity (1—adx)’R is transmitted out 
to the right and adds to the initially reflected 
light. The portion kdx(1—adx)R which is 
returned undergoes an exactly similar set of 
steps. Summing up the infinite series of rapidly 
diminishing contributions: 

kdx+ R(1 —adx)*(1+ Rkdx+ R?(Rkdx)*+ ---) 
R(1—adx)? 
= kdx+—-—————_ = R(x+dx), (2.1) 

1— Rkdx 


the new reflection coefficient R(x+dx) is ob- 
tained. Subtracting R from both sides, and 
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Fic. 6. Schematic representation of light flow for derivation 
of reflection coefficient formula. 
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neglecting differentials of higher order than one, 
we obtain 


R(x+dx) — R(x) =k(1—2(a/k)R+R®)dx 





or 
dR 
=k(1—2(a/k)R+R?). 
dx 
Thus: 
dR 
kx = [ ' 
J [1—2(a/k)R+R?] 
1 R-—a/k 
kx = ——————- tanh“! ——+const., 


L(a/k)?—1]! L(a/k)?—1]! 


(2.2) 
a 
R= di k)?—1 |! tanh((a?—k?)'x+const.). 


The constant is determined by setting R=0 
when x=0, and by further reduction we obtain 


(Rk, p) sinhpx 





= , (2.3) 
coshpx+(a/p) sinhpx 
where 


p=(a?—k?*)}. 


By similar consideration of the multiple reflec- 
tions in the transmission case, and by use of 
formula 2.3, the expression for the transmission 
coefficient is found to be 


1 





T (2.4) 


~ coshpx+ (a p) sinhpx 


where p has the same meaning as before. 

The constant p is an important characteristic 
of the material and represents the degree of ab- 
sorption. For very large x, 2.4 reduces to 


7 =constant e~’* 


the usual exponential function, and p plays the 
role of a linear absorption coefficient. At large x, 
R approaches a constant, and only after this 
condition is reached will T be exponential. The 
physical reason for this is clear; if a thin layer 
is placed behind another layer, the incident light 
on the new layer is enhanced by the multiple 
reflections between the two, and only if the 
original layer is so thick as to have attained the 
limiting reflection coefficient will the trans- 
mission be exponential. The form of these func- 
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Fic. 7. Reflection and transmission coefficients as a 
function of thickness (x) of diffusing material. 


tions is shown in Fig. 7 for values of k and a 
typical of CaWo,. The values of k& and p (or k 
and a) suffice to describe completely the light 
yield from the front and back of the screen, as 
well as the light absorbed, if the thickness and 
distribution of light sources within the screen 
are given. These coefficients may be experi- 
mentally determined by measuring R and T. 


Determination of the Optical Constants 


Values for p, k and acould be obtained in prin- 
ciple by measuring T or R as a function of thick- 
ness, by virtue of Eqs. (2.3) and (2.4). In carrying 
out such an experiment one must make sure, 
however, that the amount of light from the 
source and the action of the detector for the 
light are not affected by introduction of the 
sample. For example, a ground-glass diffuser is 
not a suitable source because it is itself reflecting; 
introduction of the sample will enhance the 
incident light intensity to a value higher than 
that measured without the sample. For the 
same reason the detector must be non-reflecting; 
a photographic plate, for instance, is not suf- 
ficiently absorbing to avoid errors due to mul- 
tiple reflections. 

Since the procurement of such ideal light 
sources and detectors constitutes a serious ex- 
perimental difficulty, it seemed preferable to 
eliminate the error altogether by means of a 
method which does not require measurement of 
light intensity in the absence of the sample. A 
tungsten bulb, suitably housed, was used to 
illuminate one side of the layer; the brightnesses 
of the two sides were measured by means of a 


Macbeth Illuminometer and the ratio of the 


’ 
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brightnesses was taken to be R/T. This was 
repeated with samples of varying thickness. 

By definition, R and T are the results of a 
scattering process which does not alter the 
angular distribution of the light. Therefore, 
Eqs. (2.3) and (2.4) must not be applied 
indiscriminately to the present case where a 
nearly parallel light beam is incident on the 
front of the sample and, upon entering, is 
gradually scattered in all directions, approaching 
an equilibrium distribution in direction only as 
a limiting case. The angular distribution of the 
light thus varies with depth of penetration, 
whereas Eqs. (2.3) and (2.4) were derived 
from an assumed stationary equilibrium for the 
angular distribution; so the latter is presumed 
to remain unaffected by further scattering and 
not to vary with depth. 

While this poses a very serious difficulty, 
generally speaking, for the measurement of 
arbitrary samples, the calcium tungstate screen 
has the fortunate quality that its scattering 
power outweighs its absorbing power by an order 
of magnitude; consequently determination is 
possible by examining sufficiently thick samples. 
If the conversion from parallel light to light 
with completely random orientation takes place 
in a relatively thin front layer, further re-dis- 
tribution will proceed in accordance with the 
underlying assumptions of Eqs. (2.3) and 
(2.4). 

An empirical criterion of this condition is 
afforded by inspection of the spatial distribution 
of the light emerging from the front and back 
of the screen. For a very thin layer it is clear 
that while the reflected light has a nearly random 
distribution, the transmitted light will contain 
large portions of unscattered light; thus a dif- 
ference in the angular distribution of the light 
from the front and the back indicates that 
equilibrium has not yet’ been reached. In the 
thinnest sample used a rough check of the 
angular distributions showed they were alike and 
nearly Lambertian. Furthermore, the reflectivity 
of the thinnest sample was calculated to be 85 
percent. Thus addition of further layers can 
change the reflected light by 10 percent at most 
and the brightness of the front side is always a 
good measure of the incident light quantity. For 
this reason, even though the values of k and a 
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may be somewhat in error due to our use of a 
parallel beam, this does not apply to the value 
for p, which can be computed from the logarith- 
mic derivative d log7/dx for large thicknesses x. 

Samples were prepared by soaking an Eastman 
Ultra-Speed front screen in alcohol and thus 
removing its fluorescent layer. This layer con- 
tained 0.0445 gm/cm? of CaWQ, per cm’. It is 
self-supporting because the internal binder is not 
soluble in alcohol, whereas the cement used to 
mount it on the card is soluble. One-inch squares 
were cut from the layer while still wet, and were 
superimposed to make up samples of 1, 2, and 
4 layers. These were cemented over holes in 
metal sheets and allowed to dry. 

It might be supposed that the protective 
coating on the layers (a film of cellulose plastic) 
would violate the assumption of homogeneity 
of the material. However, consideration of the 
multiple reflections between the crystalline 
material boundary and film-air interface show 
that for the index of refraction of this film the 
effect is completely negligible for the thicknesses 
used. 

Four points on the R/T curve were deter- 
mined, at zero, 1, 2, and 4 layers. Values of R/ T 
are shown in Table I. From Eggs. (2.3) and (2.4) 
we see R/T should be of the form &/p sinhpx. 
Values of k and p were chosen to fit the ex- 
perimental points, namely, k=5.55/layer and 
p=0.30/layer. The third column gives the cal- 
culated points for comparison. Since a? =k?+ p’, 
a is also 5.55/layer within the experimental error, 
with the understanding of course that 


(a? — k*)}=0.30/layer. 


From the fit obtained, the values of k, a, and p 
are estimated to have probable errors of less 
than 3 percent. In the following treatment these 
values have been used, and the thickness x is 
always measured in terms of one layer as the 
dimensional unit. 








TABLE I. 
x (layers) R/T (exp.) 18.5 sinh.30x 
0 0 0 
1 5.43+.05 5.63 
2 11.8 +.13 11.8 
4 27.6 +.4 27.9 


JOURNAL OF APPLIED PHYSICS 











~_ — 


~ = =e Le & 








After these measurements had been evaluated 
and the results had been used to compute the 
conversion efficiency and the curves of light 
yields, it was decided to repeat the measurements 
of the optical constants of the screen using a 
diffuse source for the incident light. The sample 
was mounted over an aperture in an Ulbricht 
photometric sphere arranged so that the bright- 
ness of both sides of the sample could be ob- 
served with a Macbeth Illuminometer. These 
measurements gave a good fit for the ratio 
R/T=k/psinhpx, using k=6.0 and p=0.33. 
Adoption of these new values would change the 
curves of light yields which are presented below. 
Since these are meant to be illustrative, and 
quantitative only pending measurement of the 
optical constants of an individual sample, it did 
not seem worth while to recompute them. How- 
ever, the new constants were 
putation of conversion efficiency. 


used for com- 


The Light Distribution in the Screen 


Having derived expressions for transmission 
and reflection of a diffuse layer, one is in a 
position to compute, for a given distribution of 
light generation throughout the screen, the frac- 
tional amounts emerging from the screen in both 
the forward and backward directions. 

Since the x-ray intensity decreases exponen- 
tially with depth of penetration, light production 
also follows an exponential law, being highest at 
the back where the x-rays are incident. Consider 
the element of thickness dx located at a distance 
x from the back of the screen whose thickness is «€ 
(Fig. 8). This element is between two layers, of 
finite thickness x (layer 1) and e—x (layer 2), 
respectively. Of the light dL generated in the 
element, 3 starts out toward layer 2, while 3 
starts toward layer 1. These are subject to mul- 
tiple reflections as shown in Fig. 8, the lower half 
being concerned with the light starting toward 
2 and the upper with that starting toward 1. 
The fraction f out the front due to dL is given 
by the sum of the contributions as before, and 


1 7:(1+R)) 
=— ——__—, (2.5) 
2 1—RiR2 


The coefficients R;, Re and 7:2 are, of course, 
functions of x and e—x as in formulas (2.3) and 
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Fic. 8. Schematic representation of light flow from a 
luminous element within the screen. 


(2.4). It is now necessary to multiply f by the 
light produced in the element dx and integrate 
over x from 0 to ¢€ to obtain the total light out. 
If the absorption coefficient of the x-rays is yu, 
we have by definition that the x-rays absorbed 
in a layer dx =Iydx where I is the intensity of 
x-rays at the point x. Let p =conversion efficiency 
of x-rays into light; the light produced in dx is 
then pludx. But IJ=IJce-**, where Jo is the 
incident x-ray intensity, and dL = pylpe~**dx (2.6) 
The forward light yield is then: 





dx. (2.7) 


pul o “ T e-2) (1+ Riz) 
L=— f he . 
0 (1 — Riz) R¢e-z)) 


2 


Before performing this integration we shall find 
it very convenient to derive a relation between 
the transmission coefficient of the screen as a 
whole (7.) and the coefficients of its parts. In 
Fig. 9 is shown the schematic diagram from 
which we find: 


T.=1,T2(1+RiR2+ (RiR2)?+---), 





(2.8) 
T.=T,T2/(1—RiR2). 
Substituting in formula (2.5) we get 
T.f1 R, 
j-(—+- ). (2.9) 
2A\T, Ti 
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7, is a constant independent of x, and from 2.3 
and 2.4 we have 1/7,=coshpx+(a/p) sinhpx, 
R,/T,=(k/ p) sinhpx. Formula 2.7 then becomes: 


pulol. P 
L= fe pe 
2 0 


at+k 
x (coshpx+ — sinhps Jas (2.10) 
p 


In preference to an expression for light out itself, 
we shall divide LZ by the total light generated to 
obtain the fraction Fr which is transmitted out 


— wi (utatk)—e™[(ut+at+k) coshpet (Lu(a+k)/p]+p) sinhpe]} 


the front of the screen. The total light produced 
is: 


pI o(1 —e#*) 
and 


bu 


(1—e-**)[coshpe+(a/p) sinhpe ] 


*f at+k 
x | e «( coshps+ — sinhpx Jas. (2.11) 
0 Pp 


Fy = 





Carrying out the integration we obtain the final 
expression : 





PF 


(2.12) 


2(1 —e-**) (wu? — p”) (coshpe+ (a, p) sinhpe) 


In an exactly similar manner the formula for the fraction out the back of the screen is found; 


ul (u—a—k)(coshpe—e~“*) + ([u(at+k) /p]—p) sinhpe} 





Fp=- 


(2.13) 


2(1—e-**)(u? — p?) (coshpe+(a@, p) sinhpe) 


Graphs of these functions are given in Figs. 
10 and 11 for various values of » and e, using 
k=5.55/layer and p=0.30/layer. Values of the 
thickness € are in terms of one-layer units, and 
u, the x-ray absorption coefficient is per laver. It 
is apparent that for zero thickness of screen half 
of the light is transmitted in each direction 
regardless of uw. At the limit u—0 light is produced 
uniformly throughout the screen, and for thicker 
screens equal amounts of somewhat less than 
50 percent are emitted in each direction, the 
remainder being absorbed in the screen. For 


thick screens and high y, the fraction emitted. 


toward the back is much more than 50 percent, 
while the light forward is much less. For p= ~, 
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Fic. 9. Light flow diagram for derivation of formula 
relating partial and total transmission coefficients. 
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all of the light is produced in a layer of zero 
thickness at the back. Half of the light is im- 
mediately emitted backwards and a fraction 
R(e) of the remaining half is reflected in addition. 
The fraction out the front is the transmission 
coefficient times the half originally directed 
toward the screen. Thus all the curves of Figs. 10 
and 11 approach different limits as p> =. 

While these curves show wide variations in 
yields, actually the ranges of uw and e€ usually 
encountered in practice are such that the frac- 
tions forward and back are fairly close to 50 
percent. At 80 kv constant potential, yu is ap- 
proximately 0.25 per layer; for a screen of «=1 
(Eastman Ultra-Speed front) the fraction for- 
ward is 47.6 percent, the fraction back 51.2 per- 
cent. At higher voltages little error is introduced 
by assuming 50 percent each way unless the 
screen is abnormally thick. At the lower voltages, 
u rises rapidly, attaining a value of 2.2 at 20 kv 
and 3 at 17.4 kv. Here the fractions out the front 
and back are markedly affected by the operating 
voltage. 


Effect of Mounting Card 


In practice a screen is never used without a 
backing of some sort. Usually this is made white, 
with the intent of recovering the light emitted 
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Fic. 10. Fraction of light generated which escapes from 
the front of the screen as a function of screen thickness (e€) 
and x-ray absorption coefficient (z). 


in the backward direction. The efficacy of this 
procedure will now be examined. Figure 12 
shows schematically a layer of fluorescent ma- 
terial backed by a card of reflection coefficient Ro. 

The ratio of the light returned to the light 
emitted backward from the screen is defined as E. 


E=RoT(1+RRo+(RRo)?+::-, 
RoT 
~1—RRo 


(2.14) 
E 


By way of example, if a screen of thickness «= 1, 
R=0.844 and T=0.15, is backed by a card of 
reflection coefficient R»p=0.80, one finds E = 0.37, 
so of all light which was previously lost out the 
back of the screen, only 37 percent is recovered 
by using a backing card of 80 percent reflec- 
tivity. This is the result of the ‘“‘whiteness”’ of the 
calcium tungstate; if the card is not of superior 
reflecting power, it is of small value. Some of this 
loss is in the layer itself, but the major portion 
goes into the card. The card upon which the 
Eastman screen was mounted had a diffuse 
reflection coefficient of 0.84. Table II illustrates 
the role played by this mounting card in the 
light losses. These data were computed for 80 kv, 
using the constants k=6.0 and p=0.33. At 
higher voltages very little change takes place in 
the light distribution, and the tabulated values 
will remain essentially constant. At lower 
voltages the losses will increase considerably, 
reaching a limit at ».=0 of 52 percent for the thin 
screen and 74 percent for the thick screen. 
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Fic. 11. Fraction of light generated which escapes from 
the back of the screen as a function of screen thickness (e) 
and x-ray absorption coefficient (x). 


The increase in light to be expected from the 
use of a card has been plotted in Fig. 13. At 
first glance these numbers may seem quite high; 
however, it must be remembered that in the vast 
majority of cases yp is 0.5 or smaller, so that the 
increase will seldom exceed 40 percent. 

Two experimental checks of the above for- 
mulas were made. In the first the light from a 
fluorescent screen was measured with and 
without a backing card. In the second, two 
screens of different thickness were measured and 
the expected relative light output calculated on 
the basis of the above formulas. Both tests 
checked calculated values within 3 percent. 

So far we have not discussed the variation 
with thickness of the amount of x-rays absorbed. 
This, of course, is the most important factor in 
determining light output for thin screens. With 
increasing screen thickness the x-rays are more 
and more nearly completely absorbed; in addi- 
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Fic. 12. Light flow diagram for derivation of effect of 
reflecting card. 
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Fic. 13. Increase in light to be gained by using a backing 
card with 80 percent reflectance, as a function of screen 
thickness (€) and x-ray absorption coefficient (x). 


tion, the effectiveness of the card backing 
decreases and losses in the phosphor and binder 
begin to gain importance. Thus the light output 
reaches a maximum at a thickness which will 
depend upon the absorption coefficient (or 
voltage) of the x-rays. In Figs. 14 and 15 is 
shown the utilization factor U, defined as the 
ratio of the light obtained from the front of the 
screen to the total light which would have been 
generated had all the x-rays been absorbed. The 
curves of Fig. 14 were calculated for a screen 
without mounting card, those of Fig. 15 with a 
mounting card having 100 percent reflection. It 
will be seen that in practical cases the utilization 
is seldom over 30 percent. 

In the case of a ‘‘back”’ screen one is interested 
in the light out the back, i.e., toward the x-ray 
‘tube. This quantity increases with the screen 
thickness, but the gain is very slight beyond a 


TABLE II. Light losses in two typical x-ray 
intensifying screens. 





Thin screen Thick screen 





(1 layer) 2.45 layer 
.045 gm/cm? .109 gm/cm? 
CaWwO,. CaWwO, 
Light lost in phosphor and binder 2.2% 9.6% 
Light lost in mounting card 28.7% 36.0% 
54.4% 


Light yield (total) 69.1% 
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certain point. A thick layer of CaWO, is much 
whiter than any available cardboard; for a 
simple example we can just set the thickness of 
the back screen at infinity and inquire how the 
utilization factor varies with uy. If ¢€ is infinite, 
Fy, in formula 2.13 is the utilization factor, since 
all the x-rays are absorbed. The formula reduces 
to: 


Us=[(ptatk)/2(p+a) ](u/utp). 


The factor is plotted in Fig. 16 using values of 
p and a for CaWQ, as before. It is seen that even 
with an infinitely thick screen the utilization is 
not very high for uw less than 0.3, i.e., for x-ray 
voltages above 80 kv. 

In the design of screens for radiographic pur- 
poses an all-important consideration is the re- 
solving power of the screen-film combination. In 
general, the thickness of the screen is largely 
determined by the stringency of this require- 
ment. No treatment power is 
attempted here, because the primary purpose of 
the investigation is to measure the true efficiency 
of the CaWQ, screen. The curves in the preceding 
portion show accurately the light produced as 
the screen thickness and x-ray absorption coef- 
ficient vary, but they form a basis for selection 
of screen thickness only on the assumption that 
resolving power is to be ignored. For this reason 
it is of littke value to extend the treatment to 


(2.14) 
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Fic. 14. X-ray utilization factors for a screen with a 
non-reflecting backing. 
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other cases, such as a front and back screen 
combination. The information presented here is 
sufficient to calculate the light yield of any such 
combination, given the absorption coefficient, 
thicknesses, and mounting card reflectances of 
the screens used. 

It should be pointed out here that if the film 
used with a screen is itself reflecting, all of the 
light will not be absorbed in the film. Measured 
white-light reflectances of x-ray screen film are 
high enough to cause errors of 8 percent with an 
Ultra-Speed front screen; however, it is quite 
possible that the reflectance is different for the 
blue and ultraviolet light of the screens. 


Ill. EFFICIENCY OF FLUORESCENCE 


To obtain the x-ray energy absorbed in the 
screen, it is necessary to have a knowledge of the 
energy spectrum of the incident radiation and of 
the absorption of energy from this spectrum by 
the screen. 

Information as to the exact continuous spec- 
trum of x-ray radiation from a thick target is not 
available in the literature. Curves originally 
published by Ulrey' have been reprinted in all 
the texts for 25 years in spite of the fact that 
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Fic. 15. X-ray utilization factors for a screen with a 
perfectly reflecting back. 


'C. T. Ulrey, Phys. Rev. 11, 401 (1918). 
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Fic. 16. X-ray utilization factor for a very thick screen 
used as a back screen, plotted against pu. 


they were known to be in serious error at the 
time of publication. Our use of these curves in 
combination with the known absorption coef- 
ficients of CaWOQ, led to obvious and large dis- 
crepancies with the measured absorption of the 
screen, and an independent determination of the 
x-ray spectrum seemed desirable. Figure 17 
shows the energy distribution from the x-ray 
tube used by us, as determined by Mr. R. L. 
Longini with a method to be published by him 
at a later date. This spectrum applies to the tube 
run at 80 kvc with } mm copper and 1 mm 
aluminum filter. On the same figure is given the 
curve of energy absorbed by the screen, calcu- 
lated from the known total absorption coef- 
ficients. The integrated areas under these curves, 
when converted to ionizing power, agreed with 
our measured absorption value. 

All x-radiation which leaves the screen is con- 
sidered not to be part of the input for the purpose 
of our computation of conversion efficiency. Our 
measurement of the undeflected x-ray beam 
leaving the screen fails to include the scattered 
and fluorescent x-radiation which leaves under 
large angles. The former can be shown to be 
negligible at the voltages used; the latter is 
important only for energies above the K-absorp- 
tion edge of tungsten. Of the 299 ergs absorbed 
in the region to the left of the K-edge in Fig. 17, 
absorption by K-electron ejection includes 190 
ergs which will be re-emitted as K-radiation. 
Of this, about 5 percent is converted into Auger 
electrons and 26.5 percent of the remainder is 
reabsorbed in our particular thickness of screen, 
so that 133 ergs escaped as x-rays. Thus the 
total x-ray input must be reduced from 714 ergs 
to 581 ergs. 
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Fic. 17. Spectral energy distributions of x-rays incident 
on the screen and x-rays absorbed by the screen. 


We are now in a position to calculate the 
efficiency of calcium tungstate. From Part I we 
have the light out of one side as 2.34 ergs/cm?/ 
sec with an incident radiation of 8.92 r/min. 


At 581 ergs/cm?’ r we have then 86.2 ergs/cm?/sec 
absorbed in the screen. Thus the ‘“‘gross’’ effi- 
ciency is 2.72 percent. From Part II, however, 
we find that the screen emits only 54.5 percent 
of the light generated, and the intrinsic efficiency 
of the calcium tungstate is finally found to be 
5.0 percent. 

This efficiency may be expected to vary only 
slightly with the wave-length of the x-radiation 
involved, since one electron ejected by an 
x-photon excites many fluorescent centers and 
dissipates its energy gradually. 

It is of interest to compute the number of 
fluorescent centers activated by a single ab- 
sorbed x-photon. Taking the average wave- 
length of the absorbed x-rays as 0.22 angstroms, 
and the average wave-length of the light pro- 
duced as 4400 angstroms, we see that energeti- 
cally it would be possible to get 20,000 light 
quanta from one x-photon. Since the material is 
only 5 percent efficient, however, on the average 
each x-photon absorbed produces 1000 light 
quanta. 
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Pressure and Oil Flow in Oil-filled Cables at Load Variations 


Hans Lottrup KNUDSEN 
Holstebro, Denmark 


(Received September 23, 1946) 


This paper presents a method for calculating pressure and oil flow in a section of an oil-filled 
cable at load variations with due regard to both the hydraulic resistance of the oil canals and 
the elasticity of the sheath and the pressure armor, if any, and the compressibility of the oil. 
The present article restricts itself to cable sections terminated with oil reservoirs in which the 
pressure is constant; in a forthcoming paper this restriction will be dropped. Formulas are 
given for pressure, oil current, and amount of oil which flows through a cross section of the 
cable for an arbitrarily prescribed oil expansion per unit length. These formulas are applied 
to the following oil expansion functions: So-1, d/dt(No-1), e~“7-1, and 1/+/t-1, which repre- 
sent the oil expansion or terms occurring in the series giving the oil expansion in some important 
cases. Generalized functions are plotted for the pressure at the midpoint, the oil current at 
the endpoints, and the amount of oil which passes the endpoints for the above mentioned oil 
expansion functions. Using these generalized functions, the calculation of pressure and oil 
flow, with the simplifying assumptions made here, is reduced to elementary operations. Nu- 
merical examples illustrate the use of the given formulas and curves. 





LIST OF SYMBOLS P(t) cm"! W Effect transferred from con- 
. ; ductors to insulation per unit 

a cm? Area of oil canal. length 

A cm" Oil volume per unit length. °C-' cm™* W sec. Specific heat of insulation. 

7 kg"? cm* Capacitance per unit length. kg cm~ sec. Hydraulic resistance in cable 

C1 kg! cm‘ Capacitance per unit length per unit length 
corresponding to the elasticity R kg cm~ sec. Hydraulic resistance of entire 
of sheath and pressure armor. ont 

™ kg? cm* Capacitance per unit length S(t) cm? sec.~! Oil expansion per unit time 
corresponding to the com- per unit length 

: pressibility of oil. . So cm? sec.~! Constant factor in expression 

Cc kg! cm5 Capacitance of the entire for oil expansion 
cab le. : ; t sec. Time. 

d cm Diameter of oil canal in ,. ; ae 
incl bl 7 sec. Time constant occurring in 
ee the expression for oil expan- 

D cm Diameter of cable under dais 
heath. : 

: s — ik V(x,t) kg cm? Pressure. 

E kg cm™? Modulus of elasticity for the ty 
heath ial. ‘ . : hi 
rn ar eebeenige ap B —" Volumetric expansion coeffi- | 

f cm Perimeter of the conductors ‘ : 

. @ ae cient of oil. 
in which short-circuit current ; cau : ‘ 
Rete B Cc Volumetric expansion coeff- 
' a ste — cient of insulation. 
F(x) 1 Function defined in (54). - a f " 
“ ; v gem Density of oil. 
I(x,t) em sec.*! Current of oil. : 
; " eee a : 6 cm Sheath thickness. 

k kg~! cm? Specific compressibility of oil. ‘ ze eg? ’ 

" ‘ 4 ; : d(x, t) c Temperature in insulation, 

Ko cm? sec. Constant factor in an expres- ; 

: , Jo i; Temperature of conductors 
sion for the oil expansion. x 
ot A : after short circuit. 
l cm Arbitrarily defined length in , etnipedy: : 
: v centistokes Viscosity of oil. 
cable cross section. 

L cm Length of cable. 

M(x,t) cm Amount of oil which flows PRESSURE AND FLOW QUALITATIVELY 
through a cross section of the REGARDED 

; cable. ; N the calculation of oil-filled cable systems it 

No cm? Instantaneous oil expansion _ : 
oot wale tones is important to have a method for computing 

p °C om W-1 Specific thermal resistance for the variations of the pressure for a certain system 
insulation. under the most severe load variations to which 
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it can be subjected in service. These commonly 
consist of (1) switching off the load and (2) short 
circuiting. In the first case the oil in the cable is 
cooled and consequently the pressure decreases 
below the normal value until an amount of oil 
equaling the contraction of the oil has flowed 
into the cable. In consequence of the diminished 
pressure the electrical strength of the cable may 
decrease so much that breakdown occurs. On the 
other hand, in the case of a short circuit the oil 
expands rapidly because of the rise in tempera- 
ture, and the pressure in the cable will for some 
time increase and later decrease as the surplus 
of the oil flows towards the oil reservoirs of the 
system. In the case of a short circuit the pressure 
in the cable may increase so much that there is 
a risk of the sheath and the pressure armor being 
destroyed. 

All calculations of the pressure must be based 
upon the variations of the temperature in the 
cable with the given load variations. Several 
methods of computing the variations of the tem- 
perature and the oil expansion at load variations 
in cables placed under different circumstances 
have been described in the literature. Some of 
these methods will be mentioned later. 

When the oil expansion has been calculated, 
the variations of the pressure and the oil flow 
will be investigated. The dependence of the tem- 
perature on load variations has been treated 
thoroughly in the literature. That is not the case 
with the pressure and the oil flow. 

In many cases the increase or decrease in 
pressure for the greatest expansion of oil per unit 
time at load variations has been calculated under 
the assumption of stationary oil flow in the cable. 
For oil cable systems with different types of oil 
reservoirs, Emanueli' has developed a theory for 
the calculation of pressure and oil flow, in which 

.the pressure drop along the cable is calculated 
as if the oil flow were stationary. Shanklin and 
Buller? make the same assumption in their cal- 
culations of the pressure drop. 

H. Spanne, chief engineer at Sieverts Kabel- 
verk, Sundbyberg, Sweden, however, demon- 
strated about ten years ago in a factory bulletin, 
that the elasticity of the sheath and of the 


'L. Emanueli, Trans. A.1.E.E. 47, 186 (1928). 
2G. B. Shanklin and F. H. Buller, Trans. A.1.E.E. 50, 
1411 (1931). 
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pressure armor, if any, and the compressibility 
of the oil in an oil-filled cable greatly limit the 
pressure variations occurring at load changes. 
Spanne has also made an estimate of the pressure 
variations. 

In the _ present 
developed for 


article a method has been 
calculating the variations of 
pressure and the oil flow in an oil-filled cable with 
reservoirs at both ends, for an oil expansion given 
as a function of the time, with due regard to the 
flow resistance in the cable as well as to the 
elasticity of the sheath and the pressure armor, 
if any, and to the compressibility of the oil. Only 
those cases in which the expansion of oil per unit 
time and per unit length is the same at every 
point of the cable have been treated here. 


SIMPLIFYING ASSUMPTIONS 


The following assumptions are the most im- 
portant of those made in this paper: 


All deformations in the cable are purely elastic. 


1. 
2. No force is needed to accelerate the oil in the canals. 
3. 


The viscosity of the oil in the oil canals is constant. 
4. The pressure in the oil reservoirs is constant. 


Assumptions 1 and 2 are approximately satis- 
fied in most practical oil-filled cables. As the 
viscosity of the oil varies much with the tem- 
perature, 3, however, is not well satisfied. The 
difficulties involved in solving the problem are, 
however, greatly increased if assumption 3 is not 
made. Assumption 4 is well satisfied for cable 
systems with oil reservoirs consisting of oil tanks 
with atmospheric pressure at the surface of the 
oil. However, assumption 4 is often rather well 
satisfied in other types of oil reservoirs, too, as 
the oil expansion in the cable for a load variation 
often is small compared with the oil expansion 
corresponding to the total temperature vari- 
ation for which the system has been constructed. 
(The total temperature includes, for instance, the 
seasonal change of temperature.) This statement 
applies especially in the case of a short circuit in 
which the temperature of the conductors may 
certainly be rather high, but in which the total 
developed heat energy, and consequently also 
the oil expansion, often are small. 

An article is being prepared, in which assump- 
tion 4 is not made, due regard on the contrary 
being given to the characteristics of different oil 
reservoirs. 
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Fic. 1. Single core oil-filled cable, 77/v3kV, 1X95 mm. 


MATHEMATICAL FORMULATION 


The mathematical expressions for the flow in 
an oil-filled cable depend on how the flowing 
medium is defined. This, for example, can be 
defined as the mass of the oil, or as the volume of 
the oil. Here we shall base the calculations on a 
third and more convenient definition, defining 
the flowing medium as the volume of the oil 
reduced to constant pressure. 

The amount of medium per unit length at a 
certain point of the cable will then increase pro- 
portionally to the pressure because of the elas- 
ticity of the sheath and of the pressure armor, if 
any, as well as of the compressibility of the oil. 
We define the capacitance per unit length, c, as 
the ratio 6f the increase in amount of medium 
per unit length to the increase of pressure. The 
capacitance c is then composed of two parts, one 
part ¢,, arising from the elasticity of the sheath 
and the pressure armor, and another part ¢2, 
arising from the compressibility of the oil. 

When no pressure armor is present, the 
capacitance c,; for a cable with a sheath of thick- 
ness 6, diameter D, and modulus of elasticity E is 


arD* 
AGE 


C1 


(1) 


It is a trivial matter to develop the corre- 
sponding formula for cables with a pressure 
armor. 
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C2 is calculated from the compressibility of the 
oil, k, and the amount of oil per unit length, A. 


Co=kA. (2) 
c is thus 
C=C1+ Co. (3) 


The current of the medium is defined as the 
amount of oil which passes a section of the cable 
per unit time. The resistance per unit length is 
defined as the ratio of the pressure drop per unit 
length to the current. 

The flow at the moderate currents in question 
will be laminar, and r can consequently be 
expressed as follows,’ 


yv 
r = B—, (4) 
14 


where y is the density, v the kinematic viscosity 
of the oil, 7 a length which increases proportional 
to the linear dimensions of the section of the oil 
canal, and B a constant which merely depends 
on the shape of the oil canal and the way in 
which / has been defined. If / stands for the inner 
diameter of oil canals with a circular section one 
gets 


128 
B= —-10-7=4.15-10-7. 
9.81-7 





For a cable with a length L the total capaci- 
tance is C=Le and the total resistance R=Lr. 
The time constant RC for the oil flow in the cable 
is a measure of the rate at which disturbances in 
the state of equilibrium of the cable are smoothed. 
In order to give an idea of the order of magnitude 
of the time constant it has been computed below 
in a practical case for two different tempera- 
tures. 


Example 1 


The time constant at the temperatures 20°C 
and 40°C shall be calculated for an oil-filled 
cable 77vV3kV, 1X95 mm?, without pressure 
armor and with the principal dimensions as 
shown in Fig. 1. 


D=3.3 cm, 
d=0.6 cm, 


3 Handbuch der Physik (1927), Vol. 7, p. 103. 


wn 
+ 
~I 








6=0.2 cm, 
A=3.3 cm’, 
L=10° cm, 
E=1.7-10° kg cm~?, 
k=8-10- kg cm’, 
44.8 centipoise for 20°C 
y= 
16.4 centipoise for 40°C, 


r:3.3' 





C1 


Riis —=8.3-10-* kg~! cm’, 
4-0.2-1.7-108 


c2=8-10-5-3.3=2.6-10-* kg~! cm‘, 
c=1.09-10-* kg~' cm‘, 
128 vv 
r=——- 10°? -—__ =3.2-10°*-yp 
).81-2 0.64 
1.43-10-4 kg cm~* sec. for 20°C 
7 5.3-10-5 kg cm~* sec. for 40°C, 
((1.43-10~*) -(1.09-10-*) -10" 
= 1560 sec. for 20°C, 


RC=rcL* =< 
(5.3-107*) -(1.09- 107%) - 10" 





= 580 sec. for 40°C. 


This example shows that the time constant for 
the oil flow is not negligible compared to the 
intervals of time commonly necessary for 
smoothing the temperature in a cable under load 
variations. This means that the elasticity of the 
sheath and the compressibility of the oil will 
limit the pressure variations in the cable with 
changes in load. 

The differential equations and the boundary 
conditions for the oil flow will now be given. 
Because we have assumed that the cable has 

constant parameters, without making the treat- 
"ment of the problem less general we can suppose 
that the static pressure in the cable and conse- 
quently also the pressure at the endpoints of the 
cable are zero at any moment. Over the dynamic 
pressure distribution, calculated according to 
this assumption, we can thus superimpose the 
static pressure distribution in the cable. In this 
paper only the dynamic pressure distribution is 
considered. For a given oil expansion per unit 
length of the cable per unit time, S(t), we shall 
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Fic. 2. Coordinate system along a cable section terminated 
with oil reservoirs. 


calculate as functions of time and position the 
pressure V, the current J, and the amount of oil 
which flows through a cross section of the cable, 
M. 

A coordinate axis is placed along the cable 
with the origin at the midpoint (see Fig. 2). V and 
I then satisfy the following partial differential 
equations. 


OV 10°V 
c—=- -+S, (5) 
Oot or Ox? 
10V 
I= --—, (6) 
r Ox 
—L/2 
V=0 for x= 0<t<a, (7) 
V=0 for t<0, —L/2<x<L/2. (8) 


Equations formally identical with the equa- 
tions given above are obtained for several 
problems within other domains of physics. 
Suppose we have a metallic wire, isolated against 
heat, the length being L, with the heat conduc- 
tion resistance per unit length r, the heat capaci- 
tance per unit length c, and with the endpoints 
at the temperature zero. If an electrical current 
producing the effect per unit length S(t) is 
flowing in the wire, the temperature -V and the 
heat current J satisfy Eqs. (5)—(8). 

There also exists an electrical analogy to the 
problem in .question. If current sources of the 
strength per unit length S(‘) are supposed 
homogeneously distributed along a_ telegraph 
cable of length L, without inductance and leakage 
but with resistance per unit length 7, capacitance 
per unit length c (Thomson cable), and with the 
ends short-circuited, the potential V and the 
current I will satisfy Eqs. (5)—(8). 

In the oil-flow problem the mass of the oil 
corresponds to the inductance of the cable in 
the electrical analogy. The mass of the oil has 
no corresponding meaning in the problem of heat 
conduction. 


JOURNAL OF APPLIED PHYSICS 





by 
col 


tre 
re 
th 
mi 
els 


SO 
ay 


ec 
m 


re 
pl 


ed 


ie 
il 
e, 
le 


id 


al 


S 





It is convenient to solve first the problem given 
by Eqs. (5)—(8) for S(t)=5So-1, where So is a 
constant. 


PRESSURE AND FLOW WITH CONSTANT OIL 
EXPANSION S)-1 


The thermal analogy to this problem has been 
treated by Straneo,! Schaufelberger,’ and Wein- 
reich® for use in measuring the electrical and 
thermal constants of metallic wires. The authors 
mentioned above solve the problem by use of 
eigenfunctions following classical methods. 

The problem, however, is very conveniently 
solved by Heaviside’s operational calculus. The 
appropriate method has been described e.g., by 
Dahr’? in his treatment of the telegraphist’s 
equation for transmission lines in which electro- 
motive forces are induced. (See also Jeffrey’s 
treatment of related problems.*) For convenient 
reference the application of the method on the 
present problem is briefly outlined below. 

Introducing the operators 


LV (x, t) =v(x, p), (9) 


£I(x, t) =i(x, p), (10) 





Eqs. (5)—(8) are transformed into the following 
equations: 





1 d*v 
cpy=- So, (11) 
r dx? 
1 dv 
(=———, (12) 
r dx 
—L/2 
v=0 for x= (13) 
L/@. 


First the more simple problem of a point 
source 1 acting at the point x=¢€ is solved. For 
the pressure and the current in this problem the 
differential equations that are obtained from 
(11) and (12) by setting So=0 are valid, while 
to the boundary condition (13) are added the 
conditions 


v(¢ —0) =v(E+0), (14) 
a(£—0)+1=2(E+0). (15) 


As the solution of this problem one gets the 
following expression. for the pressure, which in 
this special case will be called w(x, &), 


L 





(pRC)! sinh(pRC)! 


1 : &s 
sinh(—-+~) (prc) sinh -—+— )(pRO) 
2 & 2L 


=w,(x, &) for ee 
(16) 


L 





(pRC)' sinh(pRC)} 


1 i 
| sinh( --+=)ipro sinh(--+—) (RC) 
| 2 L 2L 
| 
) 
| 
| 
| 


=we(x, ~) for <es 


The expression for v, in the case of a constant, homogeneusly distributed source, is now obtained 
by integration of Sow over the entire length of the cable 


L/2 


2 zr L/2 
v(x, p) =Sof w(x, pde=S4 f We(x, E)dE+ f w(x, pag] 
—L/2 > J, 


:. 2 ; £ 
sinh (—+—) (eRC)'—sinh( —-+— )(pRC)! 
1 a £ s & 


—L/2 





= Sul) — 


5 W. Schaufelberger, Ann. d. Physik 7, 589 (1902). 
6 W. Weinreich, Zs. Math. (Leipzig) 63, 1 (1914). 


(17) 
PRC sinh(pRC)! 


4P. Straneo, Roma. Rend. Acc. dei Lincei (Ser. 5) 7, 206 (1898). 


7K. Dahr, A Course of Integrational and Operational Calculus (V. Pettersons Bokindustriaktiebolag, Stockholm, 


1935), p. 44. 


8H. Jeffreys, “Operational Methods in Mathematical Physics,’’ Cambridge Tracts in Math. and Math. Ph. 23, 


49 (1931). 
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Fic. 3. Generalized functions used in the calculation of the pressure at the midpoint of a cable section with a time constant 
for the oil flow RC. 


Vi’(t/RC) applies in the case of an oil expansion S(t) =1-1. 
V2'(t/RC) applies in the case of an oil expansion S(t) =(d/dt)(1 +1). 
Vs/(t/RC) applies in the case of an oil expansion S(t) =e~™ RC), 
Vi(t/RC) applies in the case of an oil expansion S(t) =1 (t)3 +1. 


The time function for operators of this type 
can be obtained either by expanding the operator 
in an infinite series of terms containing exponen- 
tial functions and transforming this series, term 
by term, or by applying Heaviside’s expansion 
theorem. The first method corresponds to 
expressing the solution in reflected waves as 
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demonstrated e.g., by Wagner,’ whereas by the 
last method the solution is expressed in eigen- 
functions. Here we shall use the last method. 
According to Heaviside’s expansion theorem 
the time function H(t) that corresponds to an 


*K. W. Wagner, Operatorenrechnung (Johann Am- 
brosius Barth Verlag, Leipzig, 1940), p. 196. 
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operator f(p)/g(p), where f(p) and g(p) are 
functions of p, which satisfy certain conditions, 
and where g(p)=0 has the simple roots pn, 
n=1,2, 3, -+-,is 











f(p) f(0) f(Pnjern 
H(t) = £-'—- =—_+ . (18) 
g(p) (0) Png’ (Pn) 
If we write the operator v (17) as follows 
f(P) 
v= . (19) 
g(P) 


where 


1 
f(b) =SorL?— — | sinh(pRC) 
(pRC)! 


1 x 
— sinh(- +") (pRC)} 
2L 


:. & 
+sinh( —- + )ivro'| (20) 
2 LL 


g(p) =(pRC)' sinh(pRC)}, (21) 


v=f(p)/g(p) satisfies the conditions for the use 
of the formula (18). The roots of the equation 


g(p) =(pRC)} sinh(pRC)'!=0, (22) 
different from zero are 
(p»aRC)'=+jnx, n=1,2,3,---, (23) 


or 
Pa=—-—, n=1,2,3,-°-. (24) 


The time function V = £~'v is then 


1 isx rx? 
V(x, t) = Sui4|—— ( ) 
8 2\L 


4 «2% (—1)" x 
—-—)> cos—(2n+1)x 
m® n=0 (2n+1)? L 





RC 


(2n+1)*3t 
I}: (25) 


xexo| _ 


For 4=0 the term containing the sign of sum- 


mation expresses the Fourier series of }—}(x«/L)* 


in the interval —L/2<x<L/2. V can thus be 
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expressed in the following way 


4 » (—1)" 
V(x, t)=SorL?— > — 
m® n=0 (2n+1)? 


(2n+1)?x?t 
x | 1 exp{ -——"II. (26) 
RC 





x 
cos—(2n+1)x 
L 


For x =0 this leads to 


4 - (-— 1)” 
V(0, t) = Sor L?— >> ———— 
Tr? n=0 (2n+1)* 


(2n+1)?x°t 
x | 1 exo | ; (27) 
' RC 


[V(O, t)/SorL?|=Vi' has been plotted as a 
function of ¢/RC in Fig. 3. 
I is calculated from (6) and (26) 


4 » (—1)" x 


I(x, )=SoL— > sin—(2n+1)x 
wr’ n—0(2n+1)? L 


| (2n+1)?x*t 
Es exo] | . (28) 
| RC 


For x=+L/2 this leads to 


L 4 @ 1 
(+=, t) = Sob — 
2 1? n=—0 (2n+1)? 


(2n+1)?3*t 
x |1 exp] - |} (29) 
RC 











[1(+(L/2), t)/ Sol j=)’ has been plotted as 
a function of ¢/RC in Fig. 4. 

The amount of oil which flows through a cross 
section of the cable, M, is calculated by inte- 
grating the expression for J, given above. 


‘ t 
M(x, t) -{ I(x, \)dX= SorcL* | — 
; | LRC 


4 ee (—})* 


x 
-> - —sin—(2n+1)r 
rin () (2n+1)4 L 


{ (2n+1)*3*t 
x|1 -exp| - —- || . (30) 
RC 
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Fic. 4. Generalized functions used in the calculation of the oil current at the endpoints of a cable section with a time 
constant for the oil flow RC. 


1;'(t/RC) applies in the case of an oil expansion S(t) =1-1. 
I2/(t/RC) applies in the case of an oil expansion S(t) =(d/dt)(1-1). 
1,/(t/RC) applies in the case of an oil expansion S(t) =e—@!/R€ -1, 
I(t/RC) applies in the case of an oil expansion S(t) = 1/(t))- “1. 
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Fic. 5. Generalized functions used in the calculation of the amount of oil which passes the endpoints of a cable section 
with a time constant for the oil flow RC. 


M;'(t/RC) applies in the case of an oil expansion S(t) =1-1. 
M?2'(t/RC) applies in the case of an oil expansion S(t) =(d/dt)(1-1). 
M;/(t/RC) applies in the case of an oil expansion S(t) =e~@/R€ -1, 
M4’(t/RC) applies in the case of an oil expansion S(t) =1/(t)! +1. 








For x=+L/2 this leads to In consequence of the linearity of the partial 

L t i‘ 1 differential equations for the flow, the pressure, 

u(a=, ') = Sul] —-— > — current and amount of oil which flows through 
2 2RC x‘ n=0 (2n+1)! 


a cross section of the cable for an arbitrary oil 

(2n+1)?xt l expansion S(t)-1 can be calculated from the cor- 

x1 —exp| - "||. (31) responding solutions, given above, valid for an 

RC oil expansion of the form So-1l. The linearity 

[M(+(L/2), t)/+SorcL*]= My’ has been plotted of the differential equations also causes the 
as a function of ¢/RC in Fig. 5. action (the pressure, the current, or the amount 
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of oil) of an oil expansion given as a sum of time 
functions, to be obtained as the sum of the 
actions of the individual terms of which the oil 
expansion consists. 


First it is necessary to mention some methods 
described in the literature for calculating the oil 
expansion in oil-filled cables for different types 
of load variations. Then the pressure, current, 
and amount of oil which flows through a cross 
section of the cable will be calculated for some 
important time functions of the oil expansion. 


OIL EXPANSION FOR DIFFERENT TYPES OF 
LOAD VARIATIONS 


During the last fifty years a great many papers 
dealing with temperature variations in cables 
subjected to load variations have been published. 
For the purpose of calculating the oil expansion 
in cables subjected to ordinary load varia- 
tions the methods suggested by Shanklin and 
Buller’ and Miller and Wollaston" seem to 
be very useful. Using either of these methods 
the oil expansion is obtained as a sum of ex- 
ponential Soe"! +1. 
previously been mentioned, pressure, current, 
and amount of oil which flows through a cross 
section of the cable are calculated by summing 
the quantities corresponding to the individual 
exponential terms. In a later part of this paper, 
therefore, the pressure, current, and amount of 
oil for an oil expansion Soe~'/7-1 will be cal- 
culated. 

The above mentioned methods for calculating 
the variations of temperature and oil expansion 
in oil-filled cables at load variations are rather 


terms From what has 


inconvenient when the load variation in question 
is a short circuit. The temperature function in 
the case of short circuit is better computed on 
the basis of other simplifying assumptions, as the 
interesting part of the temperature variation in 
the case of a short circuit occurs in a relatively 
thin layer of the insulation near the conductor, 
in which the short-circuit current is flowing. 
Herlitz” has worked out a theory for temperature 
during and after a short circuit in conductors in 





10G. B. Shanklin and F. H. Buller, Gen. Elec. Rev. 34, 
523 (1931). 

“kK. W. Miller and F. O. Wollaston, Trans. A.I.E.E. 
52, 98 (1933). 

7. Herlitz, Teknisk Tidskrift, Elektroteknik 71, 204 
(1941). 
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contact with paper-oil insulation, whereby con- 
sideration has been given to the heat capacitance 
of the conductors and to the temperature coef- 
ficient of the conductor resistance. The oil- 
expansion computed on the basis of this theory 
is expressed by a function involving exponential 
functions and error integrals. 

For use in the present paper we shall develop 
a very approximate and very simple expression 
for the oil expansion in an oil-filled cable in the 
case of short circuit in conductors in contact 
with paper-oil insulation. We hereby make the 
simplifying assumptions (1) that the conductors, 
in which the short-circuit current flows, have a 
plane surface, (2) that the insulation can be 
regarded as semi-infinite, (3) that the short 
circuit is instantaneous, and (4) that the tem- 
perature in ‘the conductors remains constant 
after the short circuit has occurred. The surface 
of the conductors per unit length of the cable is 
called f, the specific heat-conduction resistance 
of the insulation, p, the specific heat of the 
insulation with respect to a unit volume, q, and 
the temperature, which the conductors attain at 
the short circuit, Jo. If x denotes the distance 
from the surface of a conductor to a point of 
the insulation, the temperature field J in the 
insulation is expressed by 


x(pq)? as 
I(x, =o 1-0 —) |. (32) 
2(t)! 


The heat effect per unit length of the cable, P, 
transferred from the conductors to the insulation 
is then 


f 00 1 q\' 1 
P= —"( ) * foe( ) A, 
p\dt7 0 (mr)! p/ (t)} 


If we denote the effective volumetric expansion 
coefficient of the insulation by 6’, the oil ex- 
pansion per unit time per unit length of the 


(33) 


cable, S, is 


B’P(t) Ko 
S(t) =——_=—-: 1 (34) 
q (t)? 
where 
1 B' fio 
a= : ames (35) 
(wr)? (pq)? 


18H. S. Carslaw, Introduction to the Mathematical Theory 
of the Conduction of Heat in Solids (Dover Publications, 
New York, 1945), second ed., p. 35. 
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Pressure, current, and amount of oil which 
flows through a cross section of the cable for an 
oil expansion according to the very approximate 
method of computation set forth here will be 
calculated in a following paragraph. The pressure 
so obtained will be larger than the actual pres- 
sure, especially in cables with small conductor 
areas, since the temperature of the conductor 
actually decreases after the short circuit instead 
of being constant as here supposed. 

In a future paper, the pressure, current and 
amount of oil for the case of a short-circuit 
current flowing in conductors in contact with 
paper-oil insulation will be worked out for an 
oil expansion calculated by a more exact method 
of computation resembling that given by Herlitz. 

The expansion of oil in canals, if any, in the 
conductors is difficult to calculate, since the 
heat transfer from the conductors to the oil in 
the canals takes place both by conduction and 
by convection. One is safe in assuming that the 
oil in a canal of any conductor has at any 
moment the same temperature as the con- 
ductor. In the case of an instantaneous short 





circuit the oil in a conductor canal will then be 
subjected to a sudden temperature rise resulting 
in a sudden expansion of the oil. If, as above, #o 
denotes the temperature rise of the conductors 
because of an instantaneous short circuit, and if 
a denotes the volume of the conductor canal per 
unit length, and 6 the volumetric expansion coef- 
ficient of the oil, the oil expansion per unit length 
per unit time, S, for the oil in the conductor canal 
in question is, using the terminology of the 
operational calculus, expressed as follows 


d 
S(t)=—(No-1), (36) 
dt 


where 


No=av8. (37) 


Pressure, current, and amount of oil will be 
computed for an expansion of oil of this kind in 
the following paragraph. The quantities so ob- 
tained turn out to be more convenient for cal- 
culating pressure, current, and amount of oil for 
an arbitrary expansion of oil than the expressions 
obtained above for an oil expansion So-1. 


PRESSURE AND FLOW FOR AN INSTANTANEOUS OIL EXPANSION (d/dt)(No-1) 


Pressure, current, and amount of oil are calculated by the expressions earlier obtained for an oil 
expansion So-1 by differentiating these with respect to time, and writing No instead of So. In this 


way one obtains 








No4 «= (-—1)" x (2n+1)?x*t 
Vix, b=—- dy — cos—(2n+1)z exo| - -- -| (38) 
c mn=02n+1 L RC 
For x =0 this leads to 
No4 x (—1)" (2n+1)?xt 
V(0, .)=—- > — exp| — ————| (39) 
c mn=0 2n+1 RC 
[cV(0, t)/ No ]= V2" has been plotted as a function of t/RC in Fig. 3. 
No x x (2n + 1)°3rt 
I(x, t)=—4 © (—1)" sin—(2n+1)2 exp| ~ “oT. (40) 
rcL n=0 , RC 
For x=+L/2 this leads to 
L No (2n+1)?x*t 
i( + 1)=+ <4 ¥ ex - ———| (41) 
2 rcL n-0 RC 
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[rcLI(+(L/2), t)/+No]=J,' has been plotted as a function of ¢/RC in Fig. 4 








x 4. (-1)" x (2n+1)*x*t 
M(x, t)=NoL{——-— > sin- (2n+1)x exp| — ——|| 
Lx* n=—0 (2n+1)? L RC 
4a (-1)" x | (2n+1)?3*t | 
=N,L— >> —— sin—(2n+1)z-: 1 -exo| - — ; (42) 
wr? n—0(2n+1)? L | RC 


For x=+L/2 this leads to 
L Ge 1 (2n+1)2x% 
u( +> ') =+N,L— >}, ———i1 —exp| || (43) 
2 x? no (2n+1)? rc J 
[M(+(L/2), t)/+NoL]=M,’' has been plotted as a function of t/RC in Fig. 5. 





PRESSURE AND FLOW FOR AN ARBITRARY OIL EXPANSION S(t)-1 


If any action (pressure, current, or amount of oil) of an oil expansion (d/dt)(1-1) is A(t), then 
because of the linearity of the partial differential equations, valid for the oil flow, the action G(t) 
of an oil expansion S(t) is expressed by Duhamel’s integral 


G(t)= f S(A)A(t—D)ad. (44) 


Using this equation in connection with the above obtained expressions for pressure, current, and 
amount of oil for an oil expansion (d/dt)(N-1), (38), (40), and (42), one obtains 


14 « (—1)" x (2n+1)?x*t] (2n+1)*9*d 
V(x, )=-- > —— = (2n+1)x exp| — ~~ \f S(d) exp| ———_—__— Ja. (45) 
0 








c wn-02n+1 RC RC 
1 (2n+1)?x°t (2n+1)?x?Xd 
I(x, t)=—4 > (—1)" dot “(2n-+1)x exp| —— ~“\f S(A) exp] ~ aon (46) 
rcL n=0 L RC 0 RC 
4 #® (—1)" 


M(x, )=L— > — 
aw? n=0 (2n+1) 





x 
sin—(2n+1)x 
2 L 


(2n+1) rt (2n-+1)?3*d 
x{ f S(A)dxA— exp| — ——lf 50s) exp| —— al. (47) 


PRESSURE AND FLOW FOR AN EXPONENTIALLY DECREASING OIL EXPANSION Sye~/?-1 





Pressure, current, and amount of oil for the oil expansion Sge~'/?-1 is calculated by substituting 
S(t) = Soe"? in Eqs. (45), (46), and (47). 





4 @ (—1)*" 
V(x, t)=SorL?- > - deers “(2n+1)m 
ax n=0 (2n+1)((2n+1)?x?—(RC/T)) L 


x 








RC t (2n+1)?x*t 
exp( -— --) —exp| _--—— ||: (48) 
T RC RC 
For x =0 this leads to 


4 w (—1)" RC t (2n+1)?*3*t 
V(0, t)=SorL?- > - fe »(-=— —exp| -——_—|I. (49) 
a n=0 (2n+1)((2n+1)*x?—(RC/T)) T RC C 
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[ V(0, t)/SorL? ]= V3’ has been plotted in Fig. 3 as a function of t/RC for a=RC/T =0.02, 0.05, 


0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50, and 100. 


va (—1 " 
I(x, t))=SoL4 ¥ 
<0 (2n+1)’x?—(RC/T)L 





For x=+L/2 this leads to 





L * 1 
i *-—, ') =+5,)L4 > 
2 n=0 (2n+1)*2?—(RC/T) 


x 
sin—(2n+1)z 





RC t of (2n-+1)2x% - 
os Fe —— oe ll: (50) 


RC ¢t (2n+1)?xt 
jexr( a --) —exp| - "||. (51) 
I RC RC 


[J(+(L/2), t)/+SoL ]=I;' has been plotted in Fig. 4 as a function of t/RC for the same values 


of a=RC/T as given above. 


(5) 





oF RC 1 4 «@ 
M(x, t) = Swrek|— -41 —exp( —- -.)|-= =. —— 
LR T RC mr? n=0 (2n+1)?((2n+1)?x?—(RC/T)) 


_ RC t (2n+1)?xt 
Xin Qn+-1)x] exp( -————) ex ——"]]. (52) 


For x=+L/2 this leads to 


| 3 
u( > ) =+SorcL’ 





T RC RC 


4 » 1 





RC t 
Fm PSE, 
2RC T RC m* n=0 (2n+1)?((2m+1)?*x?—(RC/T)) 





RC 1 (2n+1)?x*t)| 
x| exe( -—= — exp -—" \|\- (53) 


[ M(+(L/2), t)/+SorcL* |= M;’ has been plotted 
of a=RC/T as given earlier. 
A practical example will demonstrate the use of 


Example 2 


For a certain cable the oil expansion per unit 
length per unit time for a sudden switching on 
the load is 


S(t) = Si [e~ 4000) _ () Ge— (¢/400) _ Q) De (t/ 100) 7 


where S; is a constant. It is desired to calculate 
the variation of pressure at the midpoint of the 
cable, partly for switching on the load at the time 
t=0, after which the load remains switched on 
an infinite time, partly for switching on the load 
at the time t=0 and switching off the load at the 
time t= 200 sec. The computation of the pressure 
will be carried out for a time constant of the 
cable RC=0, 1000, and 4000 sec. 

For RC=0 the pressure V in switching on the 
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T RC RC 


in Fig. 5 as a function of t/RC for the same values 


the curves. 


load is 
(V(t)) rc-o= § Sir L?Le~ / O™ 
— 0.6e-‘¢/ 40) — 0), 2e- (4/100) 7 


[(V(t))rc-0/§SirL*] has been plotted in Fig. 6. 


For the case where the load is switched off 
after 200 sec., the pressure W is 


(V(4)) remo for t<200 me. 
(W(t) rc=0= 4 (V(t) rc-0— (V(t — 200)) re=o 

for ¢>200 sec. 

[(W(t)) rc-0/$SirL? | has been plotted in Fig. 7. 


For RC=1000 sec. one obtains for the ex- 
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ponential terms in the expression for S 


RC 1000 


apa 95, 
T; 4000 
RC 1000 
—=—=2.5, 
T, 400 
RC 1000 
— =——-=10. 
T; 100 


The pressure V in switching on the load is, 
accordingly, 


t 
(Vb) nenssonwe.= SyrL'} v'(—_)) 
1000 —T 
t 
- (a) 
10007 J». 
t 
o(r(L)) | 
1000 — 


The are taken from Fig. 3. 
L(V (2) rc=1000 see./$57L?] has been plotted in 
Fig. 6. 

The the load remaining 
switched on 200 sec. is calculated in the same 
way as in the case RC=0 treated above. 

[(W(t)) rc=1000 see./#5irL? | has been plotted in 
Fig. 7. 


functions V;’ 


pressure W for 


06 





a7 


o 
10 100 1000 


10000 
J 


Fic. 6. Relative pressure at the midpoint of a cable in 
switching on a load giving an oil expansion S(t) = S;[e~*/49 
—0.6e~*/4°° —0.2e~*/'°] for different time constants of the 
oil flow. 
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For RC = 4000 sec. one obtains in a similar way 


t 
( V(t))re $000 see. = Sir ( vi(——)) 
4000 eat 
t 
-06(vs(.)) 
4000 — 
l 
-0.2(v(—_)) } 
4000 a=40 


L(V (t)) re=s000 see. $5irL*? ] has been plotted in 
Fig. 6 and [(W(t)) rc=s000 sec./ §S17L* | in Fig. 7. 

Figure 6 shows that in the present case the 
variation in the computed pressure, resulting 
from taking into account the elasticity and the 
compressibility of the cable, is but little even 
with a time constant for the oil flow so large as 
4000 sec. If all the terms in the expression for S 
had the same sign, the variation in the calculated 
maximum pressure, however, would have been 
considerable, with the time constants for the oil 
flow occurring in the present Mostly, 
however, S is composed of exponential terms with 
different signs in such a way that S at the time 
t=0 assumes a small value, thereafter increases 


case. 


to a maximum, and then decreases to zero in 
resemblance to the case occurring in the exam- 
ple. With the values of the time constant for the 
oil flow, and the time constant in the first ex- 
ponential term in the expression for the oil 
expansion practical the 
maximum pressure in switching on a constant 
load calculated for RC=0 seldom diverges any 


occurring in cases, 











10 <a 100 (000 70000 


Fic. 7. Relative pressure at the midpoint of a cable in 
switching on for 200 sec. the same load as occurring in 
Fig. 6 for different time constants of the oil flow. 


JOURNAL OF APPLIED PHYSICS 











0 = = oe 2 





J 4 5 


Fic. 8 F(x) =e *Sf,7etdé together with the approximations F(x)=x and F(x)=1/2x valid for small, respectively 
for large values of x. 


considerable amount from the pressure calcu- 
lated for the actual value of RC. 

This is not the case in the calculation of the 
maximum increase of pressure and the following 
maximum decrease of pressure occurring in 
switching on a constant load for a time interval 
that is small relative to the time constant for 
the oil flow in the cable. Figure 7 shows thus 
that the maximum pressure for the load in 
question with RC = 1000 sec. comes to 76 percent, 
and with RC=4000 sec. only 34 percent of the 
maximum pressure with RC=0. 

The following decrease of pressure is reduced 
too by the elasticity and the compressibility of 
the cable. 


PRESSURE AND FLOW FOR THE OIL 
EXPANSION K,/(t)}-1 


Pressure, current, and amount of oil for the 
oil expansion Ko/(t)'-1 is calculated by sub- 
stituting S(t)=Ko/(t)' in Eqs. (45), (46), and 
(47). 

The results in the present paragraph are most 
simply expressed by the, function F(x), defined 
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in the following way: 


F(x) =e- of efdé. (54) 


This function has been tabulated by Miller 
and Gordon." F(x) has been plotted in Fig. 8 
for 0<x<5 together with the functions y=2x 
and y=1/2x, which approximate the function 
F(x) for small respectively large values of the 
argument x. The advantage in using the func- 
tion F(x) instead of the function 


ne f edi, 
0 


as used e.g., by Dawson,'® is that F(x) is limited 
and more convenient for interpolation than f(x)." 
Pressure, current, and amount of oil can now 


4 W.L. Miller and A. R. Gordon, J. Phys. Chem. 35, 
2785 (1931). 

15H. G. Dawson, Proc. Lond. Math. Soc. 29, 519 
(1897-98). 
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be expressed in the following way: 


» 2 « (—1)" 
(x, V(%,)= Ki“ ) L- > a 
c mr? n=0 (2n+1)? 


x t 3 
cos An+1)xF( (2n-+1)x(——) ). (55) 
L RC 


For x =0 this leads to 
r\' 8 «@ 1 
V(0, )=K,(") L—> —— 
c 3 n=0 (2n+1)” 


tr} 
#(antiye(—)'). 
RC 


[V(0, t)/Ko(r/c)'L = V4’ 
a function of t/RC in Fig. 3 


(56) 
has been plotted as 


Ko 8 « (—1)" 
I(x, ))=—— 
(rc)) x n- 2 n-+1 


j 

x sin- - ~(2n-+1)rP( 2n+1r(— ) ). (57) 
L RC 

hg/em* 


2 
| 
| 
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Fic. 9. Pressure V at the midpoint of a cable, occurring 
at a short circuit in the cable as specified in the text, for 
20°C temperature of the oil in the oil canal. V= V2+ V4, 
where V2 corresponds to the expansion of the oil in the 
canal, V4 to the expansion of the oil in the insulation. 


For x=+L/2 this leads to 


L Ky 8 oo 1 
fad sjoate ts 
2 (rc)! w n—0 2n+1 


tr} 
xP(2n+1)x(—) ). (58) 
RC 
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(1(+(L/2), t)/+(Ko/(re)!)]=I' has 
plotted as a function of ¢/RC in Fig. 4. 


xfty\! 8 = (—1)" 
M(x, t) =Ke(re)'1-|2~(—_) —-—>) - 
LX\RC 4 n=0 (2n+1)? 


x t \3 
Xsin—(2n-+1)rF( (2n-+1)e(— ) )} (59) 
L RC 


been 
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Fic. 10. Pressure V at the midpoint of the same cable 
and for the same short circuit current in this cable as in 
Fig. 9 for 40°C temperature of the oil in the oil canal. 
V = V2+ V4, where V2 corresponds to the expansion of the 
oil in the canal, V4 to the expansion of the oil in the 


insulation. 
For x=+L/2 this leads to 
t ; 
u( 2— t )- +K (rc) 14 (—-) 
RC 
8 oe 1 


r( 2n+1)x(- -) )| (60) 
RC 


[M(+(L/2), t)/+Ko(rc)'L]=M,’ 
plotted as a function of t/RC in Fig. 5 
The use of the formulas worked out in this 
paragraph is explained by the following example. 


1 n=0 (2n-+1)8 


has been 


Example 3 


It is desired to calculate the maximal pressure 
rise in the single-core cable described in Ex. 1 for 
a cable length of 1 km in the case of a short 
circuit of 1 sec. duration, during which the con- 
ductor and the oil in the conductor canal are 
heated 10°C. It is supposed that the heat effect 
of the secondary current in the sheath is neg- 
ligible. The pressure will be calculated for the 
temperature of the oil in the canal 20°C and 
40°C. 

Besides the data for the cable given in Ex. 1 
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the following data are needed: 


p=550°C cm W-, 
8 =0.0007°C-*, 


g=2°C—' cm W sec., 
6’ =0.0004°C—", 


From the data given in Ex. 1 are further cal- 
culated 


f=n-1.3=4.1 cm, 
a= (2/4) -0.6?=0.28 cm’. 


The increase in pressure at the midpoint of the 
cable V.(0, ¢) caused by the expansion of the oil 
in the canal is according to (37) and (39) 


1 t 
V.(0, )=a08-V'( ) 
c RC 


= 0.28-10-0.0007 -— 


1 ( t ) 
———— |,’ — 
1.09-10-8 RC 


t 
= .svs(—) kg cm~? 
RC 


where V2'(t/RC) is taken from Fig. 3. 

The increase in pressure at the midpoint of the 
cable V4(0,¢) caused by the heat transferred 
from the conductor to the insulation is according 
to (35) and (56) 


B’fdoL(r)! t 
V,(0, )-————_vi( ) 
(xpqc)} RC 
0.0004 -4.1-10-10°V,/(t/RC) 
. (x-550-2-1.09-10-*)3 











= 
(5.3-10-4)! 


t 
10.1 vi(—) at 20°C 
RC 
t 
savi(—-) at 40°C 
RC 


where V4'(t/RC) is taken from Fig. 3. 

Vo, Vs, and V= V2+ V, for 20°C temperature 
of the oil in the oil canal have been plotted in 
Fig. 9, and V2, V4, and V= V2+ V, for 40°C oil 
temperature in Fig. 10. 
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From Fig. 9 and Fig. 10 is read the maximum 
pressure 


5.4 kg cm~ for t=95 sec. with 20°C 
3.8 kg cm~ for t=30 sec. with 40°C, 


Y on —_ 


In working out the theory for the oil expansion 
after a short circuit we made the simplifying 
assumption that the temperature of the con- 
ductors remained constant after the short circuit. 
It is interesting to make an estimate in the 
present case of how much the conductor tem- 
perature has decreased from the moment of 
short circuit to the moment of occurrence of 
maximum pressure. 

With the idealizing assumption mentioned 
above, namely, that the conductor temperature 
remains constant, the effect per unit length of 
the cable, P, transferred to the insulation was 
given in (33). The amount of heat per unit 
length, Q, transferred to the insulation from the 
moment of short circuit to the time ¢ is then 


t 2 q\} 
a= f POydr=—foo(*) (t)*. 


Hereby we calculate the amount of heat per 
unit length in the present case transferred to the 
insulation from the moment of short circuit to 
the moment ¢;, when the maximum pressure 


occurs 
2.3 
10:(—) 
550 


(95)!=26 W sec. cm™! for 20°C. 
ene 15 W sec. cm! for 40°C. 


2 
(t;) =——-4.1- 
— (x)! 


These transferred amounts of heat have to be 
compared with the amount of heat per unit 
length, Q;, accumulated in the conductor and in 
the oil canal. The specific heat for copper in 
respect to a unit volume is 3.4°C~' cm~ W sec. 
Neglecting the amount of heat accumulated in 
the oil canal we get 


Q,=0.95 -3.4-10=32 W sec. cm™. 


The amount of heat transferred from the con- 
ductor to the insulation from the moment of 
short circuit to the moment, when the maximal 
pressure, calculated according to the simplifying 
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assumption mentioned above,. occurs, is con- 
sequently not negligible compared with the 
total amount of heat accumulated in the con- 
ductor. The preparation of a more exact theory , 
as mentioned above, seems to be justified. 


The work here presented has been carried out 
at Sieverts Kabelverk, Sundbyberg, Sweden. 
1 thank Sieverts Kabelverk and especially Mr. 
H. Spanne, chief engineer at Sieverts Kabelverk, 
for permission to publish this paper. 
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In this paper the field of a point charge inside a hollow, 
infinitely long circular cylinder is studied. The case of an 
axial point charge is treated in detail. Three different 
methods are developed. The first method shows how to 
calculate the induced charges at the surface of the cylinder 
without explicit knowledge of the potential itself. The 
surface charge-density function is obtained as the solution 
of a Fourier-type integral equation. Then the potential 
caused by these charges is calculated. The second method 
works in the opposite direction. Here the potential is 
obtained as solution of a boundary value problem, followed 
by the calculation of the corresponding charges at the 
surface of the cylinder. 

The integral, obtained for the surface charge density, is 
transformed by contour integration. Although the resulting 
series is very useful for numerical purposes, a stronger 
method is necessary, in order to calculate the charge 
density just opposite the point source. Figure 1 shows the 
calculated values of the induced surface charge density. 
Several approximations for the charge-density function are 


1. INTRODUCTION 


REVIOUS works of Schelkunoff,! Weyrich,* 
Buchholz,’ and Oberhettinger,* among 
others, already contain complete information 


with regard to fields in circular pipes, as produced 


by Hertzian dipoles. Therefore, the generation 

'S. A. Schelkunoff, ‘‘Modified Sommerfeld’s integral 
and its applications,” Proc. I.R.E. 24, 1388 (1936); and 
“Transmission theory of plane electromagnetic waves,” 
ibid. 25, 1457 (1937); with bibliography. 

* R. Weyrich, ‘Ueber einige Randwertprobleme, insbe- 
sondere der Elektrodynamik,” J. reine angew. Math. 172, 
133 (1934). 

*H. Buchholz, ““Gekoppelte Strahlungsfelder im kreis- 
zylindrischen Hohlleiter,’’ Ann. d. Physik 39, 81 (1941). 

*F. Oberhettinger, ‘Ueber ein Randwertproblem der 
Wellengleichung in Zylinderkoordinaten,” Ann. d. Physik 
43, 136 (1943). 


562 


considered, in connection with recent work of Weber (cf. 
Table I). 

The same is done for the field inside the cylinder. 
Various formulas are given which allow of numerical 
calculations. In Table II some calculated values of the 
potential are shown. 

The third method is based on the theory of Fourier- 
Bessel-Dini series. The potential is developed in terms of 
discrete normal solutions of the potential equation in 
cylindrical coordinates. The coefficients in this develop- 
ment can be derived from the behavior of the potential in 
the immediate neighborhood of the primary source. 

Furthermore, it is emphasized that the study of the 
above potential problem can serve as a guide in questions 
of wave propagation in hollow circular cylinders. In this 
connection the third method is shown to be extremely 
useful, as it enables us to calculate directly the fields in the 
far zone from that in the immediate neighborhood of the 
exciting source. This new method is demonstrated in case 
of acoustic waves inside a cylinder, caused by a harmoni- 
cally vibrating point source. 


of special types of wave patterns, caused by 
given current and charge distributions on the 
exiting antennas, can be dealt with quantita- 
tively. The calculations in this domain of electro- 
magnetics, however, are far from elementary, 
such that a non-mathematically minded engineer 
might hardly recognize the usefulness of those 
theoretical treatments. Hence, every attempt to 
diminish the necessary amount of mathematics, 
may be worth while to consider. Until now a 
considerable amount of complex function theory 
is required. We believe, however, that this can 
be greatly avoided if one bases oneself on the 
theory of Fourier-Bessel-Dini series. This method 
is undoubtedly within reach of every engineer, 
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because the analogous conventional theory of 
the ordinary Fourier series is common knowledge 
nowadays. 

In this paper we will not consider the vector 
wave propagation as occurring in electromagnetic 
theory. Those topics will be treated elsewhere. 
The new method will be demonstrated in the 
case of scalar wave propagation only. 

The major part of the present article, however, 
will be devoted to the study of the electrostatic 
field of a point charge inside the cylinder. Conse- 
quently, only a minor part of the mathematics 
developed here is claimed to be useful in wave 
guide theory. 

Beginning with the simplest problem, we con- 
fine ourselves to the case of a point charge 
situated on the axis of the cylinder. Later on 
we will briefly indicate how to extend the theory 
to excentric point charges. 


2. STATEMENT OF THE SIMPLEST PROBLEM. 
CALCULATION OF THE CHARGES INDUCED 
ON THE SURFACE OF THE CYLINDER 


Let a be the radius of the cylinder, and let the 
point charge have a value +1. Furthermore, 
let z, p, ¢ denote the right-handed cylindrical 
coordinates of a field point P. Henceforth the 
origin may coincide with the given point charge, 
and the z-axis be taken along the axis of the 
cylinder. 

The electric force E can be derived from a 
potential V according to E=—gradV. It is 
evident that the z-axis is an axis of symmetry, 
so that V depends only on z and p. In addition, 
V is even in z. 

To calculate E and V, we first follow a rather 
unconventional method, by calculation of the 
charges, induced on the wall of the cylinder, 
without explicit knowledge of the potential itself. 
Afterwards the field is readily found from the 
original point charge and these induced charges. 
The surface density of the induced charges will 
be denoted by n(z). 

We may remark that V is the so-called Green 
function of the first kind. It is a regular solution 
of the potential equation, except for a singularity 
in the origin, so that the function V(P)—1/r 
remains finite for r-0; furthermore, V=0 on 
the boundary, and both V and gradV tend to 
zero, if P tends to infinity inside the cylinder. 
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As is well known, there is one, and only one, 
function V(z, p) fulfilling all these conditions 
simultaneously. Once this function is obtained, 
the charges on the boundary of the cylinder can 
be computed from 


1 0V 
n(z)=——, (p=a). (1) 
4r Op 


Let us now show how 7(z) can be found without 
explicit knowledge of the potential V. Be g(z, p) 
a given regular solution of the potential equation 
inside the cylinder, and suppose g and dg/dz to 
be bounded in that domain. By Green’s theorem, 
we at once express the value of g in the origin 
as a surface integral over its boundary values, 
namely 


1 OV 
g(0, 0) = —— [e—ae- - fats a)n(z)do. (2) 
An Op 


In this do means the surface element of the 
boundary of the cylinder; the integration in (2) 
extends over this boundary. 

Omitting the left-hand integral in (2), the 
remaining equality can be interpreted as an 
integral equation for n(z), which may be readily 
solved, if a convenient choice is made for the 
rather arbitrary function g. We take 


g(z, p) =e—*To(Ap) (3) 


where J is the Bessel function of purely imagi- 
nary argument and zeroth order. The function 
(3) is a regular potential, with g and dg/dz 
bounded inside the cylinder, if the constant X is 
restricted to real values. 

Substitution of (3) in (2) leads to 


i= -{f ade f e~ = To(Xa) n(z)dz, 
0 —x 


whence it follows 





* 1 
f z-*n(z)dz= — , (—-2<A< @). (4) 
2ral (da) 


—x 


Thus the integral equation (2) is reduced to one 
of the common Fourier type (4). Its solution can 
be written down immediately. By inversion of 
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(4), one obtains for the required charge density 


1 . dd 
f ei . ‘ 
4ra = To(\a) 


x 


= | 
7 2r’a* J, 


Starting from this integral expression, we can 
develop solutions more useful for numerical 
purposes. For non-vanishing 2 we may transform 
the above integral into a series by means of the 
conventional contour integration. In the cut 
plane z=0, however, stronger methods must be 
used, as will be seen below. 





* cos(z/a)l 


—dt. 
T(t) 


(5) 


Consider the following integral in the complex 
t-plane: 


1 etttla 


ni Jw To(t) 





dt, 
2ni 


the path of integration W leading from —R to 
R along the real axis, and back from R to —R 
along a semi-circle of radius R and center 0. 
Avoiding values of R for which the semi-circle 
passes through a zero of Jo(t), we can prove that 
for positive values of z, the contribution of the 
semi-circle tends to zero, when R tends to 
infinity. Let k,, ke, denote, in ascending 
order, the positive zeros of the Bessel function 
Jo(k). Then 72k, are the simple poles of the 
integrand in the upper half of the ¢-plane. The 
residue theorem now gives 


1 x eiztia - 
eae ——dt= 
2m1 —« T(t) 


e~knz/a 
+ pepmmomragey 
n=1 [9'(tk,,) 
We have Jo’ (ik,) =iJ;(k,), and so for the charge 
density 
1 2 e~knz/a 


n(z) = i ’ 
2a? n=1 J1(Ry) 





(s>0). (6) 


Only for z=0 does this series not converge in 
the ordinary sense. Unless z is very small, it is 
very useful for numerical computation. A small 
number of terms are already sufficient in the 
range 22a. 

The total charge, induced on the cylinder 
between 220 and , can be found by integra- 
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tion: 


a 


dra f n(ijde=-> 


z n= 


eh nz/a 


» (220). 
1R»J1(R») 


The latter series does converge at z=0. The total 
amount of induced charge on the upper half of 
the cylinder should be —3; we have thus an 
opportunity to verify a well-known formula in 
Bessel-function theory, originally credited to 
Nielsen : 

c 1 


a at ia 1 
a. 


} =} 
" tR,Ji(R») 


For small values of z the convergence of the 
series (6) is too slow for practical evaluation. 
Hence we have to develop other expressions for 
n(z), which do allow of numerical computation 
near the ‘‘cut plane’? z=0. This can be accom- 
plished by Watson’s method for an analogous 
problem.® 

Therefore we consider the following integral 

1 cosAw 


2m1 





-- dw 
w—t)Io(w) cos(rw/ 8) 

wherein \, 6, ¢ are real parameters; 6>0, 
cos(t/8) #0. The path of integration is a circle 
of radius R>+?, and center 0, not passing through 
any zero of the denominator of the integrand. 
When one adds the condition 8|\| <7, the inte- 
gral tends to zero, if R increases to infinite values. 
Thus the sum of the residues of the integrand 
is zero also. The poles of the integrand are all 
simple ; they lie at 


w=t, (a) 
w=+tk,, (n=1, 2, ---), (b) 
w=+(v+ >), (v=0, 1, 2, F -+). (c) 


The corresponding residues are easily found to 


be 








cosAt 
: (a’) 
I(t) cos(at/B) 
coshXk n 
- , (b’) 
J i(k») cosh(rk,/B) LR» +t | 
B (—1)v cos(Asvy+h8/2) (’) 
Bhi . 





x Io(8v+8/2)[(v+3)BFt] | 


5G. N. Watson, “The use of series of Bessel functions 
in problems connected with cylindrical wind-tunnels,”’ 
Proc. Roy. Soc. 130, 29 (1930/31). 
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Fic. 1. Curve showing the surface charge density || on 
the wall of the grounded cylinder, induced by a unit-point 
charge situated in the axis at z=0; z denotes the distance 
along the cylinder, measured from the source containing 
plane; a is the radius of the cylinder. 


The sum of all residues together being zero, one 
has 


cosMt 2k, coshrk , 1 


a _—_—— 
n=1 J,(k,) cosh(rk,/B) kp»? +0? 
(2v+1)8 


(Bv+B/2)?—-? 





Io(t) cos(at/ B) 
B «» (—1)’ cos(ABvy+A8/2) 
+—- 2 | 
Io(8v+8, 2) 





T v=0 


Evidently, the right-hand side of (7) is nothing 
else but the expansion in rational fractions of 
the left-hand side. 

We now multiply both sides of (7) by cos(t/8), 
thus neutralizing the poles in (c) 
from t=—2% to t=, 
formulas® 


, and integrating 
Using the following 


* cos(mt/B) T 
_o k,*+# ky, 








(8) 
J. cos(t/B) B) , (—1)’29r 
{=-_-—_, 
(av+8 +B/2)?—#2 (2v+1)6 
we easily obtain 
* cosM 2 CosrAs(v+ 4) 
[anne 
0 Io(t) v=0 Io(Bv+ 6/2) 
oe coshnXk », 
+2 5 — (9 


n=1 Ji(Rn)[1 + e2rkn/B] 
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Hence, for the induced charge density 
B «x cos(Bv+6/2)z/a 
—a?n(z) =— 
Io(Bv+ 8/2) 





27? »=0 
“ » coshk ,2/a 
a (10) 
oe n=1 Ti(Rn)[1e2ten/ 87 





It may be emphasized that our proof of (10) 
holds for B>0 and —2za/B<z<2xa/B. The first 
series on the right of (10) is convergent for all 
real values of z. It may be considered as an 
approximation of the integral in (5), obtained 
when the integration interval is divided into 
equal parts of length 8, and the integrand in 
each separate interval is replaced by its central 
value. The second series should then be inter- 
preted as a correction term. The latter, however, 
only converges for —22a/B<z<22a/B. By ana- 
lytic continuation, it can be shown that (10) 
holds for this extended range, twice as large as 
might be expected from the proof of (10), as 
given above. Moreover, for the special value 
B=ar/\|z|, the first series vanishes identically, 
whilst the ‘correction term’’ reduces to the 
series (6). Consequently, Eq. (10) is more basic 
than (6). 

Equation (10) is very suitable for numerical 
calculations. The constant £ is still arbitrary, 
except for —2mra/B<z<2za/8. Let us consider 
the extreme case z=0 in detail. We already saw 
that in this case (6) cannot be used on behalf of 
its divergence. Equation (10) can be used with 
any positive 6. We calculated —2z2°a’n(0) in 
seven decimals with B=1, B=3, respectively. In 
the first case we had to take eighteen terms of 
the first series and only one of the second. The 
result was found to be 


2.0832296 + 0.0000033 = 2.0832329. 


With 6=3 we had to consider six terms of the 
first, and three terms of the second series, 
obtaining in this case 


2.0053020 +-0.0779307 = 2.0832327, 
in close agreement with the earlier value. Thus 


we can state 


» dt 0.10554 
f = 2.083232---, = —. (11) 
o Io(t) a’ 








n(0) = — 
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Numerical computation of the surface charge 
density is now easily established. When 22a, 
already the first four terms of (6) are sufficient 


the correctness of the fourth 
decimal. Inserting known data from existing 


Bessel-function tables, (6) becomes 


to guarantee 


—a’n(z) =0.30657 exp( — 2.402/a) 
— 0.46774 exp(—5.52z, a) 
+0.58636 exp(—8.652/a) 


— 0.68466 exp(—11.79z/a). (12) 


In the remaining range 0< z<a, Eq. (10) may 
be used with 8=2. Then one obtains 


—a’n(z) = 0.08003 cos(z/a) +0.02076 cos(3z/a) 
+0.00372 cos(5z/a) +0.00060 cos(7z/a) 
+0.00009 cos(9z/a)+0.00001, cos(112/a) 


+0.00032 cosh(2.42/a). (13) 


The data of Fig. 1 were computed by means of 
Eqs. (12) and (13). The curve shows a rapid 
decrement of the surface charges with increasing 
distance from the cut plane. 

By expanding both sides of (9) into a power 
series of A, and identifying the corresponding 
coefficients, we obviously obtain 


va {?" - (By+B 2)?" 
f —dt=B >> - ‘ 
0 TI, (t) v=0 I)(Bv+8 2) 








wo R,*” 
+(—1)'2x > ov ARSE SEE 
n=1 J1(R,)[ 1 +e?" 5] 





The numerical values of these integrals can be 

found in the same way as shown above for r=0. 

Then for n(z) a power series in 2?/a? can be de- 
n I 

rived, convergent for |z| <a. The expression (13), 

however, is much more useful for numerical pur- 

poses in the corresponding range of z-values. 


3. FORMULAS FOR THE POTENTIAL AND THE 
ELECTRIC FIELD 


The potential V(P) can be considered as the ' 


sum of the potential 1/r, due to the unit charge 
_in the origin, and the potential Vo(P), due to the 
induced charges on the surface of the cylinder: 


1 1 n(¢) 
V(P) =-+ VoP)=—+ f —§do. 
r r ri 


In this r, 71, denote the distances from P to the 
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origin, and the surface element, respectively. 
Again, Vo(P) depends only upon z and p. 

At first we restrict ourselves to points on the 
axis of the cylinder. Then we have 


. * 9(x+u) 
Vo(z, 0) =2ra f ————_ (du. 
—« (u?+<a?)! 


Using (5), we obtain 


2 
Volz, Q)= — f 
T “oe 


On account of a well-known integral for the 
K-function, namely 


x x 


COSAZ COSAU 
inten anf ————-du. 
To(Xa). 0 (u?+a?)} 


x 


> Cosvu 
Kol da) = f — 
0 (u?+a?)} 





du, (14) 


we obtain for the potential of the induced 
charges, at points on the axis of the cylinder 


2 . Ko(Xa) 
Vo(s, 0) = — f cosdzdy. (15) 
Ty To(\a) 





In the rotationally symmetrical problem under 
consideration, the potential outside the axis is 
easily obtained, once the potential on the axis of 
symmetry has been found. We simply have to 
take the mean value of Vo(z+ip cosg, 0) over a 
full period of the angle g. Thus 


1 2r 
Volz, p) =— f Vol(z+ip cosy, O)dy. (16) 


To 


The integration over ¢ can be performed, leading 
to 





2 “ Ko(\a) ; 
Vo(z, p) = — f - To(\p) cosdzdd, (17) 
TV I o(da) 
and for the total potential we get 
1 2 * K(t) 
V(z, p) =————--— - f —-- 
(22+ p?)! Ta v4 T(t) 
X Ip(tp a) cos(/z ajdt. (18) 


The corresponding axial and radial components 
of the electric force are readily derived from 
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(18) by differentiation 


2 2 “ Ko(t) 
E, =———_——-_ - — ~ 
(s?+ p?)! ra" 0 I(t) 
XtIo(tp/a) sin(iz/a)dt, 

(19) 
p 2 * Ko(t) 

E,= ——+— f —tI,(tp ‘a) cos(tz ‘a)dt. 
(s*+ p?)} ra- 0 T(t) 


It does not seem possible to transform these 
integrals explicitly into known functions. 

Having thus found expressions for the electric 
field inside the cylinder, we can easily verify 
them afterwards. For any real value of \, the 
function 


cos(Az) - [o(Ap) 


is a regular solution of the potential equation 
inside the cylinder. Furthermore, the right-hand 
side of (17), being an absolutely convergent 
linear combination of those functions, is here a 
regular potential as well. Therefore, the function 
V of (18) is also a regular potential, with excep- 
tion of the origin, where it has just the required 
singularity. It only remains to be checked that 
V vanishes on the boundary at p=a. This, 
however, is evident from the Fourier inversion 
of (14): 


1 
(2*-+a*)! 





2 a 
f Ko(\a) cosdzd\ = 
Fo 


4. A SECOND METHOD TO SOLVE THE 
POTENTIAL PROBLEM‘ 


The first method to obtain expressions for the 
field inside the cylinder, as given above, started 
from a certain pre-calculated surface charge 
distribution on the boundary. This method is a 
little unconventional. Usually one works in the 
opposite direction: first the potential is obtained 
as the solution of a boundary value problem, and 
then the surface charges are calculated. The 
results above immediately indicate how this 
second method runs. One should try to find a 
function f(A) such that the potential of the yet 


®E. Weber, “The electrostatic field produced by a 
point charge in the axis of a cylinder,” J. App. Phys. 10, 
663 (1939). 
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unknown charges be expressible as follows 
Vals, 9) =f f0X)Ta(p) cosded, 


this being a continuous sum of elementary 
potential functions. Fortunately, the primary 
potential of the unit point charge can already 
be expressed in that way: 


1 2 " 
. J Ko(Xp) cosdzdh, (20) 
r To 


and, hence, for the total potential 
r F 2 ( 
V(z, p) -{ | f0)Ta(d0)-+-Kaldo cosAzdx. 
0 Tv 


A sufficient condition that V(z,p) be zero on 
the boundary of the cylinder is that the integrand 
vanish there ; consequently the function f(A) can 
be determined from 


2 
f( d)Io(Aa) +—K 9(da) a Q), 
T 


and then one obtains (17). 
Moreover, in combining Eqs. (17) and (20), 
we can derive the following expression 
2 *® Io(ha) Ko(Ap) — Ko(da) Lo(Ap) 
V(e0)=— f 
0 To(da) 





Tv 


Xcosrzdr, (21) 


which more directly shows the vanishing at p=a. 
Equation (21) enables us to come to the integral 
expression (5) for the charge density. Thereto- 
fore, we use (1) and the Wronskian of Jo, Ko: 


1 
Io(t) Ko’ (t) —Io'(t)Ko(t) = — a (22) 
t 


and so we are back at our former starting point. 


5. APPROXIMATE EXPRESSIONS FOR THE 
SURFACE CHARGES 


Weber,® in his theory of a point charge inside 
a cylinder, follows the second method apart from 
slight modifications of minor importance. His 
Eqs. (13), (17), and (18) are equivalent to ours 
as shown in Eqs. (18) and (19). Equation (5) for 
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the surface charge density, however, was not 
given by Weber. 

As for numerical calculations, Weber writes 
“the actual evaluation of the field distribution 
is rather difficult on account of the Fourier 
integrals which necessitate a point for point 
numerical integration in infinite limits,” and 
further “the numerical evaluation of the integrals 
is very laborious.’’ Weber, therefore, does not 
integrate point for point, but approximates the 
integrands by more tractable elementary func- 
tions, such that the remaining integrals can be 
computed explicitly. 

As a fair approximation, suggested by existing 
tables of Bessel functions, Weber uses, for in- 
stance 

K(t) t 


—tI,(t)=-e', (t>0). (23) 
I(t) 2 


From the lower integral in (19), Z, can then be 
evaluated at the surface of the cylinder. Weber’s 
result for the corresponding density of charge, 
according to n(z) = —E,/4x, may be written 


1 1 1 i-f 
—a*n(z) ~— 


Msasal’ ataarnel’ 
Ant (1+¢*)! we (1+ 7)? 
(¢=2/a). (23’) 


Furthermore, Weber derives expressions for 
the field in the cut plane as well for points at the 
surface of the cylinder, by another approxima- 
tion, different from (23). His approximation (21) 
reads 


Ko(t) 


; ~=(2r)'te-**, (t>0). (24) 
T(t) 


_ We need not, however, confine ourselves to 
the special values z=0 or p=a, as did Weber. 
When (24) is applied, we obtain 


1 ” K,(t) 
oa sf —JI(tp/a) cos(tz/a)dt 
(2x)! 0 I(t) 


x 


~Re f e*2+iD Jo(itp/a)tdt 
0 

2+7¢ 
[(2+i¢)*—p*/a?}? 


=~ 2. 
~~ © 
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TABLE I. Values of surface charge density || in units 
1/a?. First column: rigorous values; second: Weber’s 
approximation (23’). Third and fourth columns: approxi- 
mations (25’), (24’), respectively. 


z/a —a’n(z) 

0 0.1055 0.1049 0.1040 0.1285 
0.2 0.0980 0.0975 0.0969 0.1144 
0.4 0.0795 0.0795 0.0792 0.0827 
0.6 0.0581 0.0589 0.0586 0.0510 
0.8 0.0396 0.0413 0.0406 0.0280 
1.0 0.0259 0.0281 0.0269 0.0147 
®. 0.0165 0.0190 0.0174 0.0076 
1.4 0.0104 0.0129 0.0111 0.0041 
1.6 0.0065 0.0087 0.0069 0.0025 
1.8 0.0040 0.0060 ° 0.0042 0.0018 
2.0 0.0025 0.0041 0.0025 0.0016 


and then, from (18) 


1 1/8,\! 2+1¢ 
V(z, p) =—-——_-_ ( ) Re . 
(s?+p*)! atx [(2+7¢)?= p? a? ]} 





Especially on the cut plane, when differentiated 
to (—p), this formula yields 


1 8\} 6p*, a* | 
a |: en 
p? wr] [4—p?/a?} 2 | 


This result is apparently in disagreement with 
that of Weber. His relation reads 


1 | 8\? p*®/a*® ( 
E . Weber ~ — 1 + ( ) ee eee as = 0) ° 
' ol m/ [4—p?/a* }} 


The reason for this discrepancy is that Weber 
dropped one factor y in the integrand of his 
integral (22), so that actually he has used the 
following approximation 


K,(t) 
—— = (2m) *te-*# (25) 
To(t) 


instead of (24). 
In the latter case we find for general values of 
z and p 


1 1/8\3 
V(z, p) ~—_--(-) Re {(2+7¢)?— p?/a?}-* 
(2?+p?)! a\r 


leading exactly to Weber’s formula for E,. 
In either of the cases Eqs. (24) and (25), the 
corresponding surface charges, can be calcu- 
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a fh 








lated. We obtained, respectively, 


—a?n(z) = 


-| 1 

4ni (1+¢°)} 
8\} 

+(-) R 
Tv 

1 | 1 

4ni(1+¢°)3 


8\} 1 
+( ) Re- ; (25’) 
Us (1+76)3(3+79)! 


Table I shows the three different approxima- 
tions (23’), (24’), (25’), in presence of the correct 
values, as computed from (12), (13). The ap- 
proximation (24’) is obviously very poor, leading 
to values which are too large at z~0 and too 
small at s~2a. Either of the approximations 
(23’) and (25’) is excellent. That of Weber (23’) 
is a better approximation at small values of z 
than is (25’). On the other hand, for values of z 
exceeding a, the latter is much better than the 
former. Moreover, Weber’s formula yields nega- 
tive values for —a’n(z), when ¢>3.6. This is 
physically impossible, however. The surface 
charge density must always have the same sign. 
The right-hand side of (25’), on the contrary, is 
positive for large values of ¢. Therefore, as a 
whole, the approximation (25) may be considered 
as the most useful one, at least as far as we 
confine ourselves to the surface charges. 

The radial field at the surface of the cylinder, 
just opposite the point charge, is 1.3262 times 
as large as it would be under influence of the 
free-point charge, when the cylinder was re- 
moved. The corresponding factors in the three 
approximations (23), (25), and (24), are 1.3071, 
1.3183, and 1.6142, respectively. 





3(2+75) 
(24’) 


e 
(1+75)9°(3+75)9? 


—a*n(z) = 





6. FORMULAS WHICH ADMIT OF NUMERICAL 
CALCULATIONS INSIDE THE CYLINDER 


For exact numerical calculations the integral 
representation (18) is almost without any value. 
In Section 2 we derived expressions for the charge 
density which served its useful purpose. The 
field can be dealt with in a similar manner. 
The integral (15), however, is no suitable starting 
point, because its integrand shows a logarithmic 
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singularity in the lower limit. The difficulties 
owing to the many-valuedness of the K-function 
may be conveniently avoided in the following 
way. By partially integrating (15) and using 
(22), one readily obtains 


2 * sin(tz/a) 
Vo(z, 0) = -—f{ dt, 
0 


(26) 
TZ tI 9?(t) 


We now proceed in the same manner as we 
did before with (5). In this case the path of 
integration in the complex plane has to avoid 
the origin by a small indentation in the upper- 
half plane. We then obtain for z>0 


1 f” sin(tz/a) 

2 f greener 

rd, tlo*(t) 
We calculate the residue in t=ik, by putting 
t=x+1k,, developing the function into ascending 
powers of x, and taking the coefficient of 1/x in 
that expansion. The function is of the type 
f(x)/g?(x), where g(x) has a simple zero in x=0, 


and f(0) 40. The coefficient of 1/x turns out to 
be 


o eitzla 
t=3}+)= residue . 
n=1 (t=ikn) ¢192(t) 








(7 # so 
elf glee — aRn T12(Rn) - 


Consequently, for z>0 


1 22 e-bnle 
Vo(z, 0) = — +- p eae 
a n=1k,J1?(Rn) 


o 
~ 


To obtain the total potential on the axis of 
the cylinder we must add the potential of the 
unit charge in the origin. For positive values of 
z the latter is 1/z. Hence 

2 «@ e-knela 
V(z, 0) =- > — —. 
a n=1 kJ 17(Rn) 





Applying (16), we finally obtain for 2>0 


: 2 @ Jo(knp/a) 
| (z, p) =- , ® ogecg —~e—knz/a 
a n= 1 RJ 1?(Rn) 





(27) 


The series (27) converges also for x=0 with 
the only exception of the origin, where it diverges 
to infinity. In actual numerical work (27) appears 
to be useful for z/a 20.2; in the remaining range 
0<z/a<0.2 other methods should be applied. 
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One naturally asks whether or not it is possible 
to establish developments for the potential, 
analogous to that of formula (10) for the surface 
charge density. Actually, the required transfor- 
mation is readily performed for points in the 
axis of the cylinder. We only write down the 
final result: 


2 « sin(@v+8,2)z a 
V(z, 0) = -—-— + — 
rz »=0 (v+4)Io?(8v+8/2) 


11+e7-?**»/8! cosh(k,2/a) 
— (2ra/Bz) sinh(k,2/a) 


k,J:7(k») cosh?(rk,,/8) 








1 « 
--. 


a n=1 


(28) 


The first series converges for all real z, the second 
for |z| <2xa/B; B is an arbitrary positive num- 
ber. If especially 8=2z2a/z, the first series van- 
ishes identically, whereas the second one reduces 
to (27) for p=0. 

For 


actual numerical evaluation we have 


chosen 8=2. Then the series (28) becomes 


: a 
—aVo(z, 0) =—{0.79432 sin(z/a) 


~ 


+0.01782 sin(3z,a)+0.00034 sin(5z/a) 
+0.00000, sin(7z/a)+0.01014 sinh(2.42/a) } 
— 0.00323 cosh(2.4z/a). (29) 


Although convergent for z/a<z, (29) may be 


useful for z/a< 1. 

Unfortunately, there do not seem to exist 
simple developments for general non-zero values 
of p, because (16) yields for the continuation of 
(sinz/a)/(z/a) outside the axis 

a 


pla 
- f cos(tz/p)Io(t)dt, 
po 
which integral cannot be solved explicitly. 

For small values of z we can expand the initial 
‘integral (26) into powers of z/a. We then obtain 


2 





e ( = 1)’Ax Z ” 
Viz, 0)=—-- & (-) (30) 
avo (2v)! \a 
where 
1 2 ft?” 
he = ——_—— f —§dl, (31) 
(2v+1)2 0 Iy?(t) 
Ao = 0.4353382, 
de = 0.2059112, 
A, = 0.6588568. 
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TABLE II. Values of the potential (measured in units 
1/a) inside a grounded cylinder under induction of an 
axial positive unit-point charge; a=radius of the cylinder, 
z and p the usual cylindrical coordinates. Point charge 
in the origin. 





s/a\ 0 0.2 0.4 0 0.8 
0 2 4.125 1.612 0.756 0.303 
0.1 9.131 3.600 1.539 0.736 0.297 
0.2 4.137 2.670 1.357 0.681 0.280 
0.3 2.480 1.919 1.131 0.603 0.253 
0.4 1.661 1.395 0.916 0.518 0.224 
0.5 1.178 1.033 0.729 0.432 0.192 
0.6 0.864 0.776 0.575 0.355 0.161 
0.7 0.648 0.591 0.410 0.287 0.133 
0.8 0.493 0.454 0.355 0.231 0.109 
0.9 0.379 0.353 0.281 0.189 0.091 
1.0 0.292 0.272 0.219 0.147 0.071 
1.2 0.176 0.164 0.135 0.092 0.045 
1.4 0.107 0.101 0.083 0.057 0.028 
1.6 0.066 0.062 0.051 0.036 0.017 
1.8 0.041 0.039 0.032 0.022 0.011 
2.0 0.025 0.024 0.019 0.014 0.007 


Except Xo, these numerical values were given by 
Watson.°® 

We can now apply (16) with greater success 
than before. In performing the transformation 


(16) upon (30) we obtain 
Zp 
P-. =), (32) 
aa 


1 2r 
F,(u, v) =— J (u+iv cos¢y)*"dy 
0 


2n 


, 2 rs 
Vo(z, p)=-- 7 


a v=0 


(— 1)’Aoy 
(2v)! 





where 


ul 
=(u2-+0")Pa| —I, (33) 
(u?+v)! 


and P»2, is the usual symbol for the Legendre 
polynomial. For instance 


Fy= 1, 
F, =v —}v 
F.=u' — 3u?y? + 30. 


It may be 
2°+p* < 4a’. 

The development (32) becomes much more 
familiar when spherical coordinates r, #, are 
introduced, connected to the cylindrical ones by 
means of 


remarked that (32) converges for 


r=s+p"*; tand=p/s. 
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Then (32) is transformed into 


(2v)! 





Vo (z, p) = 


2 @ (—1)\»sr\”” 
(-) P»,(cosd). (34) 


a vw=0 a 


This series is convergent for r/a <2; it is nothing 
else but the development of Vo(P) into spherical 
potential functions. 

Especially in the cut plane z=0 we found 


—aV (0, p) =0.8707 +0.1030(p/a)? 
+0.0206(p/a)!+0.0057(p/a)*. 


The first three coefficients are calculated from 
the values Xo, A2, and 4, as specified above; the 
last coefficient is taken so that the total is 
1.0000, on account of V,(0, a) being —1/a, thus 
neglecting terms of higher order than (p/a)®. 

Numerical evaluation of the potential inside 
the cylinder may be readily accomplished now. 
The different expressions (27), (28), and (34) 
have reasonable overlapping regions of useful 
convergence. Numerical data for the total po- 
tential are to be found in Table II. 

The electric field itself can also be calculated 
numerically. From (27) we deduce for z>0 


2 « Jo(Rnp/a) 





e-=— -e~knzia 
a” n=1 Ji7(R») 
(35) 
2 £ Ji(Rnp a) 
E,=— F —_—_§_e- tole, 
a?n=1 J,7(R,) 


It may be noted that in Watson’s paper® two 
series S;, S. occur which are related to the series 
(35) by 


2 2 AS» 
E,=—S,;; E,=-—---—. 
a’ a Oz 





Therefore, Watson's series can be obtained from 
that shown in (27). Whereas Watson has given 
explicitly only the first four terms of the develop- 
ment corresponding to (32), we found the general 
term here expressible in spherical potential func- 
tions; cf. (34). With exclusion of the immediate 
neighborhood of the cut plane z=0, Eqs. (35) 
are very useful for numerical purposes. Around 
the cut plane we better use (34) and compute 
E,, Es instead of the cylindrical components; 
further details may be omitted, however. 
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7. A THIRD METHOD OF SOLUTION 


Either of the two distinct methods developed 
previously yields complete information con- 
cerning the field of the axial point charge in the 
cylinder. The mathematics there involved are 
quite instructive, and can serve successfully as a 
guide in problems of wave propagation inside 
the same cylinder. 

The reader is aware that even in the simplest 
potential problem, as was treated above, a con- 
siderable amount of analysis was used. In this 
respect wave guide theory is still less tractable. 
The question whether this is really necessary, 
from an engineer’s point of view, will be answered 
in due course. 

In practical wave guide applications, one is 
mainly interested in the fields far from the 
exciting antenna, where many of the higher 
modes have already been damped out. With 
this restriction in mind, and thus ignoring the 
field in the nearby zone, we have investigated 
whether it is possible to calculate the field in the 
far zone more directly, without the usual ex- 
tensive amount of complex function theory. This 
new method was found by inspection and general- 
ization of the expansion (27), obtained for the 
potential in the preceding section. Each separate 
term of that series, for non-vanishing z, is a 
solution of the potential equation, fulfilling the 
required boundary condition at the surface of 
the cylinder, and showing the suitable behavior 
at infinity. So we might have expected in 
advance the possibility of solving our potential 
problem by means of such a development, in 
terms of discrete normal solutions of the po- 
tential equation: 


V(z, p) =D CaJo(Rap/aye—*!#!/¢, 


n=1 


(36) 


When the proposed method actually works, 
it must be possible, somehow, to calculate the 
unknown coefficients c, in this development from 
the behavior of V(z, p) in the immediate neigh- 
borhood of the prescribed singularity at the 
origin, although just in the point charge, (36) is 
known to diverge. In the following we will show 
how this problem can be solved. We consider a 
more general problem, however, so that the new 
method is equally well applicable to solutions of 
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Moreover, we need not 
restrict ourselves to rotationally symmetrical 
fields. In many cases the prescribed singularity 
is of the type 


the wave equation. 


1 
td {2°+ p?—2ppo cos¢+ po?} 3, 


where R denotes the distance between the field 
point (z, p, ¢) and the source at (0, po, 0). Singu- 
larities of this character are commonly used, 
both in electrostatics and in wave guide theory. 
We may, however, admit still more general types 
of sources. 

Suppose G(z, p) is a given function with one or 
more singularities for z=0, OS pSa. Let G be 
even in z, and let it have continuous derivatives 
of the second order, outside the singularities. In 
the previous case of an axial point charge, 
G= (22+ p?)74. 

Furthermore, we suppose that f,(z) is a 
given system (n=1, 2, ---) of continuously 
differentiable functions with the property f,’(0) 
~0. For instance, f,(z) =exp(—k,2/a). 

The problem under consideration can now be 
stated in the following general form: How to 
find the coefficients c, in (37), when the following 
properties of the function F(z, p) are known: 

(i) F is an even function of z. 

(ii) F—G has continuous derivatives of the second order, 
both outside and inside the singularities of G. 


(iii) For z#0, F can be expanded into a uniformly (in p) 
convergent series of the form 


«o 


F(z, p) = = CnJm(Rn,mp/a)fn(| 2), (37) 
n=1 


where ky,» denotes the mth positive zero of the Bessel 
function J,,(k), and m is a given integer. 


The answer to this question is readily given if 
.appropriate use is made of the well-known 
Fourier-Bessel-Dini series theory, for details of 
which the reader may be referred to Watson’s 
standard work on Bessel functions.’ 

On account of the conditions stated under (ii) 
"we infer that (8/dz)(F—G) can be developed 
into a Fourier-Bessel series. This series is uni- 
formly convergent in any closed interval, con- 
tained in 0<p<a, for any value of z, whether 





7G. N. Watson, A Treatise on the Theory of Bessel 
Functions (Cambridge, 1944), Chapter XVIII. 
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outside or on the cut plane z=0. Thus we have 


0 x 
— => bn(2)Im(Rn, mp/a). 


n=1 


The coefficients },(z) can be determined in the 
conventional manner, using the orthogonality of 
the set of functions p)-Jn(Rn,mp/a) upon 0 pKa. 
One finds’ 


b,(2) _ 





2 °@ 
f p—(F—G)Jm(Rn, mp/a)dp. 
A? Jens 1?(R n,m) 0 Z 


_ 


As either of the two functions F, G is even in 2, 
with continuous derivatives, we have 6,(0)=0. 
For points above the cut plane (z>0) the integral 
may be broken up into two other integrals, 
namely 

2 0 


= — pF(z, p) I m(Rn, mp/a)dp 
7 T m417(R nw, mn) dz Jo 


2 ¢ 0G 
= bees pe agi al f p JulBautt a)dp. 
A7Insi'(Rnm) %o 92 





b,(z) 


From (37), which, for fixed z#0 is nothing else 
but the Fourier-Bessel series of the function 
F(z, p), we obtain 


f pF(z, Pp) J m( Rn. mp a)dp . 


0 
9 


a? 

=Cnf n( lz \—J m41?(R n,m) 
2 

Therefore, if z>0 


b (2) =Cnf n' (2) 


2 * 0G 
ee ee f p Jai Rn, mP, a)dp. 
a? J m4 s*(R no) 0 Oz 


Now let z tend to zero through positive values; 
then 6,(z)—-0, and f,’(z)—f,’(0)4#0. Conse- 
quently, the above integral must tend to a limit 
value also, and the coefficient c, is equal to 





2 ¢ 0G 

Cr= lim J p— 

a? fn’ (O)Fm41?(Rn, m) 20+ 0 Oz 
ym | ew a)dp. (38) 


The actual evaluation of the limit in (38) 
usually does not lead to difficulties. Let us take 
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as an example the potential problem in its 
simplest form, whereby the point charge is 
situated on the axis. Then (38) may be applied 
for m=0, G=(z?+ p?)-}. In the limit, the main 
contribution to the integral (38) is attributed to 
the values of the integrand in the neighborhood 
of p=0. Therefore we may substitute unit value 
for the Bessel function J». Consequently for the 
limit value 


© @a a f* pdp 
im f p—(o Ppt)-ldp=lim— f — 
0 0 


Oz Oz 2*+ p*)! 





0 
= lim— { (2?+a*)!—2} = —-1. 
Oz 
Furthermore, f,’(0) = —Rn,o/a= —k,/a; thus, 


2 


¢,= ites! ~s 
a*k nJi?(R al 


in agreement with (27). 

Until now we were only interested in bound- 
ary-value problems of the first kind, with vanish- 
ing V at the surface of the cylinder. In boundary- 
value problems of the second kind, the normal 
component of grad V vanishes there; this means 
0V./dp=0 in case of the cylinder. In this second 
problem the potential can also be developed in 
terms of Bessel functions. The difference is that 
the argument ky,»p/a is changed into kn, mp/d, 
where kn,» Means the mth positive zero of Jm’(x). 
Completely analogous to (37), (38), under similar 
conditions, we now have the following device: 

Given 


(2, p) = Do ¥nJm(Kn,mp/a)fn(|2z|), (37’) 


n=1 


‘behaves like’’ G(z, p) for z—0, then the coeffi- 
cients y, can be determined from 


2 
7 
a*f,’(0) 1 —m?, an iJa tan 





* @G 
Xlim f p—Jm(Kn,mp/a)dp. (38’) 
0 


z—0* Oz 
We would like to emphasize that in practical 
applications we may forget all the precise 


mathematical conditions and restrictions. Then 
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there only remains the very simple procedure to 
obtain the coefficients in the postulated expan- 
sions (37) and (37’) from the Eqs. (38) and (38’), 
respectively. In this connection the actual evalu- 
ation of those limits is of secondary importance. 
This third powerful method is incomplete, in 
so far that it does not very well provide numer- 
ical calculations in the nearby zone, because the 
convergence of the series (37) and (37’) is usually 
too slow there. This, however, might always be 
overcome by suitable extrapolation of the data 
outside the cut plane. On the other hand, it has 
many advantages. For instance, unlike the other 
two methods, it is equally well applicable to 
either the potential or the wave equation; the 
only difference occurs in the system of functions 
f(z). Nor does it matter whether or not there is 
symmetry around the axis of the cylinder. 


8. FIELD OF AN ECCENTRIC POINT CHARGE 


In case of an eccentric point charge, situated 
at (0, po, 0), the potential V(z, p, yg) will be ex- 
pressible as a Fourier series with respect to ¢, 
containing only cosine terms. Thus* 


wo 


@ 
V(z, p, ¢) = DX €mcosme > Cam 


m=0 n=1 
XI m(Rn,mp/aje*nmizi/e, (39) 


The coefficients Cn, m have to be determined from 
the knowledge that V behaves like (z?+> ? 
—2ppo cose+ po”) in the neighborhood of the 
source at (0, po,0). Obviously the method of 
Section 7 is applicable, if the functions G, F, are 
defined as follows: 





1 ad cosm gd y 
Geeo=—f — —, (40) 
24 J _, (27+ p?—2ppo cosy+ po’)! 
1 Tr 
F(z, p) =— J V(z, p, ¢) cosmed y 
2a J_, 
=D Ca, mI m(Rn, mp/a)e*amlel/e, (40’) 


n=1 


It remains to show how the limit in (38) can 
be actually evaluated in case of the special 
G-function (40). This can be accomplished by 








* eo=1; €, =2 if m>0. 
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means of a useful symbolic formula, namely 


lim(dG /dz) = —(1/ po) 5(p— po), (41) 


z-0" 


‘ 


where 6 denotes a delta-function, ‘‘vanishing 
everywhere, except for p=po, such that its 
integral over p from — ~ to + ~ has unit value.” 
As G is even in z, one has, outside the singularity, 
dG/dz—0 if z—-0; for p=po, this derivative tends 
to infinity. Therefore, the existence of something 
like (41) seems reasonable. Our final result is 
completely rigorous, although (41) has only a 
symbolic meaning. For further details the reader 
may be referred to the appendix. 

When (41) is applied to (38), we easily obtain 


2 Im(Rn,mpo/@) 


o 
Cn,m ’ 


a k,. ~~ » +17(R n, ma) 





and thus for the potential of the eccentric charge 
2« 

Viz, p, ¢) =— DY em cosme 
a m=0 


bed J m( Rn. mPo0 2) I m( Rn. mP a) 
x e-tnmisila, (42) 
n=l Raw ¥ mri (Ram) 





As it should be, this Green function of the first 
kind for the potential equation inside the cylin- 
der, is symmetrical in p and po. For po=0, (42) 
obviously reduces to (27). 

The method of the discrete normal solutions 
has thus proved its value. We could write the 
proper solution almost at once. Of course, this is 
not the only way; it is certainly possible to 
generalize the first or second method, treated in 
the preceding sections. For instance, one will 
find for the analog of (5) 


1 a 
‘n(z, ¢) = ——— €m COSM YE 


2r7a? m=0 
as I n(tpo a) 
x f ———— cos(tz/a)dt, 
0 I, (t) 


and by contour integration, analogous to (6), 
1 £ 


n(z, ¢) = ——— DL éncosme 
2ra’ m=0 
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The same expression for the surface charge 
density is obtained from (42) if (1) is applied to it. 

Similarly, other formulas for the axial point 
charge can be generalized to an eccentric charge. 
Further details may be omitted. We better 
demonstrate the usefulness of the third method 
by a problem of wave propagation. 


9. THE GREEN FUNCTIONS FOR THE WAVE 
EQUATION INSIDE THE CYLINDER. 
RADIATION OF ENERGY BY 
ACOUSTIC POINT 
SOURCES 


We restrict ourselves to harmonic vibrations 
of frequency w/2z, and time dependence 
exp(—iwt). The velocity potential of a free-point 
source is 


Uy=e"® R; R°=2?+p?—2ppyo coset+po. (43) 


Outside the point source, the velocity potential 
{’ is a solution of the wave equation 


AU+RU =0. (44) 


There are two different types of wave func- 
tions, according to whether U or dU /dp vanishes 
at the wall of the cylinder. The Green function 
of the first kind U is that solution of (44) which 
vanishes at the surface of the cylinder, which 
has at infinity the character of plane waves, 
coming from the direction of the origin of 
coordinates, and finally such that LU — U» remains 
finite for R-0. The Green function of the second 
kind U: has a vanishing 0/dp at p=a, the other 
requirements being the same as those for U,. 
Instead of the singularity (43) we may take 
merely 1/R. 

Discrete normal solutions of the wave equation 
(44) are 


cosm g: Jm( Rn, mp/a) 

Xexp[ — (2/a)(Rn,m?—k?a?)*], (45) 
cosm ¢+ Jm( Kn, mp/ a) 

Xexpl — (2/a) (Kn, m?— Ra?) ].  (45’) 


The functions (45) vanish at p=a, those of (45’) 
have a vanishing derivative with respect to p. 
On account of symmetry around the plane g=0, 
the potentials do not contain terms with sin(mg). 
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Therefore we may expect for 240 


D a 
uy = » Em cos g =: Cc n, — = mp a) 
n=l 


m=0U 


Xexp[—(/2) /a)(Rn,m?—k*a*)?], (46) 
Ms = > Em COSM Y > Fant alGe nf/@) 
m=0 n=1 
Xexpl—() x] /a)(kn,m?—R?a*)*].  (46’) 


In these cases the functions f,(z) are, respec- 
tively, 

exp[ = (z 2)(R n,m? —k?a*)?* |, 

exp[ —(z 2) (Kn, m?— k?a*)? |. 
It is seen that the only difference between (39) 
and (46) is in the occurrence of different func- 


tions f,(s). Consequently we can at once write 
down the result 


2 J m( Rn. mPo a) 
a (R» m? — k?a?)* Tm 17(R n,m) 





Can™ 


(47) 


Obviously the Green function U, for the wave 
equation reduces to that of the potential equation 
(cf. (42)) when & tends to zero. Concerning the 
square root, we must define for ka> Rn, m 


(Rn m? —k?a?)'= —i(k*a?—k,, »?)} 


because of the time factor exp(—‘of). 
For the second Green function we obtain by 
application of (38’) 


Z SaiSiaule a) 
Yxrnm>= aie (47’) 
a (1 —m- Kn. me) (Kn. me — R°2)' Tn? (Kn, m) 





It is to be noticed that in case of resonance the 
method fails, because then some coefficient c or 
y becomes infinitely large. This resonance phe- 
nomenon occurs for special values of the applied 
frequency, namely ka=ky, m Or RA=Kn,m, Tespec- 
tively, for some values of n, m. They correspond 
to the cut-off frequencies of the circular wave 
guide. In practice actual infinities do not occur, 
because of certain corrections which have been 
neglected in our treatment, and which only play 
an important part near the cut-off frequencies. 

Once the amplitudes of the different modes of 
propagation are known, it is not difficult to 
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compute the amount of energy, radiated by the 
point source into the cylinder. At large distances 
from the cut plane z=0, the potential is deter- 
mined by the undamped waves only. So we have 


> 


O <kn.m <ka 


Xexp[ —i(wt—(|z| /a)(k2a?—Rn, m?)*) J. 


Uy = EmCn, m COSM GT m(R n, mp/@) 


Because of the fact that the system of functions 
cosm oJ m(Rn, mp/a) 


is orthogonal over the cross section of the 
cylinder, one evidently has 


1 a 
hal ffimiras-Zlewal?f i 
Qn 0 


where the surface integral is extended over the 
cross section of the cylinder, and the summation 
extends over n,m such that Ry, m<ka. The inte- 
gration over p can be performed, and after 
inserting the values of Cx, m from (47), one obtains 


1 
- ff u, "dS 
Qn 


Tm? (k n, mPo ‘a) 








=2 > . (48) 
0 <ka.m <ka( ka? = k n, wt) J on-1*(R n, ™ 
Similarly, for the other problem 
1 
— ff U2|\*dS=2 > 
2r 0 <knim <ka 
Im? (Kn, mpo/@) 
. . (48’) 


x 
(R2a? — Kn, m2) [1 — (m?/ Kn, m2) |Im?(Kn, m) 


The physical meaning of the right-hand side 
of (48’) is that it provides us with the numeric 
by which the output energy of a freely radiating 
source should be multiplied in order to obtain 
the energy output, when the same source is 
placed inside the cylinder, at a distance po from 
the axis. 

By the method of images we can also treat 
the case of a point source inside a cylinder with 
one end closed. Let d denote the distance from 
the source to the closed end. Then, by intro- 
ducing an unclosed cylinder with two sources, 
a distance 2d apart, and vibrating out or in 
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phase, we obtain for the potentials 


u,= U,(z+d, p, ¢) — Ui(z—d, p, ¢), 
us= U,(z+d, p, ¢) + U2(z—d, p, ¢). 


The far field is a result of two interfering systems 
of plane waves. It is not difficult to obtain the 
correspondingly changed energy factors. Before 
in (48) and (48’) the summation is performed, 
each separate term must be multiplied by 
suitable interference factors; these factors are 








fd ] 

2 sin*| —(k2a?— Ry, m?)* I, (49) 
La ; 
rd a 

2 cos] —(k2a?— kn.m?)! I, (49’) 
La d 


respectively. The resulting expressions, which 
we need not write down in detail, will then 
provide us with the numerical factors, by which 
the output of fhe free source should be multiplied 
to obtain its output, when placed at distances d 
and pp from the closed end and the axis, respec- 
tively, inside a one-sided cylinder. 

Questions of wave propagation in electro- 
magnetic theory can be dealt with in quite the 
same way, further details of which will not be 
given at present. 


APPENDIX 


The validity of the symbolic formula (41) remains to be 
proved. That means, in rigorous mathematical terms, we 
have to prove 

: “ IGm 
lim f “ef(o)— "dp = —f(o0), (A) 
z—Orvo Oz 


where G,, is defined by 


<i" (B) 


G 1 f cosmed¢ 
Im = 7 —s > ee 3 rn 2 
—1(z*+p* —2ppo cosy + po’)! 


2a 
Equation (A) holds true for 0<po<a, if f(p) is continuous 
upon the closed interval (0, a), but otherwise arbitrary. 
For a complete proof of (A) it is sufficient to know that 
(A) is valid in the following two special cases: 
(1) f(p) =p”, where m is the same as in Gy; 
(2) f(p) vanishes at p=po. 
The most general case is obviously a linear combination of 
these special cases, because one can always write 
— S(o0) | 


f(p) ~- 
Po 


pe” +g(p), 


where g(p) vanishes at p= po. 
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Case 1 


For z>0, from well-known formulas in Bessel-function 
theory, we have 


(2? +p? — 2ppo cose + po")! 


= | e #t Jo[ t(p? —2ppo cose+ po")! dt 
Jo 


= De, cosve | e-** J, (pol) J, (pt)dt. 
/0 


0 


Upon multiplication by cos(mg), and integration over a 
full period of yg, we obtain 


Gn = | "e zt J, (pot) Im (pt)dt. (C) 
<0 


Consequently 
ta ‘ : 
p:p” (OG, /dz)dp 
“0 
a bleed 
=f p*ldp | —te zt J, (pot) Jm(pt)dt 
0 Jo 


woo it 
=— | f-=—le “In (potdt K@ rly, (x)dx 
/0 0 


—a™ wf e~ 2! In (pot) Im +1 (at)dt. 
Let z tend to zero; then the integral above tends to 

J Inout) Jn i(at)dt = (po a"™*!) (0<po<a). 
Hence 


lim J “p-p™(dGn 0z)dp = — po” (D) 
s-+0r" 0 


which is (A) for the special! function f(p) =p”. 


Case 2 


It is easy to demonstrate that the contribution of the 
intervals outside the immediate vicinity of po is zero; that 
means 


° *Po-é a ian bad a( Jun 
lim ff " Jo— do+f ef(e)— dp} =0, (E) 
z—0t le 0 Oz pore 02 J 
where ¢ is any positive constant, such that 0< po—e<po 
+e<a. For the proof of (E) we introduce M, the upper 
bound of |f(p)| upon (0, a). Then 

| fPo-* . OGm Po-€ |OGm 
Lf. fle) de| < Mf * |" dp. 
\70 Oz | v0 


| dz | 


Furthermore, 





0Gn ss & _— cosm ede ; ; 
Oz 2rd —«(2? +p? —2ppy cose +po")!’ 
thus 
| ° or : - 
|2 =<—/ ——— de = € 
| Oz | ~ 24d —2(z?+-p* —2ppo cose +po)! ~ [2? +(e —po)* 8 
Therefore, 
IG» | Z , 
|I——| <§— if —po| >e 
| dz | S(#+e)! oe 
e 


Consequently 


| (*Po-€ 9Gmn , | Mz Py 
| te a - vx ‘ 
J, ef(e) dz dp in @+eu, pdp 
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A similar bound is found for the contribution of the interval 
(po+e, a). Upon combining them we obtain 
Mza* 


| es . OGmn 
J, pf(e)-- dot”. efle) a 7 de) <  S 2(2 24 é)i 


Let now z tend to zero, e being kept constant. Then (E) 
follows at once. (It should be noted that (E) holds for 
general functions f(p), not necessarily vanishing at po.) 

For the special function f(p), with f(po) =0, the contri- 
bution of the remaining interval pp—e< p< pote is very 
small, because f(p) is very small there. Let us introduce 
the upper bound of |f(p)| in this interval, thus 

M(e)= Max |f(p)], 
|p —pol Ke 

then obviously M(«)—0 if e—~0. Furthermore, 
0Gn| > 2 (t* dg _  dGo 
“Oz ™ aJ x (s?-+-p? — 2ppo cos¢+po2)! Oz 





' 


Therefore, 


| ~Pot€ — OGm Pot! OG 
[-° of(e)"dp| < —M() {-* p“"dp 
Oz |  po-€ 3 


je po-€ 
— M0) fo ou, 


Let now z tend to zero; then the right-hand integral “— 
to —1, as follows from (D) when the latter formula is 
applied with m=0. 
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Thus we have 


tim fle) 7odp| < M(e); 
| 20+ 


that means, on account of (E) 
~M(e)< lim ioe dex M(6). 
2z—0** 


This is true for any e>0; as M(e)—>0 for e—0, conse- 


quently 
lim {7 


and this is the required expression in case 2. So we have 

rigorously proved the validity of the symbolic formula 

(41). It may be noticed here that according to Watson® 
8G. N. Watson, reference 7, page 389, Eq. (2). 

we can deduce from (C) the following interesting formula 


2+? +p0 
~~ ome (EF 2ppo ). ) 





= —f(p0) 


where Q denotes the Legendre function of the second kind. 
The Q-function has a singularity at the point where the 
argument is equal to 1; this corresponds with z=0, p=po. 


Eindhoven, March, 1946. 
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The Elastic Constants of Materials Loaded with Non-Rigid Fillers 


JANE M. DEWEY 
General Laboratories, United States Rubber Company, Passaic, New Jersey 
(Received November 14, 1946) 


General expressions are given for the modulus of rigidity, \2, and compressibility, &, of a 
medium of Lame’s constants, A; and 2, loaded with a volume fraction, ¢, of a filler of constants 
\,’ and i’. To terms in the first power of @ 


15(A2’ — 2) (A1 + 2A2) 
ha=he( 1+: 6) 
2d2'(3A1 + 8r2 )+-Ao(9A1+ 14A2) 
. 3+4)rok' 
k=k (k’—k)o. 
+ 34 dk ° 


If the “‘filler’’ is a gas at a pressure, p, in a nearly incompressible medium, these give 


h=E(1- ~ 2). 
9P+4E 
he =Ao(1— $¢). 
3 
3p+4r2" 


, 3E ) 
a=o( ~ Op+4E” ’ 


where ¢ is the volume loading, p is the pressure within the spherical cavities in the deformed 
state. Barred symbols refer to the properties of the loaded material; unbarred, to the medium 
alone. Expressions for the displacements and stresses within the medium and the particles, 


Young’s modulus 
Modulus of rigidity 


Compressibility k=k+ 


Poisson’s ratio 





neglecting interactions between particles, are also given. 


N a letter to the Editor of this Journal,’ the 

writer gave a solution of the linear equations 
describing small elastic deformations of a medium 
containing a single spherical inclusion of any 
elastic constants. This problem had previously 
been solved for simple tension by J. N. Goodier.’ 
The solution can be extended to the determina- 
tion of the elastic constants of a medium loaded 
with small spheres by the method applied by 
A. Einstein’ to the calculation of the viscosity 
of a suspension and by H. M. Smallwood‘ to 
the problem of rigid fillers. This method gives 
the terms in the first power of the volume 
loading only, and is valid, with neglect of surface 
effects, for very dilute suspensions.*:® 

1 Jane M. Dewey, J. App. Phys. 16, 55 (1945); Erratum, 
18, 132 (1947). 

2 J. N. Goodier, J. App. Mech. 1, 39 (1933). 

3 A. Einstein, Ann. d. Physik 19, 289 (1906) and 34, 591 
(1911). 

‘H. M. Smallwood, J. App. Phys. 15, 758 (1944). 

* The numerical value of the coefficient of the terms in 
the second power of the volume loading is given by Guth 
and Gold (reference 5), but no details of their calculation 
are available. 

5E, Guth and O. Gold, Phys. Rev. 53, 322 (1938) 
(abstract). See also E. Guth, J. App. Phys. 16, 20 (1945). 
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Consider a single sphere of radius R at the 
origin of coordinates in an isotropic medium in 
which the displacement, s, at a large distance 
from the sphere is given by 


s— Axi+Byj+Cszk. (1) 


TD 


Lame’s constants of the sphere will be denoted 
by Ay’ and 2’, of the medium by \; and do. 
As is the modulus of rigidity, \; may be defined 
by Ai+9A2=k", where k is the compressibility. 

The dilatation, strain, and radial components 


of the stress tensor, T, in the medium are 


V-s=A+B+C 
— 6d2(3A1+5dr2)~'c_or*#.M4 (8, @), (2) 
S>=3(r+b_er*) (A+B+C) 
+(r+c_or~*+c_ar*) Mi (8, o), = (3) 
so=([r+2d2(3A1+5A2)~!c_or? 
—§c_sr-*]M2(0, 6), (4) 
Sg =SoM2"(0, &) M3(8, >), (5) 
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T p= (Ar + 3Ad2—GA2b_or*) (A +B+C) 
+2r2[ 1 — (9A1+10d2)(3A1+5A2)7! 
Xc_or*—4c_ar* 1M (0, o) (6) 
T 9 = 2dof 1 + (31+ 2d2) 
x (31+ 5X2) ~!c_or-3 
+ §c_ar* |M.(0,¢), (7) 
T 6 = Te M2-'(0, 6) M3(8, $), (8) 
M,(0, 6) =A sin*@ cos*¢+B sin?6 sin*¢ 
+C cos*@—3(A+B+C), (9) 
M.(6, ¢) = = =sin#@ cosé 
06 


X(A cos*@+B sin*@—C), (10) 


M;(0, ¢) =——_- —— 
2sin@ d¢ 


= —(A—B) sin@ sing cos¢, (11) 
b_2= (31 +2A2—3A1’ —2)2’) 
X (4A42+3A1' +22’) R3, (12) 
c_2= —5(As! —As)(3A1 +52) 
x [22’ (3X1 +82.) 
+r2(9A1 +1422) FRY, (13) 
C-4=9(A2’ — 2) (Ar +Az2) 
X [22 (31 +822) 
+2(9A1+14A2) FIR (14) 


The energy of deformation is given by 


2Ww= f T,,-sda, (15) 


where T,, is the normal component of the stress 
tensor and the integration is over the outer 
surface of the medium. To calculate the energy 
of deformation of an infinite medium containing 
a single imbedded sphere, it is convenient to 
take a spherical surface with the particle at the 
center. Then, 
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W= 5 ( Ente + T 056 + T ,oSs)da, 


S 


=2V[\i(A+B+C)? 
2 
+2h2(A?+ B+ C) 1+ On Pe) 


x} [b_» — $A2(3A1+5A2)~!c_2 | 
xX (A +B+C)?—2d2c_2(3A1+5r2)7! 
x(A?+B+C)}, (16) 


for an infinite volume, where V is the total 
volume of the loaded medium. 

Assuming that to the first power of the volume 
loading the added energy due to N spheres is N 
times the added energy due to a single sphere 


W/V =3{d1+(A1— 9A2) 
X [b_2— 22(3A1+5d2)~!c_» |R-89} 
X(A+B+C)?+d2[1 — (A1— FA2) 
X (3A1+5d2)~'c_»R-9 ] 
x (A2+B?+C?), (17) 


where ¢=4$27NR*/V, the partial volume of the 
filler particles. 

The strain constants, A, B, and C of this 
expression are those for the medium, as appears 
from the derivation. The observed mean strain 
constants, Ao, Bo, and Co, can be calculated to 
the first approximation by the method outlined 
by Einstein. The displacement in the x direction 
can be written 


$,=Ax+) ; As,:=Aox, (18) 
where As,; is the added displacement due to the 
ith particle. 

If X; is the x coordinate of the center of the 


ith particle and x; is the x coordinate measured 


from X; 
x=xj,+X;. (19) 


Differentiating (18) 
0 0 ‘ 
Ayp=A+—); As,;=A+>D:—As,;. (20) 
Ox Ox; 


Assume that the added displacement is given 
by Eqs. (3) to (5) less the displacement in the 
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medium at a distance from a particle 


As2i=52i5—- AX =S,;—AXx;—-AX i. (21) 

Two assumptions are made here: the first, 
that at any point the added displacement due to 
each particle is the same as if the other particles 
were not present, and second, that the displace- 
ment at a distance from a particle is described 
by the boundary condition of Eq. (1). The first 
assumption should be correct to the first approxi- 
mation in the volume loading. The second is 
approximately correct if each particle is sur- 
rounded by a region free of particles of dimen- 
sions much greater than the particle radius. 
This is true of most of the particles if they are 
randomly distributed and the volume loading is 
much less than one. 

For randomly distributed particles }AX;=0 
and 


> Aszi= do 52; —-AX;j. (22) 


We replace the summation by N times the 
mean value of As,; 


0 OAs,,; x; 
>i —As,:=N vf —dV=N vf As,,—da, 
Ox; vy Ox; Ss Yi 
(23) 


where the first integral is over the volume and 
the second over the surface. To terms in R*r~?* 
the integral is independent of 7; and since the 
limits of angle are independent of 7, 


Ayp=A+N, vf As, sin8 cosd¢da, (24) 
Ss 


where As, may be taken for a particle at the 
origin. By and Cy are obtained by interchanging 
variables. 
Carrying out the integration 
Ao=(14+3(3A1+8A2)(3A1+-5A2)-!R-) A 
+ §b_2R-* — 3(3d1 + 8d2) 
X (3A1+5dr2)~'c_»R“* (A+B+C)¢, (25) 


whence, by combining with corresponding ex- 
pressions for By and Co, 
Aot+Bot+Co=(1+b_26)(A+B+C), (26) 
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Ao? + Bo? + Co? =[1+ $(3A1+8r2) 
X (3A1 +52) “1c_2R-% ] 
x (A?+ B?+ C?) 
+3(b-2—jc_2)R* 
X(A+B+C)*o. (27) 


The constants A; and A, of the suspension are 
obtained by equating the expression for the work 
per unit volume 


W, V=$A1(Ao+Bot+Co)? 
+h2(Ao?+Bo?+Co?) (28) 


to the expression (17). 
Substituting (26) and (27) in (17) and equating 
coefficients of (A9+By+Cpo)? 


Ai =A1— (Ar +: 2A2) 


X (b2— 2d2(3A1+5d2)~!c_2) RF. = (29) 
Equating coefficients of (Ao?+By?+ Co), 
Ae =Aeol 1 —3(Ar+2d2)(3A1 +5d2)~'c_2R-*9 | 
=o} 1+15(A2’ —Az) (Ar +2A2) 
* [22 (3A1+8A2) 
+A2(9A1+14A2) J-'d}. (30) 
Whence 
1 | _ it ‘ 
ga net ee 1h Ost Od (31) 
or 
k=k+(3+4d2k)(3+4d2k’)-'(k’ —k) db. (32) 
Two special cases are of importance: 
(a) The case of a rigid filler, 
2’ >Ao, and k’<Kk. 
Then, 
k=k[1—(1+ $r2k)¢ J. (32a) 
No= Aol 1+3(Ar1 +22) (Ar +$d2)—19]. (30a) 


(b) The case of included air or other gas 
do’ = 0, =p, 


where p is the pressure of the included gas 
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(after the deformation)* 
k=k+(3+4dck)(3p+4A2)—1(1—kp)d. (32b) 


The application of this formula to data on 
elastomers containing small amounts of gas is 
difficult as p is not known and 2, and # are of 
the same order of magnitude. If aged samples of 
reasonably permeable elastomers are studied, 
p can probably be taken as 1 atmosphere. In 
GR-I or other nearly impermeable materials, p, 
for very long times after molding, is probably 
approximately equal to the gas pressure at the 
time of molding. It may, therefore, differ con- 
siderably from 1 atmosphere. For any reasonably 
low gas pressure kp<1 and 


k~k+(3+4d2k)(3p+4d2)—'9, (32b’) 
h~dS 1 —$(Ar+2d2) (Ai + 45td2) 7! J. (30b) 


The approximate formulae for Young’s modu- 
lus, #, and Poisson’s ratio, ¢, valid for a medium 
of \;> de, containing gas at a low pressure are 


E~E[1—E(9p+4E)-'9], (33) 
é~o[1—3E(9p+4E)—'9]. (34) 


The percentage changes in the modulus of 
rigidity and Young’s modulus are of the same 
order of magnitude as the volume percent of gas 
in the material. The change in compressibility of 
a gum rubber stock caused by 0.01 percent of gas 
at atmospheric pressure is of the order of the 
compressibility, however. At pressures above an 
atmosphere, the compressibility varies rapidly 
with the pressure of the included gas. 

Within an embedded sphere and at its surface, 


* The equations and boundary conditions used in de- 
termining the displacements and stresses of Eqs. (2)-(8) 
are the conditions of equilibrium for a media obeying 
Hooke’s law and a linear relation between volume and 
pressure. When differential values of the constants are 
used, these conditions are valid for media of any type, since 
a linear relation may always be assumed for an infinitesimal 
deformation. For an ideal fluid \2=0 and dp= —di (dV/V), 
which reduces approximately to \;= p for a gas. 
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V-s=3ai, (35) 
S-= (€,M,(0, 6) +a,)r, (36) 
Se=cyrM (8, >), (37) 
S4=cywM;(8, >), (38) 

Tr = 2X2'C1M1(0, &) + (3d1/ +2d2’)a1, (39) 
T6=2d2'c1M2(8, >), (40) 
T 2 =2Xs'c1M3(8, ¢), (41) 


y= (M1 +22) (42+3A1/ +29")! 
X(A+B+C), (42) 

= 15do(A1 +22) [2A2’(3A1+8A2) 
+22o(9Ai+14A2) (43) 


As would be expected, the radial components 
of stress on the surface of a particle have their 
maximum value when the filler is rigid, the case 
treated by Smallwood.‘ 

Comparison of Eq. (30) for the modulus of 
rigidity with the calculation made by Taylor® of 
the viscosity of a suspension of fluid spheres of 
viscosity in a fluid of viscosity 


Su’ +2u 
a= ge +| 
2(u’ +n) 





shows that the two relations are of the same form 
only for Ay, Aq’, uw’, and d2’ infinite. The difference 
between the two cases lies in the importance of 
surface tension in the suspension of fluid drop- 
lets. Suspension of droplets of any fluid in another 
increases the viscosity, while blowing, for in- 
stance, reduces the rigidity of a solid. Added 
work on the fluid suspension is done against 
surface tension, deforming the droplets. Mathe- 
matically, the difference appears in the omission 
of the condition that the normal component of 
stress on the droplet be single valued. This 
difference in stress on the two sides of the surface 
is balanced by surface tension. 


6G. I. Taylor, Proc. Roy. Soc. A138, 41 (1932). 
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Heat Transfer between a Fluid and a Porous Solid Generating Heat! 


STuART R. BRINKLEY, JR? 
Central Experiment Station, U. S. Bureau of Mines, Pittsburgh, Pennsylvania 


(Received September 16, 1946) 


The theory, due to Anzelius, of the transfer of heat between a fluid flowing with constant 
velocity through a porous solid is extended to include the case where the solid generates heat. 
Expressions are obtained for the temperatures of solid and fluid as functions of position and 
time, it being assumed that the heat source function is a linear function of the temperature 
of the solid with coefficients independent of position and time. The application of the theory 
to a description of the temperature in catalytic reactions is indicated. 


INTRODUCTION 


THEORY of the transfer of heat between 

a porous body and a warm fluid has been 
developed by Anzelius* and Schumann.‘ A fluid 
carrying heat is supposed to flow with constant 
velocity through a porous substance such as a 
solid body in a finely divided state which is 
initially at a different temperature than that of 
the fluid. The temperatures of solid and fluid are 
determined for different positions at later times, 
it being assumed that the volume specific heats 
of solid and fluid are constant, and that there is 
no conduction of heat in either the solid or the 
fluid. The problem is, of course, idealized by 
these simplifying assumptions. However, in his 
extensive series of measurements of heat transfer 
from a gas stream to a bed of broken solids, 
Furnas’ obtained excellent experimental con- 
firmation of the simple theory. 

In the present communication, this theory is 
extended to the case where the solid is generating 
heat. In addition to the assumptions above, it is 
further assumed that the heat source function is 
a linear function of the temperature of the 
porous solid, the parameters of the linear func- 
tion being independent of position or time. 

The theory should be applicable to a cal- 
culation of solid and gas temperatures in cata- 
lytic reactions if the reaction vessel is so designed 
that all but a negligible amount of the heat of 


! Published by permission of the Director, Bureau of 
Mines, U. S. Department of the Interior. 

2 Physical chemist, Research and Development Division, 
Bureau of Mines, Pittsburgh, Pennsylvania. 

3 A. Anzelius, Zeits. f. ang. Math. u. Mech. 6, 291 (1926). 

*T. E. W. Schumann, J. Frank. Inst. 208, 405 (1929). 

5C. C. Furnas, U. S. Bureau of Mines Bulletin 361 
(1932). 
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reaction is transported by the sensible heat of the 
gas and if either (a) the reaction is of the zeroth 
other without volume change, or (b) the flux of 
the gas through the reaction vessel is so large 
that reaction is limited to a few percent. 


THE BASIC EQUATIONS 


Consider a semi-infinite (x20) porous solid 
through which a fluid is flowing in the positive 
x-direction with volume current g. Let the tem- 
peratures of solid and fluid at a distance x from 
the origin and at time ¢ be U(x, t) and V(x, 2), 
respectively. It is assumed that the solid in a 
particular element of volume loses heat by 
transfer to the fluid and gains heat by generation 
at a rate F per unit volume of porous solid, 
where in general, 


F=F(x,t; U), (1) 


and that the fluid in the element of volume gains 
heat by transfer from the solid and by the flux 
of heat into the volume which accompanies the 
flow of the fluid. By considering the rate of 
increase of heat in the solid and fluid contained 
between the transverse planes x and «+dx, and 
assuming that the rate of heat transfer is 
governed by Newton’s law of cooling, it is easy 
to show that the transfer of heat from solid to 
fluid is governed by the relations, 


aU 
Csps—t+h(U— V)— F(x, t; U)=0, 
ot 
(2) 


r 


OV 
Cypp-—+h(V — U)+cspjyg—=0, 
Ot Ox 


where c, and c; are the heat capacities per unit 
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mass of solid and fluid, respectively, p, and p; are 
masses of solid and gas contained in unit volume 
of porous solid, respectively, 4 is the coefficient 
of heat transfer between solid and gas per unit 
volume of porous solid, and where it has been 
assumed that the quantities c,p,, C;py, and g are 
constants. 

In the present communication, the solution of 
Kqs. (2) is obtained for the case where the heat 
source function is a linear function of the tem- 
perature of the solid, 


F(x, t; U)=h(a+Bsv), (3) 


with a and 8 constant. If the heat is generated by 
chemical reaction, one will usually have 


€ 
Fre-s ermewnrd 14 —(T- To)+ ee | 
. RT? 


0 


from which the accuracy of the linear approxima- 
tion, Eq. (3), may be estimated. Here, ¢€ is the 
activation energy of the reaction, R is the gas 
constant, and T is the absolute temperature. 

For convenience, the temperature of the inlet 
fluid is taken to be zero, and the temperature of 
the porous solid at the initial instant of time on 
the same temperature scale is denoted by U5. 
The boundary conditions for Eqs. (2) are there- 
fore 


U=Uo, x2¢t, (4) 


The first condition is a statement of the fact 
that the porous solid is at its initial temperature 
at points which have not yet been reached by the 
fluid. 

With the notation a=h/cypyg, b=h/c.p.g, and 
the change of independent variable 


gé=ax, 1r=b(gt—x), (5) 
Eqs. (2) and (3) become 


dU /dr=V+(8B—1)U+ae, 
dV /de=U-V, (6) 


with the boundary conditions, U=U(é, 7), 


V=V&, 7), 
U(é,0)=Uo, V(0, 7) =0. (7) 


Equations (6) and (7) lead immediately to a 
pair of ordinary differential equations for 


U(O, r) and V(é, 0). 
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d 
—U(0, r)+(1 — Bp) U(0, T) =a, U(0, 0) = Uo, 


aT 


d 
nm 0)+V(é,0)=Uo, V(0,0)=0. 


(8) 


The subsidiary boundary conditions are obtained 
by solution of Eq. (8). 


Q 
U(O, T) = Uce 8-)* -——_[e 8-1], 
1—, 


‘a . (9) 
V(é, 0) = Uo(1—e-). 


Equation (6) may be simplified in a manner 
analogous to that employed by Schumann‘ by 
introducing new independent variables, T and A, 
where 


P=(U+V)et or, 


(10) 
A=(U— V)et(-#)r, 
One obtains 
O1'/d7+0A/d07 = —A+2ae?* 1, 
ar /at—9A/at=P'+A, tt) 
and by further differentiation 
8A /atar=A. (12) 


The boundary conditions for [ and A are 
evidently 


r(é, 0) = Uo(2eE—1), A(E, 0) = Uo, 
(13) 


a 





r(0, 7) =A(E, 0) = Uot [e-#)r— 4 |, 


1—, 


The characteristics of Eq. (12) are the straight 
lines parallel to the axes through the point of 
interest (£, 7), and the boundary values of A are 
specified along a pair of characteristics (the 
axes). The solution of Eqs. (12) and (13) is thus 
a characteristic initial value problem and is 
straightforward by the method of Riemann.® 
The Reimann function is Jo[2((—&’)(7—7’))*], 
in which ~’ and 7’ are the variables of integration 
around the boundary of the region determined by 
the characteristics and the axes, and Jo is the 
zeroth-order Bessel function of the first kind for 


6 R. Courant and D. Hilbert, Methoden der Mathema- 
tischen Physik (Julius Springer, Berlin 1924), Vol. II, pp. 
311-317. 
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imaginary argument. The solution for A is 


A= Uol of 2(ér)! | 


taf e'-8)r’ of 2(E(r—7’))) ]dr’. (14) 
0 


This result may be rewritten in the convenient 
form, 


a 
A= Upl{ 2(ér)! | 4+—¢L (1 — 8), €/ (1 —8) J, (15) 
1 


in which the function g(x, y), defined by 


o(x, y) =e e~*’ Io. 2(x'y)! |dx’, (16) 
0 
has been introduced. This function, originated by 
Laplace,’ has been discussed by Thomas,* who 
has obtained the following useful properties: 
o(x, y) + ely, x) =e7+4— Tf 2(xy)!], 
(0/dx) ¢(x, y) = o(x, y) +Jo[.2(xy)!], 
(17) 
(0/dy) g(x, vy) = ¢(x, vy) —(0/dy) Io 2(xy)!], 








Tables of a related function, 


g(x, y) =e-*"-“ 9(x’, y*), (18) 


have been constructed® and are to be published. 
By means of Eqs. (17), Eq. (15) is readily shown 
to satisfy Eqs. (12) and (13), and this solution 
is unique.® 

can be 
obtained by the integration of either of Eqs. 


An expression for the quantity I 


(11), employing Eqs. (15) and (17). Expressions 
for the temperatures of solid and fluid can then 
be obtained by the combination of this result 
with the solution for the variable A, Eq. (15). 
However, the temperature of the fluid is obtained 
more directly from the relation, 
dV /d&=Ae*-"1-8) r (19) 


the substitution of the 
definition of A in the second of Eqs. (6). The 


which results from 


temperature of the solid is then obtained at once 


by the combination of this result with the 
¢(0,y)=0, g(x, 0)=e7—1. expression for A. 
THE TEMPERATURES OF SOLID AND FLUID 
If the solution for A be combined with Eq. (19), there is obtained 
OV oo a | 
—= } Uol of 2(ér)!]+ eLr(1—8), &/(1—8) Jpe FO. (20) 
ae 1—8 


Integrating this expression between limits as an ordinary differential equation, and recalling that 


V(0, r) =0, one obtains 


a 


g 
V = Uog(é, roe“ Peper f e [ r(1—8), ¢’/(1—8) ]dé’. 


1—p 


(21) 


0 


The integral of Eq. (21) can be evaluated by the method of integration by parts, employing the 


properties of g(x, y) which are summarized by Eqs. 


(17). The result of this process is 


_ 1B 
f e~*’ of r(1—8), &’/(1—B) Jd?’ =——-| 9 (E, r) —¢ LE/ (1-8), 7(1—B) J+ Le -® —eF Jem ® Je, (22) 
0 B 


if 80, and 


é 
J eo r(1—B), &/(1—B) Jd’ = [Eg(7, &)+7e(E, 7) — (Er) ViL2(Er)* Jers, 
0 





(23) 


7H. Bateman, Partial Differential Equations of Mathematical Physics (The Cambridge University Press, New York 


1932), p. 128. 
8’ H. C. Thomas, J. Am. Chem. Soc. 66, 1664 (1944). 


®S. R. Brinkley, Jr. and R. F. Brinkley, unpublished tables. 
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if 8=0. J, is the first-order Bessel function of the first kind for purely imaginary argument. Therefore, 
the temperature of the fluid is given by 


a a 
V(E, n= (Uet®) ole re-E-1-Or 4 e688) — gf E/(1—B), r(1—B) Je O-r—1}, (24) 
B 


if 8#0, and by 

V(E, 7) = {| (Uotar) 9(&, 7) +akg(r, §) —a(ér)i[2(Er) Jje*, (25) 
if 8=0. The temperature of the solid is immediately obtained by combination of these results with 
the expression for A. Then, 


aj | 


a 
ULE, N= (vot Jie ol, seo} —— | of r(1—B), E/(1—B) eT, (26) 
gli—p 


B 
if 840, and 


U(E, 1) — Uo= | Lalé+1) — Uole(r, #)+2r¢(E, 7) alr) WiL2(Er) Je, | (27) 


if 8=0. The rate of heat transfer from solid to fluid per unit volume of porous solid is equal to 
hA(é, r)c&-O-#)r, 


The application of these results to particular systems will be the subject of subsequent reports 
from this laboratory. 


THE STEADY STATE 


The steady state values of the temperatures of fluid and solid are obtained from Eqs. (24) to (27) 
by proceeding to the limit of infinite 7. The function g(, ©) is finite for finite £, and the function 
y(t, &) tends to e§*" as 7 tends to infinity. Then, one obtains 


a 


V(E, ©) = (a/B) (eFC — 1}, UE, @) =— 





eft! (1-8) 1} (28) 
aa 


V(é&, ©)=aé Ul, ~)=a(E+1), (29) 


if 80, and 


if 8=0. These steady state solutions have been previously obtained and employed for the calculation 
of catalyst and gas temperatures in a typical operation of the Michael process for the synthesis of 
hydrocarbons from carbon monoxide and hydrogen." 


S. R. Brinkley, Jr., 110th Meeting, American Chemical Society, Chicago, Illinois, Sept. 1946. 


VOLUME 18, JUNE, 1947 585 








Letters to the Editor 








Concerning Estimates of the Minimum Sparking 
Potential Based upon the Cathode 
Work Function 


DonaLD H. HALE AND W. S. HuxForpD 
Department of Physics, Northwestern U niversity, Evanston, Illinois 
April 4, 1947 


E prrew and LaRocque! have recently suggested a 
method for estimating the minimum sparking potential 
for cold cathode tubes. They measured the minimum 
sparking potentials for several cathodes in argon and 
calculated the corresponding value for y,,. These values for 
ym were then plotted against the corresponding values of 
the work functions and the minimum sparking potentials. 
By use of these two curves an attempt is made to estimate 
ym and the minimum sparking potential in argon for any 
cathode, provided only the work function for the cathode 
material is known. Applying this method the authors 
obtain an approximate and reasonable value for the 
minimum sparking potential of a sodium cathode in argon. 

This method of determining y, from sparking data is 
not reliable. It has been pointed out that this procedure 
cannot be expected to yield good values for y in all cases.” 
This is shown by the results obtained by these authors 
since their values of y vary widely from those found by 
other workers. Schéfer? and Engstrom‘ found that the 
values of y for a barium cathode in argon increased with 
increasing values of E/p starting with an E/p of about 40. 
The authors found just the opposite. Their values of + 
decrease with increasing values of E/p. Further, their 
value of y at an E/p of about 40 is three times greater 
than that measured by Schéfer. Since they worked at 
higher gas pressures than did earlier workers it may be 
that in their work the mechanism of electron emission by 
photons at the cathode was enhanced. However it is 
difficult to see how this could account entirely for their 
high values of y. 

The agreement between the calculated and experimental 
values of the minimum sparking potential for a sodium 
cathode in argon may be fortuitous. Bowls® measured + 
for a sodium cathode in nitrogen and Ehrenkrantz® meas- 
ured the sparking potentials. Bowls found that the sodium 
cathode gave distinctly lower values of y, for a wide range 
of values of E/p, than did a platinum cathode. Yet, 
Ehrenkrantz found that the sodium cathode showed a 
lower minimum sparking potential than did the platinum. 
Hale’ and Ehrenkrantz found a similar situation with 
regard to sodium cathodes in hydrogen. The authors 
explain the lower y for sodium on the grounds that hydro- 
gen or nitrogen combines with sodium to give a surface of 
high work function. Such an assumption would explain 
the low values for y but it cannot explain, on the basis of 
the author’s argument, the observed low value of the 
minimum sparking potential. Actually it has been pointed 
out that sodium surfaces in hydrogen or nitrogen show 
low photoelectric thresholds.’ Here then is a case of a 
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cathode material which has a low work function, which 
shows low values of y, in at least two gases, and yet has a 
low minimum sparking potential. 

There are other interesting examples. The work functions 
for nickel and platinum are 5.01 and 6.30 volts, respec- 
tively.’ Yet the minimum sparking potential for argon 
with cathodes of these materials is 195 volts in each case.!° 
The work function for iron is 4.72 volts.* However, the 
minimum sparking potential for argon with an iron cathode 
is 265 volts. This is distinctly higher than the minimum 
sparking potential for platinum in argon. Yet platinum has 
the higher work function. 

It is probably correct to say that, in a very general way, 
materials with low work functions will show high values 
of y. However, exceptions are known and any prediction 
of the minimum sparking potential which is based upon 
the work function of the cathode may be badly in error. 
Much work remains to be done before this point is clear. 

We wish to call attention to what appears to be a mis- 
statement. The authors state that Bowls found a lower y 
for platinum in a large region of E/p than he did for 
sodium. It appears that they intended to refer to a lower y 
for sodium which was what Bowls observed. 

1H. Jacobs and A. P. LaRocque, J. App. Phys. 18, 199 (1947). 

2?L. B. Loeb, Fundamental Processes of Electrica | Discharge in Gases, 
pp. 348 and 416. 

3 R. Schéfer, Zeits. f. Physik 110, 21 (1938). 

*R. W. Engstrom, Phys. Rev. 55, 239 (1939). 

5 W. E. Bowls, Phys. Rev. 55, 293 (1938). 

®F. Ehrenkrantz, Phys. Rev. 55, 219 (1939), 

7D. H. Hale, Phys. Rev. 55, 815 (1939), 

5 L. B. Loeb, reference 2, p. 395. 

» A. L. Hughes and L. A. DuBridge, Photoelectric Phenomena, p. 75 


10M. J. Druvesteyn and F. M. Penning, Rev. Mod. Phys. 12, 114 
(1940). 





Measurement of the Dynamic Stretch-Modulus 
and Hysteresis of Tire Cords 


W. JAMES Lyons AND IRVEN B. PRETTYMAN 
The Firestone Tire and Rubber Company, Akron 17, Ohio 
April 17, 1947 


HE modulus of elasticity (Young’s) of synthetic fila- 

ments and films has been measured by Meyer and 
Lotmar! and Ballou and Silverman,? who employed sound- 
velocity methods which gave the dynamic value of the 
modulus at acoustical frequencies. These methods, how- 
ever, are not applicable to the measurement of concomitant 
hysteretic losses. Wakeham Honold’® studied the 
elastic modulus and hysteresis effects in cotton, rayon, 
and Nylon tire cords under dynamic conditions, but at a 


and 


relatively low frequency: 1 cycle/sec. 
The exploratory studies reported herein were undertaken 
to develop a method for measuring not only the dynamic 


250v.DC. 
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Fic. 1. Schematic diagram of cord stretch-vibrometer. 
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TABLE I. Representative measurements of dynamic properties of conventional tire cords, obtained on cord vibrometer 


at constant impressed force-amplitude. 


Static properties 
(conditioned at 65% 
R.H., 70°F) 


Breaking Ultimate 

Cord sample load elongation 

Fiber Construction (Kg) (Percent) 
Cotton 11/4/2 8.5 14.3 
Viscose rayon 1100/2 6.8 15.9 
Nylon 210/3/3 14.6 19.5 


modulus of filaments, yarns, and cords in the sonic or near- 
sonic range, but also their internal frictional (hysteretic) 
properties. Measurements of the modulus and _ friction 
have been successfully made on cotton, rayon, and Nylon 
cords forcibly vibrated longitudinally at frequencies above 
100 cycles/sec., by electromagnetic means. 

The apparatus is shown schematically in Fig. 1. A small 
light-weight coil is placed coaxial with the core of a strong 
electromagnetic field coil, so that the former is in a powerful 
magnetic field when a direct current passes through the 
field coil. The small coil is rigidly mounted on a plastic 
rod which passes through a hole at the axis of the core. 
The rod is supported by pieces of the test cord attached co- 
axially at each end. Static loads at the outer ends of the 
two pieces of cord balance each other and allow the cords 
to be placed under tension without displacing the rod and 
small coil. When an oscillating current is sent through this 
coil, the rod vibrates longitudinally, causing small periodic 
extensions and contractions in the pieces of test cord. In 
the initial trials, a small mirror was rotatably connected to 
the rod at the end opposite the vibrating coil, so that by 
means of an optical-lever system, amplitude of vibration 
could be detected and measured on a ground-glass scale. 
It was subsequently found that a micrometer microscope 
could be advantageously substituted for the optical-lever 
system, with improved accuracy; and this method was 
used to obtain the results reported herein. In the present 
experiments the vibrating coil was energized by a 50-watt 
power amplifier which was controlled by a regular com- 
merical calibrated audiofrequency oscillator unit. 

The principle underlying the present method is ex- 
pressed by the equation of motion of the vibrating element 


under impressed force: ‘ 
m(d?s/dt®)+-(n/q)(ds/dt)+(E/q)s = F coset, (1) 
where m=vibrating mass, s=displacement from the 


equilibrium position; »=coefficient of internal friction 
(hysteresis); E =dynamic stretch modulus; g = shape factor 
of cord; and F=amplitude of applied force having fre- 
quency w/2m. Equation (1) has been applied to forced 
compressive vibrations in flat cylinders of rubber by 
Gehman, Woodford, and Stambaugh,‘ and to shear vibra- 
tions in similar specimens by Dillon, Prettyman, and 
Hall.® The design of the driving mechanisms used by these 
groups provided the prototype for the corresponding 
elements of the present equipment. 

It can be shown readily that at resonance, to good 
approximation: 


E=mqwo = m(L/2A)(2xrfo), (2) 
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Dynamic properties 
(at approximately 14% R.H., 78°F) 
Dynamic Internal Hysteresis 
modulus, E friction, 7 index, 2x/fn 
(101° dynes/cm?) (10° poises) (10!° dynes/cm*) 


_ Resonant 
frequency, fo 
(cycles/sec.) 








169 rE 4.0 0.42 
148 11.2 $.3 0.49 
116 6.1 2.4 0.17 








where A =cross-sectional area of cord or yarn (exclusive 
of interstices between fibers); L=free length of each 
matched piece of cord, and fo=resonant frequency. 

Equation (2) has been applied to observations of reso- 
nant frequencies made with the cord stretch-vibrometer 
on three kinds of tire cord, with the results summarized 
in Table I. To identify the cords the conventional static 
properties have been included in the table. All cord samples 
were vibrated under a tensional load of 1.2 kg, with 
L=25.4 cm (10 in.). The effective mass of the vibrating 
system, m, was found to be 17.0 g. The cross-sectional area 
of each sample was determined from the weight of a 200-cm 
specimen of the cord and the bulk density of the component 
fiber. 

Internal friction, 7, may be readily derived from obser- 

vations on the stretch-vibrometer by means of the equa- 
tion: 
(3) 
the maximum (resonance) strain-amplitude. 
Values of this quantity, as well as of an index of hysteresis 
losses 27fn, are entered in Table I. The latter quantity 
has been shown, in the case of rubber, to be substantially 
independent of frequency. 


y= Fq/Smo, 


Sm being 


1K. H. Meyer and W. Lotmar, Helv. Chim. Acta 19, 68 (1936). 

2J. W. Ballou and S. Silverman, Textile Research 14, 282 (1944), 
J. Acoust. Soc. Am. 16, 113 (1944). 

3H. Wakeham and E. Honold, J. App. Phys. 17, 698 (1946). 

4S. D. Gehman, D. E. Woodford, and R. B. Stambaugh, Ind. Eng. 
Chem. 33, 1032 (1941). 

5J. H. Dillon, I. B. Prettyman, and G. L. Hall, J. App. Phys. 15, 
309 (1944). 





On the Analysis of Internal Flow Machinery 


SABINUS H. CHRISTENSEN 
Carrier Corporation, Syracuse, New York 
April 11, 1947 


HIS is a brief report of a new method of applied re- 

search which has opened an area of analysis hitherto 
inaccessible to internal flow machinery having moving 
passages (centrifugal and axial-flow compressors, and 
turbines). 

The new technique is based on the rejection of the 
idealization in the gas-dynamic analysis of internal flow 
machinery that the fluid may be considered to be a com- 
pressible ideal gas, i.e. frictionless. It has been found by 
exhaustive analytic study of test data that the use of the 
frictionless gas concept fails to explain performance within 
very wide limits. 

It is well known that the inefficiency of internal flow 
machinery is caused by the dissipation of available energy 


587 














which is reflected in increase of entropy through the 
machine. The basic domain for research on internal flow 
machinery is then recognized to be that of making a de- 
tailed explanation of entropy increase caused by friction, 
separation, and turbulence. Taking the detailed change of 
entropy into full account now explains such hitherto 
inaccessible phenomena as quantitative values of shock 
loss at the entrance to rotating grids, and the theory of 
the maximum flow in a centrifugal compressor (as have 
been done by the writer in unpublished company reports). 

The new method is directed towards the local evaluation 
of entropy in rapidly moving or rotating passages. The 
method utilizes the device of speeding up the moving 
passage (or the equally effective method of using a very 
heavy gas with slower passage speeds) until somewhere 
within the passage the relative velocity is locally sonic. 
The exact equations of state, energy, and continuity are 
then satisfied by a definite entropy value at the blocking 
point. If a study is made of the energy vs. area gradient 
along the passage it will be seen that the blocking point 
is fixed by the passage design, thus providing the fourth 
simultaneous equation which solves the problem of deter- 
mining the exact thermodynamic state inside a rapidly 
moving passage. It is believed that this technique is far 
superior to the less reliable method of trying to measure 
the state of the fluid by instruments which are forced 
to transmit their data through commutator rings and 
packings. 

Instead of working with entropy change, re-arrangement 
of the gas-dynamic equations shows that the isentropic 
relations apply if the passage area is modified by the factor 
exp(—JAS/R), J being Joule’s mechanical equivalent of 
heat, AS the increase in entropy, and R the gas constant. 





Objective Aperture System for the 
Electron Microscope 


Cecit E. HALL 
Massachusetts Institute of Technology, Cambridge, Massachusetts 
April 17, 1947 


HE conventional objective aperture used in magnetic 

electron microscopes consists of a circular diaphragm 
a»out 25 microns in diameter, fixed with respect to the 
lens and coaxial with it. It prevents electrons scattered 
through angles larger than about 10~* radian from reaching 
the image plane, thereby reducing background intensity 
and enhancing contrast in the image. An objective aperture 
is necessary for the examination of thick objects producing 
a high scattered intensity, such as replicas, sections, and 
the like, but it is of lesser value in the examination of 
specimens of suitable thinness and contrast. Electron 
‘microscopists are familiar with the many practical diff- 
culties associated with the use of apertures. They are 
difficult to make, install, and align, and must be replaced 
frequently, because they accumulate films of non-conduct- 
ing material which become charged and produce distortions 
in the image. It is therefore desirable to have an objective 
aperture system that can be installed quickly and renewed 
as occasion demands. 
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Fic. 1, Objective aper- 
ture system. 














An objective aperture system having the advantage 
that it is attached to the specimen cartridge rather than to 
the lens, is shown in Fig. 1. It has been used successfully 
in this laboratory with an RCA type B electron microscope. 
The device consists essentially of a second screen mounted 
and aligned so that each opening of the specimen screen 
has below it a corresponding opening serving as an aperture 
stop. Aperture screens may be detached readily, if they are 
not required for the specimen under examination, and may 
be renewed frequently from standard material. A is the 
lower part of the specimen cartridge to which is fitted the 
cap, B, containing the specimen. The tip, C, containing the 
aperture screen is mounted in a friction cup so that it may 
be turned about the axis and rocked slightly to give small 
lateral displacements. Aperture screens, which must be 
non-magnetic since they are in the lens field, are punched 
from 200-per-inch copper Lektromesh.! 

The procedure in using the device is as follows: Clip the 
cartridge to the stage of a 100-power optical microscope 
with cross-hairs in the eyepiece. Focus on the object screen 
and rotate the eyepiece until the cross-hairs are aligned 
with the screen wires. Focus on the aperture screen and 
rotate the tip C until the screen wires are square with the 
cross-hairs. Refocus on the object and center an open 
square on the cross-hairs with the mechanical stage. Focus 
on the aperture mesh and rock the tip, C, until an open 
square is centered on the cross-hair intersection. This 
procedure usually takes one or two minutes. Possibly a 
smaller object-to-aperture distance than the three milli- 


Fic. 2. Collodion-silicon replica 
of glass grating, shadowed with 
chromium; (a) without an aperture, 
(b) with mesh aperture. 3600. 
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Fic. 3. Dark-field electron micrograph of carbon particles 
shadowed with uranium. 11,666. 


meters shown in the diagram would provide adequate 
contrast in the image and at the same time simplify the 
alignment of the two meshes. 

Two micrographs of a collodion-silicon replica of a glass 
grating are shown in Fig. 2, to demonstrate the difference 
in contrast with and without the aperture screen; (a) was 
taken without an aperture, while (b) was taken with one. 
The replica was shadowed with chromium. Contacts were 
made from the original micrographs on the same Eastman 
Medium lantern slide plate and a contact print made from 
this on F-3 Kodabromide. Contrast, of course, depends on 
the photographic processing, but an effort was made to 
have the prints represent a reasonable comparison of the 
original plates. Although improved contrast in the image 
is most pronounced with thick specimens, improvement has 
also been noted in thin specimens containing small par- 
ticles of low scattering power. Open areas in standard 
Lektromesh screens are square, but no distortions have 
been noted in the image which could be attributed to this 
unusual shape. Also, it is possible that electrons may be 
so widely scattered that they will reach the image plane 
through an opening other than the one directly below the 
point under observation. In the mesh normally employed 
the solid parts are three times as wide as the open squares, 
and since scattered intensity falls off rapidly with angle, 
the proportion of such widely scattered electrons appears 
to be inappreciable. 

The multiple aperture was originally devised to enhance 
contrast in images of thick specimens, but it has proved 
to be very useful also for dark-field electron microscopy. 
If alignment of the two meshes is not perfect, and it very 
rarely is, a part of each object opening will be seen in 
dark field. An example of a dark-field micrograph made in 
this manner is reproduced in Fig. 3. The specimen consisted 
of carbon particles? on a collodion film, shadowed with 
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uranium to a calculated thickness of about 13A. Shadow 
areas are dark because they contain no uranium, hence 
scatter least. The uranium film appears as a background of 
discrete spots ranging in diameter from 50 to 200A. The 
spots are quite well resolved for dark field, indicating that 
they are produced by electrons that have been scattered 
coherently and are therefore confined to a narrow angular 
aperture and velocity range. The dark-field micrograph 
thus reveals the presence of small ordered regions in the 
film that are not perceptible as such in bright field. Similar 
structure has also been observed with thin films of several 
other metals and compounds. 

Images of the carbon particles shown in Fig. 3 are charac- 
terized by bright contours and dark centers, resulting from 
the fact that at small angles there is a relatively high in- 
tensity from electrons scattered close to the surface. Elec- 
trons scattered within the particles are spread over a 
much larger solid angle so that a smaller proportion enter 
the limiting aperture. An interesting feature, not visible 
in bright field, is the occurrence of numerous bright spots 
at the edges of the particles. This appearance is not associ- 
ated with the uranium since the spots also appear in un- 
shadowed specimens. In the micrograph shown, all such 
spots appear on edges away from the aperture, hence 
represent pencils that have been deflected under the 
particles. However, the occurrence on one edge only, is 
dependent on alignment of the electron microscope and on 
the objective current. They have been observed on edges 
toward or away from the aperture and on all sides simul- 
taneously. The reason for enhanced scattered intensities 
is not evident at present. Possibly they are in the nature of 
crystalline reflections owing to regularities near the surface, 
although the images are considerably more diffuse than 
those produced by diffractions from common crystalline 
substances. Whether the enhanced scattering is to be 
accounted for on the basis of crystallinity or not, it must 
be noted that near the edges, electrons are incident and 
emergent at such small angles that refraction effects may 
be expected to be of significant magnitude. The observed 
deflection of electrons under the particles could be ac- 
counted for on the basis of refraction resulting from an 
inner potential of the order of 10 volts. 

These observations demonstrate the utility of the mul- 
tiple aperture which is, to be sure, a compromise for con- 
venience rather than the equivalent of a clean, well- 
centered aperture of the conventional type. Probably a 
completely satisfactory solution of the aperture problem 
will require the redesign of present objectives to allow for 
insertion and alignment of apertures fixed to the lens, 
without the difficulties and hazards that presently attend 
this operation. As an alternative, the multiple aperture 
has some advantages, and furthermore, the ease with which 
it may be used for dark-field observations may stimulate 
this interesting and somewhat neglected aspect of electron 
microscopy. 


1C. O. Jelliff Manufacturing Corporation, Southport, Connecticut. 
2 P33, Thermatomic Carbon Company, Inc. 
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New Books 








A Shorter History of Science 
By Sir WitttAM CeciL DAMPIER. 
Company, New York, 1944. 
When Sir William Cecil Dampier’s larger book, A History 
of Science and Its Relations with Philosophy and Religion, 
came from the press, we recognized that a great service had 


The 


Macmillan 


been rendered: science had been given its “‘setting in other 
modes of thought.’”” Could Dampier streamline the larger 
book and produce A Shorter History of Science without 
devitalizing the story of the growth of science which made 
the larger book famous? In the condensation, the story 
has been more simply told, but even so its vitality has been 
preserved. The connections of science with other human 
activities have been reduced to the more direct impacts, 
but these fewer situations have been used expertly to show 
the relations of science to philosophy and religion. Un- 
doubtedly the shorter history, like its monumental pro- 
totype, will become known, not simply for the facts and 
dates which it contains, but especially for the story it 
tells of the expanding mind and spirit of man. 

“Science has two streams,” we read in the first chapter, 
“corresponding to two sources, the first a gradual invention 
of tools and implements whereby men earn their living 
more safely and easily, and second the beliefs they form to 
explain the wonderful universe around them. The first may 
perhaps better be called technology, for its problems are 
too difficult for early theory, and only in later stages does 
it become applied science; the second, which in historic 
times grew into a pure search for knowledge, is the main 
subject of this book.” 

As one reads the book’s eleven chapters 
such as: the origins, Greece and Rome, 
Newton, the eighteenth century, recent biology, the new 
physics and chemistry, and the stellar universe—one turns 
back to the Preface for this summary: ‘‘The Greek 
Atomists, besides speculating on the structure of matter, 


titles 
and 


with 
Galileo 


developed therefrom a mechanical theory of life. Con- 
versely, the philosophy of Plato, as modified by Aristotle, 
laid too much stress on innate ideas and logical deduction 
to make a favorable background for the beginnings of 
inductive experimental science. To Newton and his im- 
mediate followers the Heavens declared the Glory of God, 
but Newton’s work produced a very different effect on the 
minds of Voltaire and other eighteenth-century sceptics. 
Darwin's revival of the old theory of evolution on the new 
basis of natural selection not only suggested an alternative 
origin for mankind, but spread evolutionary doctrine far 
beyond the limits of biology. The recent revolution in 
-physics has shaken the evidence for philosophic deter- 
minism which the older synthesis seemed to require. Such 
broad effects must be noticed, but the more technical 
aspects of philosophy may be passed by.” 

Exponents of general education, who insist that all 
courses should be integrated into the whole life pattern 
of the learner, will like the spirit of A Shorter History of 
Science. Its author thinks that “early specialization is 
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dangerous,”” and writes, ““Those older schoolboys whose 
chief subjects are scientific should look at them also from 
a humanist standpoint, and realize their setting in other 
modes of thought, while those studying literature need 
some knowledge of science before they can be said to be 
well educated. For both groups, I believe that the history 
of science, the story of man's attempts to understand the 
mysterious world in which he finds himself, makes the best 
way of approach to common ground.” 

Cart F. EyrinG 

Brigham Young University, Provo, Utah 


Applied Atomic Power 


By E. Smiru, A. Fox, R. Sawyer, AND H. AusTIN. 
Pp. 227, 649} in. Prentice-Hall, New York, 1946. 
Price $4.00. 


This book with a rather ambitious title is easy and 
interesting reading uncomplicated by mathematics. Under- 
graduates will find everything simple except perhaps the 
descriptions of nuclear physics work at Los Alamos in 
Appendix I, which had best be read after reading Part IIT. 

Part I contains a description of 32 uranium-bearing 
minerals and a history of radioactivity. Part II is intro- 
ductory physics, mostly off-the-slide-rule stuff. Part III 
is an abstract of the Smyth Report with two noteworthy 
portraits of Fermi and Dunning, and a description of 
the three production piles at Hanford. The description of 
the Clinton Engineer Works is more sketchy because of 
security reasons. 

Part IV, comprising the meat of the book, describes in 
considerable detail possible methods of using atomic power. 
The author, Sawyer, projects a 35,000 hp AK gas turbine 
using a closed cycle with a helium-carbon dioxide mixture. 
The inert helium makes possible a simple uranium heater. 
Many arrangements are discussed using fission heat in 
place of one fuel or another. It is emphasized that atomic 
energy will never replace other fuels but will supplement 
them. For some uses it is suggested that water be used as 
a moderator-coolant which is vaporized and then used to 
run a turbine. Atomic energy locomotives, ships, and 
rockets are considered. 

Industrial benefits from gaseous-diffusion research are 
the subject of Part V. 

Appendix III lists nuclear and atomic masses of isotopes 
from hydrogen to gold-197. A seven-page index is provided. 

DUANE STUDLEY 
Colorado Springs, Colorado 


Analytical Experimental Physics (Revised Edi- 
tion) 
By Harvey B. LEMON AND MICHAEL FERENCE, JR., 
University of Chicago. Pp. 588+XVI, Figs. 584, 
Plates 66, 2430 cm. University of Chicago Press, 
1946. Price $8.00. 

This college textbook first published in 1942 was hailed 
as the “Colossus of General Physics, one of the seven 
wonders of the textbook world.’’ Now it has been revised, 
corrections have been made, rough spots have been 
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smoothed out, figures, plates, and pages have been added. 
It is indeed a colossus in quantity and quality. 

Little revision has been made in the chapters on me- 
chanics, while in the section on heat qualitative discussions 
of the quantum theory of specific heats and of the proper- 
ties of dilute solutions have been added. The chapter on 
alternating currents now includes a section on parallel 
circuits containing resistance, inductance, and capacitance. 
Much has been added to the chapter on conduction through 
gases so that it now includes a discussion of the motion of 
charged particles in electric and magnetic fields, of rela- 
tivistic variation of mass with velocity, of the cyclotron, 
and of the betatron. As is to be expected the sections on 
nucleonics and electronics have been greatly expanded. 
A discussion is given of some nuclear reactions of tracers, 
nuclear fission, transuranic elements, nuclear chain re- 
actions, and the atomic bomb. The electronic chapter now 
includes an analysis of the equivalent circuit of a triode, 
a discussion of gas-filled tubes, pentodes, photoelectric 
cells, electronic oscillators, radio transmitting and receiving 
apparatus, microelectric waves with radar transmitting 
and receiving systems, television, and the electron micro- 
scope. An interesting plate is added on radar, showing 
typical sending and receiving apparatus, a “P.P.I. scope” 
presentation of a convoy leaving New York harbor, and 
the trace in an oscilloscope of the reflection of radar waves 
from the moon. Some additions have been made to the 
chapter on wave motion and only minor changes in 
the chapters on sound and light. A useful set of trigo- 
nometric and logarithmic tables, together with a table of 
important physical constants, has been appended. 

The book is what its title implies, both analytical and 
experimental in its approach to physics. Physical concepts 
are carefully interwoven with experimental data. The 
mathematics employed in the analysis is largely algebra 
sparingly. The 
presentation is excellent and of the kind to appeal to the 


and trigonometry with calculus used 
abler students. Though the book is more inclusive and 
more difficult than almost any other college one-year 
general text in physics, the authors have starred about one 
hundred sections which may be omitted so as to bring it 
within the scope of most college courses. There is, however, 
a section on the equivalent circuit of a triode which would 
appear to be beyond even the abler students with their 
present knowledge of. physics. The topics are discussed in 
such a manner as to leave the student with a clear idea 
of the fundamental concepts and how these were de- 
veloped. 

Since the book is so large and heavy it might have been 
more practical to publish it in two smaller volumes. 
Whether this be good or bad advice the authors are to be 
congratulated on the challenge which they have given to 
students. Any student who masters this material can be 
sure of having a solid: foundation in physics. 

R. J. STEPHENSON 
The College of Wooster 
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New Booklets 


Allied Radio Corporation, 833 West Jackson Boulevard, 
Chicago 7, Illinois, has announced the publication of a new 
164-page catalog for 1947, covering radio and electronic 
products. Special emphasis is placed on equipment for 
industrial maintenance, research, and production require- 
ments as well as for the needs of government agencies. 
10,000 items. Free on request. 








The Superpressure Division of the American Instrument 
Company, Silver Spring, Maryland, has issued its new 
88-page catalog No. 406, entitled Superpressure Catalytic 
Hydrogenation Apparatus. Inquire on company letterhead. 


United States Testing Company, Inc., 1415 Park Ave- 
nue, Hoboken, New Jersey, has revised its booklet Testing 
Terms for Better Understanding. The book describes some 
of the important tests of merchandise—tells how tests 
are made, what the results mean, what are standard and 
minimum specifications. 


The Scientific Personnel Branch of the Navy Depart- 
ment, Washington 25, D. C., issues bi-monthly a Scientific 
Personnel Bulletin, listing employment opportunities in the 
Navy Department for scientists and technicians. 76 pages 
in the January-February, 1947, issue. 


Eberbach and Son Company, Ann Arbor, Michigan, 
featured in its March issue of The Announcer of Scientific 
Equipment an article on the background and development 
of the modern analytical balance. 16 pages. Available on 
request. 


Fischer and Porter Company, Hatboro, Pennsylvania, 
has released a new bulletin including dimensional and 
operational details of an improved continuous liquid 
specific gravity indicator. Address Department 6F-B of 
the company. 


National Carbon Company, Inc., 30 East 42 Street, 
New York 17, New York, has published a new bulletin, 
Catalogue Section M8808, describing standard, seven-tube 
heat exchangers of “Karbate” impervious graphite for 
use under highly corrosive conditions. 


The Winter 1946-47 number of Interchemical Review, 
published at 432 West 45th Street, New York 19, New 
York, has the following table of contents: 

Maleic and Fumaric Resins. ...Charles S. Rowland. 


Market Research in Industrial Fields... .] John P. Duane, 
Oilcloth Pioneers in Maine. ...Charles R. Bragdon. 


The Gaertner Scientific Corporation, 1201 Wrightwood 
Avenue, Chicago 14, Illinois, has recently published a 
new bulletin, 156-74, describing ‘its optical benches and 
accessories. 26 pages. 
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Here and There 


New Appointments 








Eugene W. Boehne and Lan Jen Chu were recently 


appointed associate professors in the Department of 
Electrical Engineering at the Massachusetts Institute of 


Technology. 


[wo new appointments to the staff of the National 
Bureau of Standards were recently announced. Albert S. 
Cahn, Jr., will assist in the research involved in the design 
of two high speed electronic computing machines to be 
built for the Bureau of the Census and the Navy. P. J. 
Selgin will work on the development of electronic ordnance 
for the military services in the Ordnance Development 
Division. 


The U. S. Weather Bureau has announced the appoint- 
ment of Ross Gunn as the director of a newly organized 
Division of Physical Research, which will undertake studies 
of the atmosphere and basic physical processes that are 
related to the weather and the earth. 


F. C. Henriques, Jr., has joined the staff of Tracerlab, 
Inc., at Boston, Massachusetts, as director of the Radio- 
chemical Division. 


Awards 


Latest Franklin Institute awards for 1947 are as follows: 

The Franklin Medal, highest honor of the Institute, to 
Enrico Fermi, physicist at the Nuclear Research Institute 
in Chicago, and to Sir Robert Robinson, professor of 
chemistry at Oxford University, England. 

The Newcomen Medal to Everett G. Ackart, recently 
retired chief engineer of the E. I. 
and Company. 


du Pont de Nemours 


The Clark Medal to Edward G. Boyer, manager of the 
Philadelphia Electric Company’s gas department. 

The Ballantine Medal to George C. Southworth of the 
Bell Telephone Laboratories. 


New Chairman for A. I. P. 


Harrison, dean of the School of Science, 
Massachusetts Institute of Technology, has been elected 
chairman of the American Institute of Physics, succeeding 
Paul E. Klopsteg, director of Northwestern 
University Technological Institute, who has served as 
chairman since 1940. 


George R. 


research, 


Western Branch for Airborne Instruments Laboratory 


; A branch of the Airborne Instruments Laboratory, Inc., 
Mineola, New York, was recently established at Burbank, 
California. This will make possible a closer liaison between 
the main laboratory and west coast aircraft manufac- 


turers in need of a consultation and development service. 
Airborne Instruments Laboratory is the peacetime out- 
growth of the NDRC laboratory of the same name which 
operated under the direction of Columbia University during 
the war. Hector R. Skifter is president of the Laboratory 
and Robert D. Martin is supervising the west coast unit. 


Electrochemical Symposium at Pittsburgh 


Che annual spring symposium of the Pittsburgh section 
of the Electrochemical Society was held May 23 and 24 
at the Mellon Institute. The subject was “High Tempera- 
ture and Hetrogeneous Reactions.”’ 


Electrochemical Congress at Louisville 


A congress sponsored by the Electrochemical Society 


met in Louisville, Kentucky, April 9 to 12. A day’s session 
was devoted to electronics, one to the corrosion of metals 
at elevated temperatures, and a third to the production 
of the best steels in the electric furnace. The Young 
Authors’ Prizes were awarded to N. A. Neilsen of Du Pont 
and to Burke Cartwright of the U. S. Bureau of Mines. 


Engineering Foundation Issues Report 


The Engineering Foundation, 29 West 39th Street, New 
York 18, New York, made public in March its 32nd annual 
report to the board of trustees. The report stated, “It is 
clear that industry has, in general, become more research 
minded and large expansions of industrial research are to 
be expected. In some cases this will be by individual com- 
panies, in others, by associations of various types.’’ Six- 
teen projects were sponsored by the Engineering Founda- 
tion during the past year, six of which were carried over 
from the preceding year. Among the projects were soil 
mechanics, fluid mechanics, materials of engineering com- 
prising metal researches, welding, lubrication, and rolling 
friction. In addition to sponsoring engineering research, 
the Foundation also supports agencies or activities having 
as their objectives the advancement of the engineering 
profession. 


Changes in Electrical and Photometric Units 


In pursuance of decisions of the International Com- 
mittee on Weights and Measures, the National Bureau of 
Standards will imtroduce as of January 1, 1948, revised 
values of the units of electricity and of light. While the 
definitions of the units and the methods of fixing their 
magnitudes will be different from the present practical 
systems, the changes in magnitude will be so small as to 
affect appreciably only measurements of high precision. 
In certificates for standards and instruments issued by the 
Bureau during 1947, values will be given in both the old 
and the new units. 

A complete description of the new units has appeared in 
the Review of Scientific Instruments for May, 1947. 


JOURNAL OF APPLIED PHYSICS 





